Saved XDMF file immediately crashes Paraview when opened

dolfinx
, ,
1

Hi, I’ve been using dolfinx.io.XDMFFile to visualise my solutions since I started to use dolfinx. But recently I updated dolfinx to around the commit Add test of symmetric tensor elements (#2880) · FEniCS/dolfinx@dd38c63 · GitHub And with that update, saving a Function with xdmf seems to cause paraview to crash immediately when opening. I have an MWE that saves an xdmf as follows:

from mpi4py import MPI
from dolfinx.fem import Function, functionspace
from dolfinx import io
from dolfinx import mesh

msh = mesh.create_unit_square(MPI.COMM_WORLD, 10, 10)

V = functionspace(msh, ('Lagrange', 1))
u = Function(V)
u.name = 'u'
with io.XDMFFile(MPI.COMM_WORLD, './test.xdmf', 'w') as xdmf:
    xdmf.write_mesh(msh)
    xdmf.write_function(u)

but trying to open test.xdmf immediately crashes Paraview with no error output. Here is a google drive of the output files.

For your information I’m using Paraview 5.10.0 and the script was run on CentOS 8 with a dolfinx installed via conda (following the instructions for C++ core build and pip install).

2

Maybe try to update paraview to 5.11.0 or 5.12.0? I can open the xmdf files on your google drive, as well as those produced by running again your code, with paraview 5.11.0 installed by debian apt.

As expected, those files contain a square domain and a zero solution.

3

Okay, seems to work with 5.11.0 on CentOS and 5.10.1 on Windows. Just not 5.10.0.