eonc::rng Namespace Reference

Functions
double	random (long newSeed=0)
double	randomDouble ()
double	randomDouble (int max)
double	randomDouble (long max)
double	randomDouble (double max)
long	randomInt (int lower, int upper)
double	gaussRandom (double avg, double std)

Function Documentation

gaussRandom()

double eonc::rng::gaussRandom	(	double	avg,
		double	std )

Definition at line 80 of file RandomNumbers.cpp.

                                                  {
  double r = 2, v1, v2, l, result;
  while (r >= 1.0 || r < 1e-300) {
    v1 = 2.0 * randomDouble() - 1.0;
    v2 = 2.0 * randomDouble() - 1.0;
    r = v1 * v1 + v2 * v2;
  }
  l = v1 * sqrt(-2.0 * ::log(r) / r);
  result = avg + std * l;
  return (result);
}

random()

double eonc::rng::random

(

long

newSeed = 0

)

Definition at line 16 of file RandomNumbers.cpp.

                                   {
  static long seed = -1;
  if (newSeed) {
    seed = -newSeed;
  }
  int j;
  long k;
  static long seed2 = 123456789;
  static long iy = 0;
  static long iv[NTAB];
  double temp;
  if (seed <= 0) {
    if (-(seed) < 1)
      seed = 3;
    else
      seed = -(seed);
    seed2 = (seed);
    for (j = NTAB + 7; j >= 0; j--) {
      k = (seed) / IQ1;
      seed = IA1 * (seed - k * IQ1) - k * IR1;
      if (seed < 0)
        seed += IM1;
      if (j < NTAB)
        iv[j] = seed;
    }
    iy = iv[0];
  }
  k = (seed) / IQ1;
  seed = IA1 * (seed - k * IQ1) - k * IR1;
  if (seed < 0)
    seed += IM1;
  k = seed2 / IQ2;
  seed2 = IA2 * (seed2 - k * IQ2) - k * IR2;
  if (seed2 < 0)
    seed2 += IM2;
  j = int(iy / NDIV);
  iy = iv[j] - seed2;
  iv[j] = seed;
  if (iy < 1)
    iy += IMM1;
  if ((temp = double(AM * iy)) > RNMX)
    return RNMX;
  else
    return temp;
}

randomDouble() [1/4]

double eonc::rng::randomDouble

(

)

Definition at line 62 of file RandomNumbers.cpp.

{ return (random()); }

randomDouble() [2/4]

double eonc::rng::randomDouble

(

double

max

)

Definition at line 74 of file RandomNumbers.cpp.

{ return (dmax * randomDouble()); }

randomDouble() [3/4]

double eonc::rng::randomDouble

(

int

max

)

Definition at line 64 of file RandomNumbers.cpp.

                                    {
  double dmax = double(max);
  return (dmax * randomDouble());
}

randomDouble() [4/4]

double eonc::rng::randomDouble

(

long

max

)

Definition at line 69 of file RandomNumbers.cpp.

                                     {
  double dmax = double(max);
  return (dmax * randomDouble());
}

randomInt()

long eonc::rng::randomInt	(	int	lower,
		int	upper )

Definition at line 76 of file RandomNumbers.cpp.

                                            {
  return lround((upper - lower) * randomDouble() + lower);
}