From 6d69eb58b6a4436054e50e15075c8859627c2cee Mon Sep 17 00:00:00 2001
From: artem-ogre <artem.ogre@gmail.com>
Date: Fri, 24 Feb 2023 15:22:23 +0100
Subject: [PATCH] #123 Make it possible to add non-boundary edges by inserting
 such edges twice into the triangulation

---
 CDT/include/Triangulation.hpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/CDT/include/Triangulation.hpp b/CDT/include/Triangulation.hpp
index ee99ec30..0c96fc7e 100644
--- a/CDT/include/Triangulation.hpp
+++ b/CDT/include/Triangulation.hpp
@@ -1605,7 +1605,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)