rng.h

/*
 * See Licensing and Copyright notice in naev.h
 */
#pragma once
 
#define RNG(L,H)  (((L)>(H)) ? RNG_BASE((H),(L)) : RNG_BASE((L),(H))) /* L <= RNG <= H */
#define RNG_BASE(L,H) ((int)L + (int)((double)(H-L+1) * randfp())) /* L <= RNG <= H */
#define RNGF()    (randfp()) /* 0. <= RNGF <= 1. */
#define RNG_1SIGMA()       NormalInverse(0.158655255 + RNGF()*(1.-0.158655255*2.))
#define RNG_2SIGMA()       NormalInverse(0.022750132 + RNGF()*(1.-0.022750132*2.))
#define RNG_3SIGMA()       NormalInverse(0.0013498985 + RNGF()*(1.-0.0013498985*2.))
 
/* Init */
void rng_init (void);
 
/* Random functions */
unsigned int randint (void);
double randfp (void);
 
/* Probability functions */
double Normal( double x );
double NormalInverse( double p );