I'm trying import into deal.ii a mesh generated with Gmsh. The only
thing I know is that deal.ii doesn't accept triangular elements from
imported files.

I construct my grid in Gmsh and at the end I use the command 'Recombine
Surface' to recombine the triangular elements into quads. I save the
file in Gmsh V1.0 and I make sure that all the elements were quads, i.e.
all element types were 3 (4 node quadrangle).

Then, when I tried to import the mesh into deal.ii and it gave the
error:

terminate called after throwing an instance of
'dealii::StandardExceptions::ExcInternalError'
  what():  --------------------------------------------------------
An error occurred in line <642> of file
</home/zerpiko/FEM/deal.II/source/grid/grid_reordering.cc> in function
static void dealii::GridReordering<dim,
spacedim>::invert_all_cells_of_negative_grid(const
std::vector<dealii::Point<spacedim> >&,
std::vector<dealii::CellData<dim> >&) [with int dim = 2, int spacedim =
2]
The violated condition was: 
    n_negative_cells==0 || n_negative_cells==cells.size()
The name and call sequence of the exception was:
    ExcInternalError()
Additional Information: 
(none)
--------------------------------------------------------

I think is something related with reordering cells. There are any other
restrictions that a mesh must fulfill? Tnks.


_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to