random_distr.h File Reference

Go to the source code of this file.

Functions
template<typename T >
T	rand_uniform ()
template<typename T >
T	rand_normal ()
template<typename T >
int	rand_poisson (T lambda)

Function Documentation

rand_normal()

template<typename T >

T rand_normal

(

)

Definition at line 11 of file random_distr.h.

{
    T U0 = rand_uniform<T>()+(T)1e-9;
    T U1 = rand_uniform<T>();

    // Box-Muller transform
    return sqrt( -2 * log(U0) ) * cos(2*(T)3.141593 * U1);
}

rand_poisson()

template<typename T >

int rand_poisson

(

T

lambda

)

Definition at line 21 of file random_distr.h.

{
    if (lambda > 10) {
        T v = rand_normal<T>();
        v = (T)0.5 + lambda + v*sqrt(lambda);
        return (int)std::max((T)0,v);
    }
    else {
        T L = exp(-lambda);
        int k = 0;
        T p = 1;
        do {
            k++;
            T u = rand_uniform<T>();
            p = p * u;
        } while (p > L && k<200);
    
        // copy back
        return k-1;
    }
}

rand_uniform()

template<typename T >

T rand_uniform

(

)

Definition at line 4 of file random_distr.h.

{
    return rand() / (T)RAND_MAX;
}