Deal.II cannot use triangular grids. You can get a geo file here
http://www.dolfyn.net/dolfyn/gmsh/tutorial08.html which generates a quad mesh. praveen On Fri, Mar 11, 2011 at 8:57 PM, Marek Capek <[email protected]> wrote: > Hello, > i want to make tests on benchmark Flow past the cylinder. > In application Gmsh I generated .msh file for the problem, I generated also > some refinement, but the grid is quite coarse. > I used GridIn:readMsh, but the fem program throws following exception > > ---------------------------------------------------- > Exception on processing: > -------------------------------------------------------- > An error occurred in line <642> of file > </home/marek/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::allocator<dealii::Point<spacedim> > > >&, std::vector<dealii::CellData<dim>, > std::allocator<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)--og! > ---------------------------------------------------- > > > I have thought, that Deal.II fully supports .msh format. Could anyone give > me an advice how to preprocess grids in gmsh in order to have correct > numbering. > Alternatively, anyone having mesh in deal.ii compatible format for Flow > past the cylinder problem will be welcomed :-) > Could You give me an hint on generating non-standard grids (=not from > GridGenerator), which would be fully compatible with deal's GridIn:read_xxx? > > Thank You in advance > > Marek C > > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > >
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
