Converting Gmsh node ordering in higher-order hex elements to VTK node ordering
I'm trying to visualize higher-order elements created by Gmsh in Paraview using VTK's arbritrary-order Lagrange elements. But I'm running into trouble juggling the different node orderings. I have Gmsh exporting an ASCII VTK file with what I presume is Gmsh's local node ordering. VTK expects a different local ordering of the nodes, so I have to bridge the two. On the VTK side, I have a method PointIndexFromIJK
from vtkHigherOrderHexahedron
class that allows me to build a map from the IJK coordinates of a node to its corresponding index in the VTK local ordering. Is there a similar method on the Gmsh side---maybe one I missed when I was looking through the MHexahedronN
class---that will allow me to build a similar map between IJK coordinates and the index in Gmsh local ordering?
In short, is there anything in Gmsh that could help me unscramble my elements?