From 06e630287bc4adb5e9c3f2f17aa6744d891ba26b Mon Sep 17 00:00:00 2001
From: Laurent Rineau <laurent.rineau@cgal.org>
Date: Tue, 14 Nov 2023 17:50:45 +0100
Subject: [PATCH] another fix

Now it seems to work!
---
 .../include/CGAL/Constrained_Delaunay_triangulation_3.h    | 7 ++++++-
 1 file changed, 6 insertions(+), 1 deletion(-)

diff --git a/Triangulation_3/include/CGAL/Constrained_Delaunay_triangulation_3.h b/Triangulation_3/include/CGAL/Constrained_Delaunay_triangulation_3.h
index 84eb5424dc0..6e0ec1e9f2f 100644
--- a/Triangulation_3/include/CGAL/Constrained_Delaunay_triangulation_3.h
+++ b/Triangulation_3/include/CGAL/Constrained_Delaunay_triangulation_3.h
@@ -1803,12 +1803,14 @@ private:
       insert_new_vertex(v);
     }
 
-    std::vector<Facet> missing_faces;
+    boost::container::small_vector<Facet, 32> missing_faces;
     do {
       missing_faces.clear();
+      boost::container::small_vector<Facet, 32> internal_facets;
       for(auto f : facets_of_cavity_border) {
         if(cells_of_cavity.contains(f.first)) {
           // internal facet, due to cavity growing
+          internal_facets.push_back(f);
           continue;
         }
         const auto [v0, v1, v2] = this->make_vertex_triple(f);
@@ -1821,6 +1823,9 @@ private:
           missing_faces.push_back(f);
         }
       }
+      for(auto f : internal_facets) {
+        facets_of_cavity_border.erase(f);
+      }
       for(auto [cell, facet_index] : missing_faces) {
         facets_of_cavity_border.erase({cell, facet_index});
         {