I/O from XDMF/HDF5 files in dolfin-x

If you want to

  1. Read a mesh_0 from file
  2. Store a solution to file (f). let us call this mesh mesh_1.
  3. Read data from another file (based on the mesh_0) and use with f from mesh_1

You would have to store more information than I currently do. You would have to store the original cell index and the original geometry indices.

It is not infeasible, but it would slow down storage of checkpoints (and reading them) a bit, as one of the operations would have to deal with this reordering (reading a global array on every process).

The workflow I would suggest is shown in: