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 90 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 29 of file RandomNumbers.cpp.

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

randomDouble() [1/4]

double eonc::rng::randomDouble

(

)

Definition at line 72 of file RandomNumbers.cpp.

{ return (random()); }

randomDouble() [2/4]

double eonc::rng::randomDouble

(

double

max

)

Definition at line 84 of file RandomNumbers.cpp.

{ return (dmax * randomDouble()); }

randomDouble() [3/4]

double eonc::rng::randomDouble

(

int

max

)

Definition at line 74 of file RandomNumbers.cpp.

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

randomDouble() [4/4]

double eonc::rng::randomDouble

(

long

max

)

Definition at line 79 of file RandomNumbers.cpp.

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

randomInt()

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

Definition at line 86 of file RandomNumbers.cpp.

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