> > > > > Are you saying that you are generating meshes with curved boundaries > that > > > are > > > described by quadratic functions? How are these represented in the > output > > > file? I thought gmsh format can only represent cells that are described > > > by 4 > > > (in 2d) or 8 (in 3d) vertices, which would of course know nothing about > > > curved > > > boundaries. > > > > > > Yes. gmsh can generate higher order elements, see this > > > > http://www.geuz.org/gmsh/doc/texinfo/gmsh.html#Node-ordering > > What I meant to say is this: Are you actually writing your mesh with these > higher order cells? Because I don't think that the GridIn class can *read* > these cells from .msh files... > > Yes, I tried to read the higher order elements, which of course did not work.
I need a higher order representation of the boundary, which gmsh can already provide. Would it be possible to use the higher order elements as given in gmsh to construct the mapping in deal.ii ? I am thinking of a general solution which should work for any geometry described by a gmsh mesh file, maybe using the approach suggested by Luca. I understood his idea of creating the displacement field but how to create the mapping from the displacement field ? Thanks praveen
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
