Delete visibility scalar data if it already exists

ViennaTools · Nov 22, 2024 · 8200717 · 8200717
1 parent 476b274
commit 8200717
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/include/viennals/lsCalculateVisibilities.hpp b/include/viennals/lsCalculateVisibilities.hpp
@@ -88,8 +88,12 @@ template <class NumericType, int D> class CalculateVisibilities {
       }
     }
 
-    levelSet->getPointData().insertNextScalarData(visibilities,
-                                                  visibilitiesLabel);
+    auto &pointData = levelSet->getPointData();
+    // delete if already exists
+    if (int i = pointData.getScalarDataIndex(visibilitiesLabel); i != -1) {
+      pointData.eraseScalarData(i);
+    }
+    pointData.insertNextScalarData(visibilities, visibilitiesLabel);
 
     assert(id == numDefinedPoints);
   }