diff options
| -rw-r--r-- | libcrystfel/src/geometry.c | 35 |
1 files changed, 17 insertions, 18 deletions
diff --git a/libcrystfel/src/geometry.c b/libcrystfel/src/geometry.c index e847711d..ff6d14c7 100644 --- a/libcrystfel/src/geometry.c +++ b/libcrystfel/src/geometry.c @@ -194,47 +194,46 @@ double gaussian_fraction(double rlow, double rhigh, double R) } +static void reseed(gsl_rng *rng) +{ + unsigned long int seed = gsl_rng_get(rng); + gsl_rng_set(rng, seed); +} + + static double random_partiality(signed int h, signed int k, signed int l, int serial) { gsl_rng *rng = gsl_rng_alloc(gsl_rng_mt19937); - unsigned long int seed; double p; int i; gsl_rng_set(rng, serial); - seed = gsl_rng_get(rng); - gsl_rng_set(rng, seed); + reseed(rng); for ( i=0; i<abs(h)+1; i++ ) { - seed = gsl_rng_get(rng); + gsl_rng_get(rng); } - gsl_rng_set(rng, seed); + reseed(rng); if ( h >= 0 ) { - seed = gsl_rng_get(rng); + reseed(rng); } - seed = gsl_rng_get(rng); - gsl_rng_set(rng, seed); for ( i=0; i<abs(k)+1; i++ ) { - seed = gsl_rng_get(rng); + gsl_rng_get(rng); } - gsl_rng_set(rng, seed); + reseed(rng); if ( k >= 0 ) { - seed = gsl_rng_get(rng); + reseed(rng); } - seed = gsl_rng_get(rng); - gsl_rng_set(rng, seed); for ( i=0; i<abs(l)+1; i++ ) { - seed = gsl_rng_get(rng); + gsl_rng_get(rng); } - gsl_rng_set(rng, seed); + reseed(rng); if ( l >= 0 ) { - seed = gsl_rng_get(rng); + reseed(rng); } - seed = gsl_rng_get(rng); - gsl_rng_set(rng, seed); p = gsl_rng_uniform(rng); gsl_rng_free(rng); |