It, however, does not save functions in association with the mesh, so, when loading those functions, you would need to provide the exact same mesh with the same parallel distribution as that those functions were originally defined on.

A new checkpointing feature though will be available in a month or so, which facilitates saving functions in associations with meshes using arbitrary number of MPI processes. It is designed for a general DoF layout, so it will still save DoF values serialised.

From: firedrake-bounces@imperial.ac.uk <firedrake-bounces@imperial.ac.uk> on behalf of Andrew Hicks <ahick17@lsu.edu>
Sent: Tuesday, September 28, 2021 3:53 AM
To: firedrake <firedrake@imperial.ac.uk>
Subject: [firedrake] Writing to file in parallel

Hello all,

The website mentions a way to print data from a function that has been computed in parallel:

with u.dat.vec_ro as vu:

vu.view()

However, I want to write this data to a file, and I’m not sure how. Moreover, is there a way I can format this so that instead of printing everything in one row, it prints it in as many rows as there are coordinates of the function? So for instance, can I write a file with three rows of data for a three dimensional function space?

Thank you,

Andrew Hicks