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

Reply via email to