pvtu with polyhedra segfault
Dear developers,
Opening an issue as suggested in discussion: https://discourse.paraview.org/t/pvtu-with-polyhedra-issue/13858
The display of polyhedra in Paraview depends on the order of pieces in the pvtu file.
In the attached example parallel_polyhedron.pvtu, there are 3 hexahedra. In the first piece parallel_polyhedron_0.vtu, there is the usual hexahedron of type 12. In the following two parallel_polyhedron_1.vtu and parallel_polyhedron_2.vtu, there are two hexes stored as polyhedra of type 42 with additional “faces” and “faceoffsets” arrays.
In the original order of pieces in the example, Paraview displays a gap in the middle of the mesh, whereas there should be three tight hexahedra. If the order of the pieces in the pvtu file is changed with parallel_polyhedron_0.vtu going last, the mesh is displayed correctly.
It appears that Paraview renders the mesh correctly if the pieces with polyhedra are followed by the pieces without polyhedra. It skips rendering the polyhedra if the piece without polyhedra is followed by the piece with polyhedra. This is verified in a large mesh example with numerous pieces.
Encountered with Paraview 5.11.1 on windows and tested with the latest ParaView-5.12.0-RC2-Windows-Python3.10-msvc2017-AMD64.msi.
The next issue with the attached example is encountered when applying a “Clip” filter. It will give a warning “A cell with a non-manifold triangulation has been encountered. This cell cannot be contoured.”. Only the top cell is rendered in this case. However, if you check “Crinkle clip”, it will show the cell in the middle (this was a gap) and wouldn’t render the cell at the bottom.
This works as described with Paraview 5.11.1 on windows. ParaView-5.12.0-RC2-Windows-Python3.10-msvc2017-AMD64.msi crashes at “Clip” filter.
Regards, Kirill