[gnuastro-commits] master 3823701 2/3: Default random number generator set to gsl_rng_ranlxs1

[Mohammad Akhlaghi]

branch: master
commit 382370184ab48a0e6c5977050ec7b2b735767e91
Author: Mohammad Akhlaghi <address@hidden>
Commit: Mohammad Akhlaghi <address@hidden>

    Default random number generator set to gsl_rng_ranlxs1
    
    Until now, the three MakeCatalog, MakeNoise and MakeProfile would use
    `gsl_rng_default' as the random number generation algorithm. When
    unchanged, this was the `mt19937' algorithm. However, the GSL documentation
    suggests that "For the most reliable source of uncorrelated numbers, the
    second-generation RANLUX generators have the strongest proof of
    randomness". There are three luxary levels, so the default value is set to
    the middle one (`ranlxs1').
---
 bin/mkcatalog/ui.c | 2 +-
 bin/mknoise/ui.c   | 2 +-
 bin/mkprof/ui.c    | 2 +-
 3 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/bin/mkcatalog/ui.c b/bin/mkcatalog/ui.c
index 0cfb7bf..8c57145 100644
--- a/bin/mkcatalog/ui.c
+++ b/bin/mkcatalog/ui.c
@@ -858,7 +858,7 @@ ui_preparations_upperlimit(struct mkcatalogparams *p)
 
   /* Set the random number generator. */
   gsl_rng_env_setup();
-  p->rng=gsl_rng_alloc(gsl_rng_default);
+  p->rng=gsl_rng_alloc(gsl_rng_ranlxs1);
   p->seed = ( p->envseed
               ? gsl_rng_default_seed
               : gal_timing_time_based_rng_seed() );
diff --git a/bin/mknoise/ui.c b/bin/mknoise/ui.c
index 606e957..98adb7e 100644
--- a/bin/mknoise/ui.c
+++ b/bin/mknoise/ui.c
@@ -310,7 +310,7 @@ ui_preparations(struct mknoiseparams *p)
 
   /* Allocate the random number generator: */
   gsl_rng_env_setup();
-  p->rng=gsl_rng_alloc(gsl_rng_default);
+  p->rng=gsl_rng_alloc(gsl_rng_ranlxs1);
   p->rng_seed = ( p->envseed
                   ? gsl_rng_default_seed
                   : gal_timing_time_based_rng_seed() );
diff --git a/bin/mkprof/ui.c b/bin/mkprof/ui.c
index 8f239ce..8dbc372 100644
--- a/bin/mkprof/ui.c
+++ b/bin/mkprof/ui.c
@@ -1414,7 +1414,7 @@ ui_preparations(struct mkprofparams *p)
 
   /* Allocate the random number generator: */
   gsl_rng_env_setup();
-  p->rng=gsl_rng_alloc(gsl_rng_default);
+  p->rng=gsl_rng_alloc(gsl_rng_ranlxs1);
 
   /* Make the log linked list. */
   ui_make_log(p);