Hello,

Is it possible to build/use a mesh consisting of 20-node hexes in PETSc? 

I tried reading (-dm_plex_filename) an existing mesh, with a 20-node cube, from gmsh (using the option "SecondOrderIncomplete"); statistics show 20 nodes, however the "Entities" section still only contains 8 vertices, so that may not work here.

Using -petscdualspace_lagrange_XYZ I could not achieve this either.

I searched some old threads in the mailing list; there are a couple, but they don't specify how they build it; they just mention they are using such mesh.

Searching for "serendipity" in the documentation did not return anything either. 

Thanks,