24#ifndef BZ_RANDOM_NORMAL
25#define BZ_RANDOM_NORMAL
27#ifndef BZ_RANDOM_UNIFORM
33template<
typename T = double,
typename IRNG =
defaultIRNG,
47 const T s = 0.449871, t = -0.386595, a = 0.19600, b = 0.25472;
48 const T r1 = 0.27597, r2 = 0.27846;
60 v = 1.715527769921413592960379282557544956242L
66 T q = x*x + y*(a*y - b*x);
77 if (v*v <= -4.0 * log(u)*u*u)
87template<
typename T = double,
typename IRNG =
defaultIRNG,
100 Normal(T mean, T standardDeviation,
unsigned int i) :
NormalUnit()
Definition: normal.h:40
NormalUnit(unsigned int i)
Definition: normal.h:42
T random()
Definition: normal.h:45
T T_numtype
Definition: normal.h:38
T T_numtype
Definition: normal.h:92
T mean_
Definition: normal.h:114
Normal(T mean, T standardDeviation)
Definition: normal.h:94
T standardDeviation_
Definition: normal.h:115
Normal(T mean, T standardDeviation, unsigned int i)
Definition: normal.h:100
T random()
Definition: normal.h:107
sharedState defaultState
Definition: default.h:55
MersenneTwister defaultIRNG
Definition: default.h:120