Do not use the boost::variate_generator

Similar to the change https://github.com/PointCloudLibrary/pcl/pull/2956/files

Author: afabri

Generator/include/CGAL/random_convex_hull_in_disc_2.h

@@ -17,7 +17,7 @@
 #ifndef CGAL_RANDOM_CONVEX_HULL_DISC_H
 #define CGAL_RANDOM_CONVEX_HULL_DISC_H 1
 #include <random>
-#include <boost/random.hpp>
+//#include <boost/random.hpp>
 #include <CGAL/Polygon_2_algorithms.h>
 #include <CGAL/function_objects.h>
 #include <CGAL/copy_n.h>
@@ -59,9 +59,9 @@ void generate_points_annulus(long n, double a, double b, double small_radius,
                              GEN& gen) {  // generate n points between a and b
   if (n > 1) {
     std::binomial_distribution<long> bin_distribution(n, .5);
-    boost::variate_generator<GEN&, std::binomial_distribution<long> >
-        g(gen, bin_distribution);
-    long nb = g();
+    //boost::variate_generator<GEN&, std::binomial_distribution<long> >
+    //    g(gen, bin_distribution);
+    long nb = bin_distribution(gen);
     generate_points_annulus(nb, a, (a + b) / 2.0, small_radius, big_radius, l,
                             gen);
     generate_points_annulus(n - nb, (a + b) / 2.0, b, small_radius, big_radius,
@@ -73,13 +73,13 @@ void generate_points_annulus(long n, double a, double b, double small_radius,
         small_radius * small_radius / (big_radius * big_radius), 1);
 
     std::uniform_real_distribution<double> random_angle_distribution(a, b);
-    boost::random::variate_generator<
-        GEN&, std::uniform_real_distribution<double> > random_angle(gen, random_angle_distribution);
-    boost::random::variate_generator<
-        GEN&, std::uniform_real_distribution<double> > random_squared_radius(gen, random_squared_radius_distribution);
+    //boost::random::variate_generator<
+     //   GEN&, std::uniform_real_distribution<double> > random_angle(gen, random_angle_distribution);
+    //boost::random::variate_generator<
+     //   GEN&, std::uniform_real_distribution<double> > random_squared_radius(gen, random_squared_radius_distribution);
 
-    double alpha = random_angle();
-    double r = big_radius * std::sqrt(random_squared_radius());
+    double alpha = random_angle_distribution(gen);
+    double r = big_radius * std::sqrt(random_squared_radius_distribution(gen));
     typedef  Creator_uniform_2<double, P> Creator;
     Creator creator;
     typedef typename Creator::argument_type T;
@@ -219,11 +219,11 @@ void random_convex_hull_in_disc_2(std::size_t n, double radius, std::list<typena
       nb = static_cast<long>((std::min)(T, n - simulated_points));
     }
     std::binomial_distribution<long> dbin(nb, to_double(p_disc));
-    boost::variate_generator<
-        GEN&, std::binomial_distribution<long> > bin(gen, dbin);
+    //boost::variate_generator<
+    //    GEN&, std::binomial_distribution<long> > bin(gen, dbin);
 
     // How many points are falling in the small disc and won't be generated:
-    long k_disc = bin();
+    long k_disc = dbin(gen);
     simulated_points += k_disc;
 
     std::list<P> m;

Generator/test/Generator/random_hull_test.cpp

@@ -16,7 +16,7 @@ main( )
 
   Polygon_2 p,q;
   int n( 1000);
-  boost::mt19937 gen(0);
+  std::mt19937 gen(0);
 
   // build random hull from n random points in a disc:
   random_convex_hull_in_disc_2(n,1.0,gen,std::back_inserter(p),K(), true);