minor change

Author: francoisbonneau

bindings/python/src/stochastic/sampling/direct/double_sampler.hpp

@@ -60,6 +60,11 @@ namespace geode
             .def_static( "create_distribution",
                 &DoubleSampler::create_distribution, pybind11::arg( "desc" ),
                 "Create a distribution from a description" )
+            .def_static( "create_rad_angle_distribution_from_degree",
+                &DoubleSampler::create_rad_angle_distribution_from_degree,
+                pybind11::arg( "desc" ),
+                "Create a angle distribution in radian from a description "
+                "provided in degree" )
             .def_static( "sample", &DoubleSampler::sample,
                 pybind11::arg( "engine" ), pybind11::arg( "dist" ),
                 "Sample a value from a distribution using a RandomEngine" );

bindings/python/src/stochastic/sampling/distributions.hpp

@@ -87,6 +87,17 @@ namespace geode
                 &TruncatedGaussian::distribution_type_static )
             .def( "distribution_type", &TruncatedGaussian::distribution_type )
             .def( "string", &TruncatedGaussian::string );
+
+        // VonMises
+        pybind11::class_< VonMises >( module, "VonMises" )
+            .def( pybind11::init<>() )
+            .def_readwrite( "mean", &VonMises::mean )
+            .def_readwrite( "concentration", &VonMises::concentration )
+            .def( "is_valid", &VonMises::is_valid )
+            .def_static( "distribution_type_static",
+                &VonMises::distribution_type_static )
+            .def( "distribution_type", &VonMises::distribution_type )
+            .def( "string", &VonMises::string );
     }
 
 } // namespace geode

bindings/python/src/stochastic/sampling/random_engine.hpp

@@ -64,6 +64,9 @@ namespace geode
                 &RandomEngine::sample_truncated_gaussian,
                 pybind11::arg( "law" ),
                 "Sample a value from a truncated Gaussian" )
+            .def( "sample_von_mises", &RandomEngine::sample_von_mises,
+                pybind11::arg( "law" ),
+                "Sample a value from a Von Mises-Fisher" )
 
             // Other distributions
             .def( "sample_log", &RandomEngine::sample_log,

bindings/python/tests/stochastic/test-py-mh-fractures.py

@@ -57,7 +57,7 @@ def test_fracture_simulator():
     setA.azimuth.max_value = 10.0
 
     # positioning
-    setA.p20 = 0.1
+    setA.p20 = 0.06
     setA.minimal_spacing = 1.0
 
     runner = stochastic.FractureSimulationRunner(box)
@@ -97,10 +97,10 @@ def test_two_fracture_sets_simulator():
     setA.length.distribution_type = stochastic.DistributionType("UniformClosed")
     setA.length.min_value = 1.0
     setA.length.max_value = 10.0
-    setA.azimuth.distribution_type = stochastic.DistributionType("UniformClosed")
-    setA.azimuth.min_value = 1.0
-    setA.azimuth.max_value = 10.0
-    setA.p20 = 0.1
+    setA.azimuth.distribution_type = stochastic.DistributionType("VonMises")
+    setA.azimuth.mean = 45
+    setA.azimuth.standard_deviation = 10.0
+    setA.p20 = 0.05
     setA.minimal_spacing = 1.0
 
     # --- Object set B
@@ -112,7 +112,7 @@ def test_two_fracture_sets_simulator():
     setB.azimuth.distribution_type =stochastic.DistributionType("UniformClosed")
     setB.azimuth.min_value = 90.0
     setB.azimuth.max_value = 100.0
-    setB.p20 = 0.1
+    setB.p20 = 0.03
     setB.minimal_spacing = 2.0
 
     runner = stochastic.FractureSimulationRunner(box)

tests/stochastic/sampling/mcmc/test-mh-fractures.cpp

@@ -54,7 +54,7 @@ namespace
         setA.azimuth.max_value = 10.;
 
         // positionning
-        setA.p20 = 0.1;
+        setA.p20 = 0.05;
         setA.minimal_spacing = 1.;
 
         geode::FractureSimulationRunner runner( box );
@@ -108,7 +108,7 @@ namespace
         setA.azimuth.max_value = 10.;
 
         // positionning
-        setA.p20 = 0.1;
+        setA.p20 = 0.05;
         setA.minimal_spacing = 1.;
 
         // --- Object set
@@ -128,7 +128,7 @@ namespace
         setB.azimuth.standard_deviation = 15.;
 
         // positionning
-        setB.p20 = 0.1;
+        setB.p20 = 0.05;
         setB.minimal_spacing = 2.;
 
         geode::FractureSimulationRunner runner( box );