dolfinx icon indicating copy to clipboard operation
dolfinx copied to clipboard

Remove `MPI_All_foo` calls in `VTKHDF5.h`

Open garth-wells opened this issue 8 months ago • 1 comments

We worked hard to remove all unnecessary MPI_All_foo function call in DOLFINx. We can't let them slip back in.

garth-wells avatar Apr 12 '25 13:04 garth-wells

We would need to figure out how to distribute the various cell types in a consistent order across all processes. The main issue is that in VTKHDF5, there is no global collection of the topology types in the mesh. (It is similar to what in XDMFFile would call a mixed cell type).

For arbitrary order Lagrange meshes, the cmaps are determined by cell_types[i], cell_offsets[i+1]-cell_offsets[i] which is read in as distributed arrays across all processes.

jorgensd avatar Apr 12 '25 17:04 jorgensd