I have an FEA mesh of 3d cells, with results that are discontinuous between cells. I model this by splitting the points such that each cell has its own points.
I want to cut through those cells, yielding a single polygon per cell.
Initially I just used a vtkCutter and set GenerateTrianglesOff to produce polygons. This works somewhat but it merges back the points that were split, causing unexpected irregularity in the bands that are generated by the banded contour filter.
To address this I used a vtkNonMergingPointLocator for the vtkCutter, but now I end up with triangles even though I said that I don’t want triangles…
An image of the cut surface before the contouring shows clearly that triangles are generated…
I realize that what I’m asking is somewhat of a blend between non-merging and merging point-locator behaviors: merging while dealing with a single input cell but non-merging across input cells. Is there a way to achieve this with the current locators and cutter or do I have to roll my own, or should this be considered a bug of vtkCutter?
Here’s the complete script cutting-plane.py (2.8 KB)