Hello,
I have a question about how ParaView handles slicing of 3D datasets and whether it is possible to preserve the original surface mesh topology. I have
1. a 2D surface mesh saved as a .vtu file
This file contains the exact triangulation of a certain surface (a boundary) of a 3D mesh.
It is a clean unstructured grid with known point positions and connectivity.
2. a 3D dataset stored in an ADIOS2 BP folder
This represents a solution on the same 3D volume mesh from which the 2D surface was originally extracted.
When loading the .bp folder, I use the following pipeline:
ADIOS2VTXReader → MergeBlocks → Slice (Plane) → CleanToGrid
The slice plane is chosen such that it should intersect the 3D mesh exactly on the same surface for which I already have the 2D reference VTU.
Visually, the produced slice looks identical to the reference .vtu surface.
However, the number of points and cells in the sliced surface is different.
So the topology of the sliced surface does not match the reference VTU surface, even though both represent the same geometric surface.
This is a problem because I want to compute errors between them which fails due to different dimensions (number of points in each file).
Here is a visual representation of the pipeline and the two vtu files (one is shifted at the moment)
Attached is the state file, the reference surface (ref_solution_4.250.vtu) and the content of the solution.dp folder (md.idx, md.0, and data.0))
ref_solution_4.250.vtu (26.1 KB)
md.idx (1.1 KB)
md.0 (18.5 KB)
data.0 (2.8 MB)
state.pvsm (634.2 KB)