CGAL · MaelRL · Feb 17, 2026 · Feb 17, 2026 · Feb 17, 2026 · Feb 17, 2026
diff --git a/Generator/examples/Generator/random_points_in_tetrahedral_mesh_3.cpp b/Generator/examples/Generator/random_points_in_tetrahedral_mesh_3.cpp
@@ -32,7 +32,7 @@ typedef CGAL::Mesh_triangulation_3<
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr>               C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr>                                 Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3>                               Mesh_criteria;
 
 // (Unweighted) point type
 typedef Tr::Bare_point                                            Point;

diff --git a/Generator/examples/Generator/random_points_on_tetrahedral_mesh_3.cpp b/Generator/examples/Generator/random_points_on_tetrahedral_mesh_3.cpp
@@ -31,7 +31,7 @@ typedef CGAL::Mesh_triangulation_3<
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr>               C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr>                                 Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3>                               Mesh_criteria;
 
 // (Unweighted) point type
 typedef Tr::Bare_point                                            Point;

diff --git a/Generator/test/Generator/generic_random_test.cpp b/Generator/test/Generator/generic_random_test.cpp
@@ -53,7 +53,7 @@ typedef Mesh_domain::Curve_index                                  Curve_index;
 typedef CGAL::Mesh_triangulation_3<Mesh_domain>::type             Tr;
 typedef CGAL::Mesh_complex_3_in_triangulation_3<
           Tr, Corner_index, Curve_index>                          C3t3;
-typedef CGAL::Mesh_criteria_3<Tr>                                 Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3>                               Mesh_criteria;
 typedef C3t3::Point                                               Point;
 
 int test_triangles_2(const FT eps)

diff --git a/Installation/include/CGAL/Mesh_3/Mesh_complex_3_in_triangulation_3_fwd.h b/Installation/include/CGAL/Mesh_3/Mesh_complex_3_in_triangulation_3_fwd.h
@@ -24,14 +24,14 @@ template <typename Tr,
           typename CurveIndex = int>
 class Mesh_complex_3_in_triangulation_3;
 
-template<class Tr, bool c3t3_loader_failed>
-bool build_triangulation_from_file(std::istream& is,
-                                   Tr& tr,
-                                   bool replace_domain_0);
-
-template<class Tr, bool c3t3_loader_failed>
-bool build_triangulation_from_file(std::istream& is,
-                                   Tr& tr);
+template <class C3T3, bool c3t3_loader_failed>
+bool build_mesh_complex_from_file(std::istream& is,
+                                  C3T3& c3t3,
+                                  bool replace_domain_0);
+
+template <class C3T3, bool c3t3_loader_failed>
+bool build_mesh_complex_from_file(std::istream& is,
+                                  C3T3& c3t3);
 } // CGAL
 #endif
 

diff --git a/Installation/include/CGAL/SMDS_3/Mesh_complex_3_in_triangulation_3_fwd.h b/Installation/include/CGAL/SMDS_3/Mesh_complex_3_in_triangulation_3_fwd.h
@@ -37,12 +37,19 @@ namespace IO {
 
 namespace SMDS_3 {
 
-  template<class Tr>
-  bool build_triangulation_from_file(std::istream& is,
-                                     Tr& tr,
-                                     bool verbose = false,
-                                     bool replace_domain_0 = false,
-                                     bool allow_non_manifold = false);
+template<class C3T3>
+bool build_mesh_complex_from_file(std::istream& is,
+                                  C3T3& c3t3,
+                                  bool verbose = false,
+                                  bool replace_domain_0 = false,
+                                  bool allow_non_manifold = false);
+
+template<class Tr>
+bool build_triangulation_from_file(std::istream& is,
+                                   Tr& tr,
+                                   bool verbose = false,
+                                   bool replace_domain_0 = false,
+                                   bool allow_non_manifold = false);
 
 } // namespace SMDS_3
 } // namespace CGAL

diff --git a/Lab/demo/Lab/Plugins/IO/VTK_io_plugin.cpp b/Lab/demo/Lab/Plugins/IO/VTK_io_plugin.cpp
@@ -428,7 +428,7 @@ class CGAL_Lab_vtk_plugin :
         subdomains.push_back(si);
       }
       std::map<Facet, int> border_facets;
-      //Preprocessing for build_triangulation
+      //Preprocessing for build_mesh_complex
       //check for orientation and swap in cells if not good.
       for(std::size_t i=0; i<finite_cells.size(); ++i)
       {
@@ -445,8 +445,8 @@ class CGAL_Lab_vtk_plugin :
         }
       }
 
-      CGAL::SMDS_3::build_triangulation_with_subdomains_range(
-        c3t3_item->c3t3().triangulation(),
+      CGAL::SMDS_3::build_mesh_complex_with_subdomains_range(
+        c3t3_item->c3t3(),
         points, finite_cells, subdomains, border_facets,
         false, false, true);
 
@@ -459,9 +459,9 @@ class CGAL_Lab_vtk_plugin :
           c3t3_item->c3t3().add_to_complex(cit, cit->info());
         for(int i=0; i < 4; ++i)
         {
-          if(cit->surface_patch_index(i)>0)
+          if(c3t3_item->c3t3().surface_patch_index(cit, i) > 0)
           {
-            c3t3_item->c3t3().add_to_complex(cit, i, cit->surface_patch_index(i));
+            c3t3_item->c3t3().add_to_complex(cit, i, c3t3_item->c3t3().surface_patch_index(cit, i));
           }
         }
       }
@@ -480,11 +480,11 @@ class CGAL_Lab_vtk_plugin :
           Cell_handle nc = c->neighbor(fit->second);
 
           // By definition, Subdomain_index() is supposed to be the id of the exterior
-          if(c->subdomain_index() != C3t3::Triangulation::Cell::Subdomain_index() &&
-             nc->subdomain_index() == C3t3::Triangulation::Cell::Subdomain_index())
+          if(c3t3_item->c3t3().subdomain_index(c) != C3t3::Subdomain_index() &&
+             c3t3_item->c3t3().subdomain_index(nc) == C3t3::Subdomain_index())
           {
             // Color the border facet with the index of its cell
-            c3t3_item->c3t3().add_to_complex(c, fit->second, c->subdomain_index());
+            c3t3_item->c3t3().add_to_complex(c, fit->second, c3t3_item->c3t3().subdomain_index(c));
           }
         }
       }

diff --git a/Lab/demo/Lab/Plugins/Mesh_3/C3t3_io_plugin.cpp b/Lab/demo/Lab/Plugins/Mesh_3/C3t3_io_plugin.cpp
@@ -137,7 +137,7 @@ CGAL_Lab_c3t3_binary_io_plugin::load(
       item->setName(fileinfo.completeBaseName());
       item->set_valid(false);
 
-      if(CGAL::SMDS_3::build_triangulation_from_file(in, item->c3t3().triangulation(),
+      if(CGAL::SMDS_3::build_mesh_complex_from_file(in, item->c3t3(),
          /*verbose = */false, /*replace_subdomain_0 = */false, /*allow_non_manifold = */true))
       {
         update_c3t3(item->c3t3());
@@ -167,7 +167,7 @@ CGAL_Lab_c3t3_binary_io_plugin::load(
 
       item->setName(fileinfo.completeBaseName());
 
-      if (CGAL::IO::read_AVIZO_TETRA(in, item->c3t3().triangulation()))
+      if (CGAL::IO::read_AVIZO_TETRA(in, item->c3t3()))
       {
         update_c3t3(item->c3t3());
 
@@ -449,14 +449,14 @@ update_c3t3(C3t3& c3t3)
   c3t3.rescan_after_load_of_triangulation(); //fix counters for facets and cells
   for (Cell_handle cit : c3t3.triangulation().finite_cell_handles())
   {
-    CGAL_assertion(cit->subdomain_index() >= 0);
-    if (cit->subdomain_index() != C3t3::Triangulation::Cell::Subdomain_index())
-      c3t3.add_to_complex(cit, cit->subdomain_index());
+    CGAL_assertion(c3t3.subdomain_index(cit) >= 0);
+    if (c3t3.subdomain_index(cit) != C3t3::Subdomain_index())
+      c3t3.add_to_complex(cit, c3t3.subdomain_index(cit));
 
     for (int i = 0; i < 4; ++i)
     {
-      if (cit->surface_patch_index(i) > 0)
-        c3t3.add_to_complex(cit, i, cit->surface_patch_index(i));
+      if (c3t3.surface_patch_index(cit, i) > 0)
+        c3t3.add_to_complex(cit, i, c3t3.surface_patch_index(cit, i));
     }
   }
 
@@ -469,11 +469,11 @@ update_c3t3(C3t3& c3t3)
       Cell_handle nc = c->neighbor(fit.second);
 
       // By definition, Subdomain_index() is supposed to be the id of the exterior
-      if (c->subdomain_index() != C3t3::Triangulation::Cell::Subdomain_index() &&
-          nc->subdomain_index() == C3t3::Triangulation::Cell::Subdomain_index())
+      if (c3t3.subdomain_index(c) != C3t3::Subdomain_index() &&
+          c3t3.subdomain_index(nc) == C3t3::Subdomain_index())
       {
         // Color the border facet with the index of its cell
-        c3t3.add_to_complex(c, fit.second, c->subdomain_index());
+        c3t3.add_to_complex(c, fit.second, c3t3.subdomain_index(c));
       }
     }
   }

diff --git a/Lab/demo/Lab/Plugins/Mesh_3/Mesh_function.h b/Lab/demo/Lab/Plugins/Mesh_3/Mesh_function.h
@@ -118,7 +118,7 @@ class Mesh_function
   typedef C3t3::Vertex_handle                       Vertex_handle;
 
   typedef C3t3::Triangulation                       Tr;
-  typedef CGAL::Mesh_criteria_3<Tr>                 Mesh_criteria;
+  typedef CGAL::Mesh_criteria_3<C3t3>                 Mesh_criteria;
   typedef Mesh_criteria::Edge_criteria              Edge_criteria;
   typedef Mesh_criteria::Facet_criteria             Facet_criteria;
   typedef Mesh_criteria::Cell_criteria              Cell_criteria;
@@ -409,10 +409,8 @@ launch()
                        0,
                        &stop_);
 
-#ifdef CGAL_MESH_3_PROFILING
   CGAL::Real_timer t;
   t.start();
-#endif
 
 #if CGAL_MESH_3_MESHER_STATUS_ACTIVATED
   mesher_->initialize();
@@ -430,15 +428,16 @@ launch()
 
   // Ensure c3t3 is ok (useful if process has been stop by the user)
   mesher_->fix_c3t3();
-  std::cerr<<"Done."<<std::endl;
+  t.stop();
+  std::cerr<<"Done. (in " << t.time() << " seconds, wall-clock)" << std::endl;
 }
 
 
 template < typename D_, typename Tag >
 void
 Mesh_function<D_,Tag>::
 tweak_criteria(Mesh_criteria& c, Mesh_fnt::Polyhedral_domain_tag) {
-  typedef CGAL::Mesh_3::Facet_topological_criterion_with_adjacency<Tr,
+  typedef CGAL::Mesh_3::Facet_topological_criterion_with_adjacency<C3t3,
        Domain, typename Facet_criteria::Visitor> New_topo_adj_crit;
 
   if(((int(c.facet_criteria_object().topology()) &

diff --git a/Lab/demo/Lab/Plugins/Mesh_3/Offset_meshing_plugin.cpp b/Lab/demo/Lab/Plugins/Mesh_3/Offset_meshing_plugin.cpp
@@ -272,8 +272,7 @@ SMesh* cgal_off_meshing(QWidget*,
 
   using Mesh_domain_base = CGAL::Labeled_mesh_domain_3<GT, int, int>;
   using Mesh_domain = CGAL::Mesh_domain_with_polyline_features_3<Mesh_domain_base>;
-  using Tr = C3t3::Triangulation;
-  using Mesh_criteria = CGAL::Mesh_criteria_3<Tr>;
+  using Mesh_criteria = CGAL::Mesh_criteria_3<C3t3>;
 
   const CGAL::Bbox_3 bbox = source_item->bbox();
 

diff --git a/Lab/demo/Lab/Plugins/Mesh_3/config_mesh_3.h b/Lab/demo/Lab/Plugins/Mesh_3/config_mesh_3.h
@@ -54,7 +54,6 @@
 // MESH_3 GENERAL PARAMETERS
 // ==========================================================================
 
-//#define CGAL_MESH_3_USE_OLD_SURFACE_RESTRICTED_DELAUNAY_UPDATE // WARNING: VERY SLOW
 //#define CGAL_MESH_3_INITIAL_POINTS_NO_RANDOM_SHOOTING
 //#define CGAL_MESHES_DEBUG_REFINEMENT_POINTS
 //#define CHECK_AND_DISPLAY_THE_NUMBER_OF_BAD_ELEMENTS_IN_THE_END

diff --git a/Lab/demo/Lab/Plugins/Point_set/Surface_reconstruction_poisson_impl.cpp b/Lab/demo/Lab/Plugins/Point_set/Surface_reconstruction_poisson_impl.cpp
@@ -187,7 +187,7 @@ SMesh* poisson_reconstruct(Point_set& points,
   typedef CGAL::Labeled_mesh_domain_3<Kernel> Mesh_domain;
   typedef typename CGAL::Mesh_triangulation_3<Mesh_domain, CGAL::Default, Concurrency_tag>::type Tr;
   typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
-  typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+  typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
   // AABB tree
   typedef CGAL::AABB_face_graph_triangle_primitive<SMesh> Primitive;

diff --git a/Lab/demo/Lab/Scene_tetrahedra_item.cpp b/Lab/demo/Lab/Scene_tetrahedra_item.cpp
@@ -183,7 +183,7 @@ void Scene_tetrahedra_item::computeElements()const
     const Tr::Bare_point& pb = wp2p(cit->vertex(1)->point());
     const Tr::Bare_point& pc = wp2p(cit->vertex(2)->point());
     const Tr::Bare_point& pd = wp2p(cit->vertex(3)->point());
-    const QColor color = d->c3t3_item->getSubdomainIndexColor(cit->subdomain_index());
+    const QColor color = d->c3t3_item->getSubdomainIndexColor(c3t3.subdomain_index(cit));
     //0 - 1 : dihedral angle
     double min_dihedral_angle = 360.0;
     double max_dihedral_angle = -360.0;

diff --git a/Mesh_3/benchmark/Mesh_3/benchmark_config.h b/Mesh_3/benchmark/Mesh_3/benchmark_config.h
@@ -28,7 +28,6 @@
 // #define CGAL_MESH_3_POLYHEDRON_WITH_FEATURES
 // #define CGAL_MESH_3_IMPLICIT_WITH_FEATURES
 
-// #define CGAL_MESH_3_USE_OLD_SURFACE_RESTRICTED_DELAUNAY_UPDATE // WARNING: VERY SLOW
 #define CGAL_MESH_3_INITIAL_POINTS_NO_RANDOM_SHOOTING
 
 // #define CGAL_MESH_3_BENCHMARK_LLOYD

diff --git a/Mesh_3/benchmark/Mesh_3/benchmark_mesh_3.cpp b/Mesh_3/benchmark/Mesh_3/benchmark_mesh_3.cpp
@@ -148,7 +148,7 @@ enum Exit_code
 // the call to Mesh_3 happens here
 template <typename C3t3, typename Domain>
 Exit_code make_mesh(const Domain& domain,
-                    const CGAL::Mesh_criteria_3<typename C3t3::Triangulation>& criteria,
+                    const CGAL::Mesh_criteria_3<C3t3>& criteria,
                     const std::string& output_filename)
 {
 #ifdef _DEBUG
@@ -244,7 +244,7 @@ Exit_code make_mesh(const Domain& domain,
                     double cell_sizing, double cell_shape,
                     const std::string& output_filename)
 {
-  typedef CGAL::Mesh_criteria_3<typename C3t3::Triangulation> Mesh_criteria;
+  typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
   const CGAL::Bbox_3 bbox = domain.bbox();
 

diff --git a/Mesh_3/doc/Mesh_3/Concepts/MeshCriteriaWithFeatures_3.h b/Mesh_3/doc/Mesh_3/Concepts/MeshCriteriaWithFeatures_3.h
@@ -15,7 +15,7 @@ that describes the requirements, in terms of sizing, for the discretization of t
 \cgalRefines{MeshCriteria_3}
 
 \cgalHasModelsBegin
-\cgalHasModels{CGAL::Mesh_criteria_3<Tr>}
+\cgalHasModels{CGAL::Mesh_criteria_3}
 \cgalHasModelsEnd
 
 \sa `MeshEdgeCriteria_3`

diff --git a/Mesh_3/doc/Mesh_3/Concepts/MeshCriteria_3.h b/Mesh_3/doc/Mesh_3/Concepts/MeshCriteria_3.h
@@ -13,7 +13,7 @@ are described by the concept `MeshFacetCriteria_3`.
 The concept `MeshCriteria_3` encapsulates these concepts.
 
 \cgalHasModelsBegin
-\cgalHasModels{CGAL::Mesh_criteria_3<Tr>}
+\cgalHasModels{CGAL::Mesh_criteria_3}
 \cgalHasModelsEnd
 
 \sa `MeshEdgeCriteria_3`

diff --git a/Mesh_3/doc/Mesh_3/PackageDescription.txt b/Mesh_3/doc/Mesh_3/PackageDescription.txt
@@ -92,7 +92,7 @@ related to the template parameters of some models of the main concepts:
 - `CGAL::Mesh_vertex_base_3<GT,MD,Vb>`
 - `CGAL::Compact_mesh_cell_base_3<GT,MD,Tds>`
 - `CGAL::Mesh_cell_base_3<GT,MD,Cb>`
-- `CGAL::Mesh_criteria_3<Tr>`
+- `CGAL::Mesh_criteria_3<C3T3>`
 - `CGAL::Mesh_cell_criteria_3<Tr>`
 - `CGAL::Mesh_facet_criteria_3<Tr>`
 - `CGAL::Mesh_edge_criteria_3<Tr>`

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_gray_image.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_gray_image.cpp
@@ -21,7 +21,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain>::type Tr;
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 namespace params = CGAL::parameters;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_gray_image_multiple_values.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_gray_image_multiple_values.cpp
@@ -24,7 +24,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain>::type Tr;
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 namespace params = CGAL::parameters;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_gray_vtk_image.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_gray_vtk_image.cpp
@@ -34,7 +34,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain>::type Tr;
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 template <bool less_than_iso>
 class Is_interior {

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image.cpp
@@ -24,7 +24,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::t
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 namespace params = CGAL::parameters;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image_variable_size.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image_variable_size.cpp
@@ -25,7 +25,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::t
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 typedef CGAL::Mesh_constant_domain_field_3<Mesh_domain::R,
                                            Mesh_domain::Index> Sizing_field;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image_with_custom_initialization.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image_with_custom_initialization.cpp
@@ -30,7 +30,7 @@ using Tr = CGAL::Mesh_triangulation_3<Mesh_domain, CGAL::Default, Concurrency_ta
 using C3t3 = CGAL::Mesh_complex_3_in_triangulation_3<Tr>;
 
 // Criteria
-using Mesh_criteria = CGAL::Mesh_criteria_3<Tr>;
+using Mesh_criteria = CGAL::Mesh_criteria_3<C3t3>;
 
 namespace params = CGAL::parameters;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image_with_detection_of_features.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image_with_detection_of_features.cpp
@@ -32,7 +32,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::t
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 // To avoid verbose function and named parameters call
 namespace params = CGAL::parameters;

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image_with_features.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image_with_features.cpp
@@ -32,7 +32,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::t
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 namespace params = CGAL::parameters;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image_with_image_initialization.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image_with_image_initialization.cpp
@@ -32,7 +32,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::t
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 namespace params = CGAL::parameters;
 

diff --git a/Mesh_3/examples/Mesh_3/mesh_3D_image_with_initial_points.cpp b/Mesh_3/examples/Mesh_3/mesh_3D_image_with_initial_points.cpp
@@ -30,7 +30,7 @@ typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::t
 typedef CGAL::Mesh_complex_3_in_triangulation_3<Tr> C3t3;
 
 // Criteria
-typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
+typedef CGAL::Mesh_criteria_3<C3t3> Mesh_criteria;
 
 namespace params = CGAL::parameters;