I have a very simple PyGmsh script to produce, for example, a triangulation of the disk (attached)
I'd like to use this in "Introduduction: Example 3" (solving 2D poisson) by replacing MeshTools::Generation::build_square (mesh, n_elements_1dim, n_elements_1dim, xmin_mesh, xmax_mesh, ymin_mesh, ymax_mesh, QUAD9); with 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. Is this possible through libmesh, or something I would I have to do something in gmsh instead? Thanks for reading, and Happy Thanksgiving ------------------------------------------------------------------------------ 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