uni_rand.cpp

Digital filter designer s handbook C++ code source

//
//  File = uni_rand.cpp
//
#define IA 16807
#define IM 2147483647
#define AM (1.0/IM)
#define IQ 127773
#define IR 2836
#define MASK 123459876

float UniformRandom(long *seed)
{
 long k;
 float result;
 k=(*seed)/IQ;
 *seed = IA * (*seed-k*IQ)-IR*k;
 if(*seed<0) *seed += IM;
 result = float(AM*(*seed));
 
 return(result);
 }
 
 
double DoubleUniformRandom(long *seed)
{
 long k;
 double result;
 k=(*seed)/IQ;
 *seed = IA * (*seed-k*IQ)-IR*k;
 if(*seed<0) *seed += IM;
 result = AM*(*seed);
 
 return(result);
 }