On Mar 27, 2011, at 1:40 PM, Praveen C wrote: > 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 ?
If you have the displacement field, then MappingQEulerian does exactly what you are asking for. What is currently missing from the deal.II library is a function that reads the high order element, generates a standard mesh from it, and an additional displacement field coherent with the read mesh. If you want to give it a try, just take a look at grid_in.cc. It should be self explanatory... Luca. _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
