Mesh I/O could be broken in serial
Issue #9
new
I accidentally discovered this behavior.
Consider for example Unicorn's smCube
benchmark. It runs fine in serial an in parallel, and mesh conversion done with dolfin_post
runs without errors in both cases with command
./dolfin_post -m mesh_out.bin -t vtk -n 1 -s velocity
However, when trying to open the generated pvd
files in paraview
, the one generated from conversion of the mesh and solution output in parallel is displayed as expected, while the one generated from output from the serial run yields the following error in paraview
#!
ERROR: In /build/paraview-arsa8T/paraview-5.0.1+dfsg1/VTK/IO/XML/vtkXMLDataReader.cxx, line 449
vtkXMLUnstructuredGridReader (0x426c9a0): Cannot read point data array "U" from PointData in piece 0. The data array in the element may be too short.
It is not clear if the error is in dolfin-hpc
or in dolfin_post
.