Quad mesh nodes reordered upon reading XDMF file

Most software reorder meshes (the topology) (usually to be able to use higher order Nedelec/DG spaces, see for instance the introduction of: [2102.11901] Construction of arbitrary order finite element degree-of-freedom maps on polygonal and polyhedral cell meshes). However, in DOLFINx, we use dof-transformations (see said paper, which in theory means that you do not need to reorder the mesh topology).

The re-ordering you see in the geometry, is to increase data-locality, see: https://github.com/FEniCS/dolfinx/blob/23c76248f0686a795f170927036f1331416c1876/cpp/dolfinx/mesh/Geometry.cpp#L73-L105
@garth or @chris can probably provide you with more information regarding this.

It would be helpful to know what kind of software you would like to interface with, and how you were thinking about making this interface (if we do not consider the reordering of the geometry).

1 Like