more seed

Author: nschloe

pygalmesh/main.py

@@ -442,6 +442,7 @@ def generate_from_array_with_subdomain_sizing(
     facet_distance=0.0,
     cell_radius_edge_ratio=0.0,
     verbose=True,
+    seed=0,
 ):
     assert vol.dtype in ["uint8", "uint16"]
     fh, inr_filename = tempfile.mkstemp(suffix=".inr")
@@ -460,6 +461,7 @@ def generate_from_array_with_subdomain_sizing(
         facet_distance,
         cell_radius_edge_ratio,
         verbose,
+        seed,
     )
     os.remove(inr_filename)
     return mesh

src/generate_from_inr.cpp

@@ -47,9 +47,12 @@ generate_from_inr(
     const double facet_distance,
     const double cell_radius_edge_ratio,
     const double cell_size,
-    const bool verbose
+    const bool verbose,
+    const int seed
     )
 {
+  CGAL::get_default_random() = CGAL::Random(seed);
+
   CGAL::Image_3 image;
   const bool success = image.read(inr_filename.c_str());
   if (!success) {

src/generate_from_inr.hpp

@@ -39,7 +39,8 @@ generate_from_inr_with_subdomain_sizing(
     const double facet_size = 0.0,
     const double facet_distance = 0.0,
     const double cell_radius_edge_ratio = 0.0,
-    const bool verbose  = true
+    const bool verbose  = true,
+    const int seed = 0
     );
 
 } // namespace pygalmesh