diff options
Diffstat (limited to 'lib/rand.c')
-rw-r--r-- | lib/rand.c | 16 |
1 files changed, 8 insertions, 8 deletions
@@ -2,14 +2,14 @@ #include "fracture.h" long double rand_dist_pow(const gsl_rng *r, double beta) { - long double x = 0; + long double x = 0; - // underflow means that for very small beta x is sometimes identically zero, - // which causes problems - while (x == 0.0) { - long double y = logl(gsl_rng_uniform_pos(r)) / beta; - x = expl(y); - } + // underflow means that for very small beta x is sometimes identically zero, + // which causes problems + while (x == 0.0) { + long double y = logl(gsl_rng_uniform_pos(r)) / beta; + x = expl(y); + } - return x; + return x; } |