Save a Mesh in parallel.

Question

I have a simulation with an evolving mesh, and would like to save the resulting mesh in parallel. Is this possible?

I can save the mesh as a .pvd file, but when I load it,

mesh = Mesh('mesh.pvd')

I get the error

*** -------------------------------------------------------------------------
*** DOLFIN encountered an error. If you are not able to resolve this issue
*** using the information listed below, you can ask for help at
***
***     https://answers.launchpad.net/dolfin
***
*** Remember to include the error message listed below and, if possible,
*** include a *minimal* running example to reproduce the error.
***
*** -------------------------------------------------------------------------
*** Error:   Unable to read object from file.
*** Reason:  Cannot read objects of type Mesh from VTK files.
*** Where:   This error was encountered inside GenericFile.cpp.
*** Process: 0
*** -------------------------------------------------------------------------

Answer 1

Hi, for parallel IO with mesh, you should use the XDMF file format

File('mesh.xdmf') << mesh

Also take a look at this previously asked question.