Thanks Roy and Alexander, I appreciate the insight, I will ask those maintaining PyGmsh about the order of elements if I can't seem to find the option readily myself, and then report back here.
Disregarding order for now, if I simply run this as it is, is it assuming Tri3 elements and a number of elements corresponding to those given in the mesh? On Wed, Nov 22, 2017 at 4:36 PM, Alexander Lindsay <alexlindsay...@gmail.com > wrote: > I've never used pygmsh, but I know if you're running gmsh from the command > line, passing the option `-order 2` will create second order elements. > > On Wed, Nov 22, 2017 at 2:14 PM, Roy Stogner <royst...@ices.utexas.edu> > wrote: > >> >> On Wed, 22 Nov 2017, Zack Vitoh wrote: >> >> GmshIO gmsh_io(mesh); >>> gmsh_io.read("disk.vtu"); >>> >>> but--and I am not sure if this is possible--imposing additional nodes on >>> each element generated in the triangularization generated by gmsh with >>> nodes such that I have, say, Tri6 elements. >>> >> >> You have Tri3 elements now? After the read calling >> mesh.all_second_order(); >> will turn them into Tri6 elements. >> >> Is this possible through libmesh, or something I would I have to do >>> something in gmsh instead? >>> >> >> I'm afraid if you have a disk then this is something you *should* do >> in gmsh instead if that's possible. By the time libMesh reads a >> first-order triangulation, curvature information has been discarded, >> and we don't try to reconstruct it, so you'll still get straight-sided >> elements rather than curved-sided elements. The only way to get the >> best second-order boundary geometry is to patch it up manually, >> looping over boundary nodes and "snapping" them to the curve where you >> know they should be. >> --- >> Roy >> >> >> ------------------------------------------------------------ >> ------------------ >> Check out the vibrant tech community on one of the world's most >> engaging tech sites, Slashdot.org! http://sdm.link/slashdot >> _______________________________________________ >> Libmesh-users mailing list >> Libmesh-users@lists.sourceforge.net >> https://lists.sourceforge.net/lists/listinfo/libmesh-users >> > > ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users