Saving and loading dolfinx function in parallel

[dokken]

ADIOS2 shouldn’t be hard for IT staff at a super-computer facility to install.
Could you be more specific with what particular issues they had?

It is for instance on spack (Spack Packages) and similarly adios4dolfinx: Spack Packages

@MaartenV

feel like it could be possible with XDMF, but this does not seem to be recommended.

XDMF only stores data at mesh nodes, limiting what function spaces one can read and write data with.

@MaartenV

Alternatively, I could just write the end solution to a numpy array using np.save, but then I have problems due to inconsistent mesh partitioning.

I don’t know what you imply by “inconsistent” mesh partitioning. There are several posts on the forum on how to gather data to a single process, see for instance:

How to collect the global mesh without writing a file

Here is an example code of how to do it: # Copyright (C) 2023 Jørgen S. Dokken # # SPDX-License-Identifier: MIT # # Code to plot mesh with pyvista in parallel by gather data on a root process # # Last changed: 2023-10-23 import dolfinx from mpi4py import MPI import pyvista import numpy as np pyvista.start_xvfb(0.2) mesh = dolfinx.mesh.create_unit_square(MPI.COMM_WORLD, 10, 10) V = dolfinx.fem.functionspace(mesh, ("Lagrange", 2)) u = dolfinx.fem.Function(V) u.interpolate(lambda x: x[0] + 2*…

if you want to store everything as pure arrays (you have to take into account what function space you are using to consider this as a viable approach or not).