src/base/rng/exponential.c


1
2
3
4
5
6
7
8
9
10
11

#include "internal.h"

/* Returns a random float64 between 0 and 1 */
double
rng·exponential(double lambda)
{
    double f;

    f = rng·random();
    return -math·log(1 - f)/lambda;
}