<< "3 Quadrilateral (4 nodes, 4 edges).\n"
<< "5 Hexahedron (8 nodes, 12 edges, 6 faces).\n"
<< "15 Point (1 node, ignored when read)");
+
+
+ DeclException0 (ExcGmshNoCellInformation);
protected:
/**
* Store address of the triangulation to
in >> version >> file_type >> data_size;
- Assert (version == 2.0, ExcNotImplemented());
+ Assert ( (version >= 2.0) &&
+ (version <= 2.1), ExcNotImplemented());
Assert (file_type == 0, ExcNotImplemented());
Assert (data_size == sizeof(double), ExcNotImplemented());
ExcInvalidGMSHInput(line));
in >> n_cells;
+
// set up array of cells
std::vector<CellData<dim> > cells;
SubCellData subcelldata;
AssertThrow (in, ExcIO());
- // do some clean-up on vertices...
+ // check that we actually read some
+ // cells.
+ AssertThrow(cells.size() > 0, ExcGmshNoCellInformation());
+
+ // do some clean-up on
+ // vertices...
GridTools::delete_unused_vertices (vertices, cells, subcelldata);
// ... and cells
if(dim==spacedim)