XDMF+HDF5 file crashes ParaView non-deterministically

ParaView Support
1

Hi everyone,

My application writes the results into an HDF5 file, and also creates an XDMF file to interpret the data. When I open the XDMF file in Paraview 5.13.2 and click on Apply, ParaView crashes and exits.

Here are my HDF5 and XDMF files:

circle.xdmf (3.7 KB)
circle.hdf5 (2.1 MB)

This is a point cloud with two XDMF attributes attached.

Observations, which may help to find the problem:

  1. The HDF5 file should be correct: I read it with h5py and visualized it with Matplotlib. It behaves as intended.
  2. The same crash occurs with ParaView 5.11 too.
  3. Whatever XDMF reader I use (there are three in ParaView), the same crash happens.
  4. If you delete the two attributes (“neighbors” and “flags”) from the XDMF file, it still crashes, so I believe that the problem lies in the topology or the geometry.
  5. If I delete the two last grids (“Level_2” and “Level_3”) from the XDMF file, I can open it in ParaView
    1. You need to select the Point Gaussian representation; I don’t understand why the Point representation does not show all the points.
    2. If you insert the last two grids back into the XDMF file and reload the data in ParaView (F5), the whole grid is show. However, after a few times you change the Gaussian Radius property, ParaView crashes again.

Do you have any idea about the nature of the problem?

Thank you.

2

I found the problem: the point numbering must be contiguous and start from 0 because these indices are used in the XDMF topology of type “Polyvertex” to index the point coordinates in the XDMF geometry.