ccolin/pfe
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

remove dupes from surfacePoints list

Browse Source
This commit is contained in:
papush! 2022-03-24 13:48:18 +01:00
parent 1782687e9a
commit 226ef2dc0a
1 changed files with 10 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
src/surface_points_filter.cc
View File

@ -10,6 +10,8 @@
#include <vtkStaticCellLinks.h> #include <vtkStaticCellLinks.h>
#include <vtkIntArray.h> #include <vtkIntArray.h>
#include <set>
vtkStandardNewMacro(SurfacePointsFilter); vtkStandardNewMacro(SurfacePointsFilter);
@ -46,6 +48,7 @@ vtkTypeBool SurfacePointsFilter::RequestData(
vtkNew<vtkIdList> neighborCells; vtkNew<vtkIdList> neighborCells;
vtkNew<vtkIdList> facePoints; vtkNew<vtkIdList> facePoints;
std::set<vtkIdType> surfacePointsSet;
auto *it = input->NewCellIterator(); auto *it = input->NewCellIterator();
for (it->InitTraversal(); !it->IsDoneWithTraversal(); it->GoToNextCell()) { for (it->InitTraversal(); !it->IsDoneWithTraversal(); it->GoToNextCell()) {
vtkIdList *cellPointIds = it->GetPointIds(); vtkIdList *cellPointIds = it->GetPointIds();
@ -60,15 +63,19 @@ vtkTypeBool SurfacePointsFilter::RequestData(
neighborCells); neighborCells);
facePoints->Reset(); facePoints->Reset();
if (neighborCells->GetNumberOfIds() == 0) { if (neighborCells->GetNumberOfIds() == 0) {
surfacePoints->InsertNextValue(idA); surfacePointsSet.insert(idA);
surfacePoints->InsertNextValue(idB); surfacePointsSet.insert(idB);
surfacePoints->InsertNextValue(idC); surfacePointsSet.insert(idC);
isSurface->SetValue(idA, 1); isSurface->SetValue(idA, 1);
isSurface->SetValue(idB, 1); isSurface->SetValue(idB, 1);
isSurface->SetValue(idC, 1); isSurface->SetValue(idC, 1);
} }
} }
} }
surfacePoints->Allocate(surfacePointsSet.size());
for (const vtkIdType &id : surfacePointsSet) {
surfacePoints->InsertNextValue(id);
}
output->GetPointData()->SetScalars(isSurface); output->GetPointData()->SetScalars(isSurface);
output->GetFieldData()->AddArray(surfacePoints); output->GetFieldData()->AddArray(surfacePoints);