No, you cannot enforce that, as
- The order of vertices only makes sense in serial.
- Data locality is important for performance, and therefore it is reordered to ensure this.
A question that comes to mind:
Why do you need to preserve the vertex/node order?
You can get the mapping by looking at mesh.geometry.input_global_indices,
or: How to apply a boundary condition to a range of DOFs? - #7 by dokken