Chenchen, We only support to read in meshes that are not refined. Translated to your situation this means that you are only able to use Q1 meshes as initial grid. Of course, you can then use triangulation.set_manifold to define the correct behavior when you refine your mesh within dealii along with a Mapping of appropriate order.
Best, Daniel Am Mittwoch, 21. September 2016 17:15:05 UTC+2 schrieb Chenchen Liu: > > Hi all, > > > I have to use quadratic shape function (Q2 elements) to do my project > (2D). Previously the code works well with Q1 linear elements, when I want > to input the gird with Q2 elements generated by Gmsh, it gives the > following error: > > An error occurred in line <1418> of file <../source/grid/grid_in.cc> in > function > > void dealii::GridIn<2, 2>::read_msh(std::istream &) [dim = 2, spacedim > = 2] > > The violated condition was: > > false > > The name and call sequence of the exception was: > > ExcGmshUnsupportedGeometry(cell_type) > > Additional Information: > > The Element Identifier <8> is not supported in the deal.II library when > reading meshes in 2 dimensions. > > > And, it says the supported quad element should have 4 nodes and 4 edges. > Just want to check whether deal.ii support Q2 elements by Gmsh. > > My mesh looks like the following > > > <https://lh3.googleusercontent.com/-Lk96g0ZRgFg/V-KjPDZp3cI/AAAAAAAAAKs/NZpvWhx06Y4uL6x-6b1gLVfUAayRY5p2gCLcB/s1600/11.png> > > > > If deal.ii does not support higher order element by Gmsh, what other > format is recommended? Thanks! > > > Best, > > Chenchen > > > > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
