44 #ifndef SCI_GAUSSIAN_H__
45 #define SCI_GAUSSIAN_H__
63 Gaussian(
double mean=0,
double sigma=1,
int seed=0);
68 inline double rand() {
return sqrt(-2*log((*mr_)()))*cos(2*
M_PI*(*mr_)())*sigma_+mean_;}
71 double prob(
double x) {
return exp(-(x-mean_)*(x-mean_)/(2*sigma_*sigma_))/(sigma_*sqrt(2*
M_PI));}
76 #endif //SCI_GAUSSIAN_H__
Definition: Gaussian.h:58
#define SCISHARE
Definition: share.h:39
Musil random number generator.
Definition: MusilRNG.h:49
MusilRNG * mr_
Definition: Gaussian.h:62
double mean_
Definition: Gaussian.h:60
#define M_PI
Definition: fft.c:36
double prob(double x)
Definition: Gaussian.h:71
double rand()
Definition: Gaussian.h:68
double sigma_
Definition: Gaussian.h:61