diff --git a/Isosurfacing_3/examples/Isosurfacing_3/contouring_discrete_data.cpp b/Isosurfacing_3/examples/Isosurfacing_3/contouring_discrete_data.cpp
index 0ea0a4817ff..7b35dd37d2d 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/contouring_discrete_data.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/contouring_discrete_data.cpp
@@ -54,7 +54,7 @@ void run_marching_cubes(const Grid& grid,
   Polygon_range triangles;
 
   // run marching cubes isosurfacing
-  IS::marching_cubes(domain, isovalue, points, triangles);
+  IS::marching_cubes<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles);
 
   std::cout << "Output #vertices: " << points.size() << std::endl;
   std::cout << "Output #triangles: " << triangles.size() << std::endl;
@@ -95,7 +95,7 @@ void run_dual_contouring(const Grid& grid,
   Polygon_range triangles;
 
   // run dual contouring isosurfacing
-  IS::dual_contouring(domain, isovalue, points, triangles);
+  IS::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles);
 
   std::cout << "Output #vertices: " << points.size() << std::endl;
   std::cout << "Output #triangles: " << triangles.size() << std::endl;
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/contouring_implicit_data.cpp b/Isosurfacing_3/examples/Isosurfacing_3/contouring_implicit_data.cpp
index b0c44a42e16..fb9728dd000 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/contouring_implicit_data.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/contouring_implicit_data.cpp
@@ -70,7 +70,7 @@ void run_marching_cubes(const Grid& grid,
   Polygon_range triangles;
 
   // run Marching Cubes
-  CGAL::Isosurfacing::marching_cubes(domain, isovalue, points, triangles);
+  CGAL::Isosurfacing::marching_cubes<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles);
 
   timer.stop();
 
@@ -101,8 +101,8 @@ void run_dual_contouring(const Grid& grid,
   Polygon_range triangles;
 
   // run Dual Contouring
-  CGAL::Isosurfacing::dual_contouring<CGAL::Sequential_tag>(domain, isovalue, points, triangles,
-                                                            CGAL::parameters::do_not_triangulate_faces(true));
+  CGAL::Isosurfacing::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                                       CGAL::parameters::do_not_triangulate_faces(true));
 
   timer.stop();
 
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/contouring_inrimage.cpp b/Isosurfacing_3/examples/Isosurfacing_3/contouring_inrimage.cpp
index b596b3d6b0c..ace05f79030 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/contouring_inrimage.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/contouring_inrimage.cpp
@@ -48,7 +48,7 @@ void run_marching_cubes(const Grid& grid,
   Polygon_range triangles;
 
   // execute marching cubes
-  IS::marching_cubes(domain, isovalue, points, triangles);
+  IS::marching_cubes<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles);
 
   std::cout << "Output #vertices: " << points.size() << std::endl;
   std::cout << "Output #triangles: " << triangles.size() << std::endl;
@@ -75,7 +75,7 @@ void run_dual_contouring(const Grid& grid,
   Polygon_range triangles;
 
   // run dual contouring isosurfacing
-  IS::dual_contouring(domain, isovalue, points, triangles);
+  IS::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles);
 
   std::cout << "Output #vertices: " << points.size() << std::endl;
   std::cout << "Output #triangles: " << triangles.size() << std::endl;
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/contouring_mesh_offset.cpp b/Isosurfacing_3/examples/Isosurfacing_3/contouring_mesh_offset.cpp
index 10ef2c75807..5408e38b9c5 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/contouring_mesh_offset.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/contouring_mesh_offset.cpp
@@ -84,7 +84,7 @@ void run_marching_cubes(const Grid& grid,
 
   // run marching cubes
   std::cout << "Running Marching Cubes with isovalue = " << offset_value << std::endl;
-  CGAL::Isosurfacing::marching_cubes(domain, offset_value, points, triangles);
+  CGAL::Isosurfacing::marching_cubes<CGAL::Parallel_if_available_tag>(domain, offset_value, points, triangles);
 
   std::cout << "Output #vertices (MC): " << points.size() << std::endl;
   std::cout << "Output #triangles (MC): " << triangles.size() << std::endl;
@@ -115,8 +115,9 @@ void run_dual_contouring(const Grid& grid,
 
   // run dual contouring
   std::cout << "Running Dual Contouring with isovalue = " << offset_value << std::endl;
-  CGAL::Isosurfacing::dual_contouring(domain, offset_value, points, triangles,
-                                      CGAL::parameters::do_not_triangulate_faces(true));
+  CGAL::Isosurfacing::dual_contouring<CGAL::Parallel_if_available_tag>(
+                        domain, offset_value, points, triangles,
+                        CGAL::parameters::do_not_triangulate_faces(true));
 
   std::cout << "Output #vertices (DC): " << points.size() << std::endl;
   std::cout << "Output #triangles (DC): " << triangles.size() << std::endl;
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring.cpp b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring.cpp
index 171d5a44a8c..07c250117e0 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring.cpp
@@ -65,9 +65,9 @@ int main(int argc, char** argv)
 
   // run dual contouring isosurfacing
   std::cout << "Running Dual Contouring with isovalue = " << isovalue << std::endl;
-  CGAL::Isosurfacing::dual_contouring(domain, isovalue, points, triangles,
-                                      CGAL::parameters::do_not_triangulate_faces(true)
-                                                       .constrain_to_cell(false));
+  CGAL::Isosurfacing::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                                       CGAL::parameters::do_not_triangulate_faces(true)
+                                                                                        .constrain_to_cell(false));
 
   std::cout << "Output #vertices: " << points.size() << std::endl;
   std::cout << "Output #triangles: " << triangles.size() << std::endl;
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_intersection_oracles.cpp b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_intersection_oracles.cpp
index e8b9c34670b..2861d651a55 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_intersection_oracles.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_intersection_oracles.cpp
@@ -73,8 +73,8 @@ int main(int argc, char** argv)
     Polygon_range triangles;
 
     std::cout << "--- Dual Contouring (Dichotomy)" << std::endl;
-    IS::dual_contouring(domain, isovalue, points, triangles,
-                        CGAL::parameters::do_not_triangulate_faces(!triangulate_faces));
+    IS::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                         CGAL::parameters::do_not_triangulate_faces(!triangulate_faces));
 
     timer.stop();
 
@@ -96,8 +96,8 @@ int main(int argc, char** argv)
     Polygon_range triangles;
 
     std::cout << "--- Dual Contouring (Linear Interpolation)" << std::endl;
-    IS::dual_contouring(domain, isovalue, points, triangles,
-                        CGAL::parameters::do_not_triangulate_faces(!triangulate_faces));
+    IS::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                         CGAL::parameters::do_not_triangulate_faces(!triangulate_faces));
 
     timer.stop();
 
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_octree.cpp b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_octree.cpp
index facc4b11709..6198f5b2045 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_octree.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_octree.cpp
@@ -117,8 +117,8 @@ int main(int argc, char** argv)
   Polygon_range triangles;
 
   // run Dual Contouring
-  CGAL::Isosurfacing::dual_contouring<CGAL::Sequential_tag>(domain, isovalue, points, triangles,
-                                                            CGAL::parameters::do_not_triangulate_faces(true));
+  CGAL::Isosurfacing::dual_contouring<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                                       CGAL::parameters::do_not_triangulate_faces(true));
 
   timer.stop();
 
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_strategies.cpp b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_strategies.cpp
index b76922e8ce9..58e91a8271c 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_strategies.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/dual_contouring_strategies.cpp
@@ -69,7 +69,7 @@ int main(int argc, char** argv)
     Polygon_range triangles;
 
     std::cout << "--- Dual Contouring (QEM - unconstrained)" << std::endl;
-    IS::internal::Dual_contourer<CGAL::Sequential_tag, Domain,
+    IS::internal::Dual_contourer<CGAL::Parallel_if_available_tag, Domain,
                                  IS::internal::DC_Strategy::QEM> contourer;
     contourer(domain, isovalue, points, triangles,
               CGAL::parameters::do_not_triangulate_faces(!triangulate_faces));
@@ -91,7 +91,7 @@ int main(int argc, char** argv)
     Polygon_range triangles;
 
     std::cout << "--- Dual Contouring (QEM - constrained)" << std::endl;
-    IS::internal::Dual_contourer<CGAL::Sequential_tag, Domain,
+    IS::internal::Dual_contourer<CGAL::Parallel_if_available_tag, Domain,
                                  IS::internal::DC_Strategy::QEM> contourer;
     contourer(domain, isovalue, points, triangles,
               CGAL::parameters::do_not_triangulate_faces(!triangulate_faces)
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes.cpp b/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes.cpp
index 8cc6ac9cf48..2ca5984a58a 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes.cpp
@@ -50,8 +50,8 @@ int main(int argc, char** argv)
 
   // run marching cubes isosurfacing
   std::cout << "Running Marching Cubes with isovalue = " << isovalue << std::endl;
-  CGAL::Isosurfacing::marching_cubes(domain, isovalue, points, triangles,
-                                     CGAL::parameters::use_topologically_correct_marching_cubes(true));
+  CGAL::Isosurfacing::marching_cubes<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                                      CGAL::parameters::use_topologically_correct_marching_cubes(true));
 
   std::cout << "Output #vertices: " << points.size() << std::endl;
   std::cout << "Output #triangles: " << triangles.size() << std::endl;
diff --git a/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes_strategies.cpp b/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes_strategies.cpp
index 857d0fdf5cd..aad37317e9f 100644
--- a/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes_strategies.cpp
+++ b/Isosurfacing_3/examples/Isosurfacing_3/marching_cubes_strategies.cpp
@@ -53,8 +53,8 @@ int main(int argc, char** argv)
 
     // run marching cubes isosurfacing
     std::cout << "Running Marching Cubes with isovalue = " << isovalue << std::endl;
-    CGAL::Isosurfacing::marching_cubes(domain, isovalue, points, triangles,
-                                      CGAL::parameters::use_topologically_correct_marching_cubes(false));
+    CGAL::Isosurfacing::marching_cubes<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                                        CGAL::parameters::use_topologically_correct_marching_cubes(false));
 
     timer.stop();
 
@@ -74,8 +74,8 @@ int main(int argc, char** argv)
 
     // run marching cubes isosurfacing
     std::cout << "Running Marching Cubes (TMC) with isovalue = " << isovalue << std::endl;
-    CGAL::Isosurfacing::marching_cubes(domain, isovalue, points, triangles,
-                                      CGAL::parameters::use_topologically_correct_marching_cubes(true));
+    CGAL::Isosurfacing::marching_cubes<CGAL::Parallel_if_available_tag>(domain, isovalue, points, triangles,
+                                                                        CGAL::parameters::use_topologically_correct_marching_cubes(true));
 
     timer.stop();