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 );