24 #ifndef BZ_RANDOM_NORMAL
25 #define BZ_RANDOM_NORMAL
27 #ifndef BZ_RANDOM_UNIFORM
33 template<
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)
87 template<
typename T = double,
typename IRNG =
defaultIRNG,
100 Normal(T mean, T standardDeviation,
unsigned int i) :
120 #endif // BZ_RANDOM_NORMAL
_bz_global blitz::IndexPlaceholder< 11 > t
Definition: indexexpr.h:267
T standardDeviation_
Definition: normal.h:115
_bz_global blitz::IndexPlaceholder< 0 > i
Definition: indexexpr.h:256
NormalUnit()
Definition: normal.h:40
_bz_global blitz::IndexPlaceholder< 8 > q
Definition: indexexpr.h:264
Normal(T mean, T standardDeviation)
Definition: normal.h:94
NormalUnit(unsigned int i)
Definition: normal.h:42
Normal(T mean, T standardDeviation, unsigned int i)
Definition: normal.h:100
T T_numtype
Definition: normal.h:92
MersenneTwister defaultIRNG
Definition: default.h:120
N_length & a
Definition: tvecglobs.h:47
const T2 & b
Definition: minmax.h:48
T mean_
Definition: normal.h:114
_bz_global blitz::IndexPlaceholder< 10 > s
Definition: indexexpr.h:266
T T_numtype
Definition: normal.h:38
T random()
Definition: normal.h:45
sharedState defaultState
Definition: default.h:55
T random()
Definition: normal.h:107