#123 Make it possible to add non-boundary edges by inserting such edg…

…es twice into the triangulation
artem-ogre · Apr 13, 2023 · 212e3fe · 212e3fe
1 parent 7e2432f
commit 212e3fe
Show file tree

Hide file tree

Showing 2 changed files with 30 additions and 1 deletion.
diff --git a/CDT/include/Triangulation.hpp b/CDT/include/Triangulation.hpp
@@ -1717,7 +1717,7 @@ Triangulation<T, TNearPointLocator>::peelLayer(
     {
         const TriInd iT = seeds.top();
         seeds.pop();
-        triDepths[iT] = layerDepth;
+        triDepths[iT] = std::min(triDepths[iT], layerDepth);
         behindBoundary.erase(iT);
         const Triangle& t = triangles[iT];
         for(Index i(0); i < Index(3); ++i)

diff --git a/visualizer/data/issue-42-full-boundary-overlap-2.txt b/visualizer/data/issue-42-full-boundary-overlap-2.txt
@@ -0,0 +1,29 @@
+12 16
+0 0
+1 0
+1 1
+0 1
+0.25 0.25
+0.75 0.25
+0.75 0.75
+0.25 0.75
+0.4 0.4
+0.6 0.4
+0.6 0.6
+0.4 0.6
+0 1
+1 2
+2 3
+3 0
+4 5
+4 5
+5 6
+5 6
+6 7
+6 7
+7 4
+7 4
+8 9
+9 10
+10 11
+11 8