[Libmesh-users] Reading Mesh produced from Gmsh to use as as Unstructured Mesh--Higher Order Elements?

[Zack V. <jan...@gm...>]

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