Distributed QUEST for GPU
Functions
mt19937ar.h File Reference

Go to the source code of this file.

Functions

void init_by_array (unsigned long init_key[], int key_length)
 
void init_genrand (unsigned long s)
 
double genrand_real1 (void)
 
double genrand_real2 (void)
 
double genrand_real3 (void)
 
double genrand_res53 (void)
 

Function Documentation

◆ genrand_real1()

double genrand_real1 ( void  )

Definition at line 146 of file mt19937ar.cpp.

References genrand_int32().

147 {
148  return genrand_int32()*(1.0/4294967295.0);
149  /* divided by 2^32-1 */
150 }
unsigned long genrand_int32(void)
Definition: mt19937ar.cpp:102

◆ genrand_real2()

double genrand_real2 ( void  )

Definition at line 153 of file mt19937ar.cpp.

References genrand_int32().

154 {
155  return genrand_int32()*(1.0/4294967296.0);
156  /* divided by 2^32 */
157 }
unsigned long genrand_int32(void)
Definition: mt19937ar.cpp:102

◆ genrand_real3()

double genrand_real3 ( void  )

Definition at line 160 of file mt19937ar.cpp.

References genrand_int32().

161 {
162  return (((double)genrand_int32()) + 0.5)*(1.0/4294967296.0);
163  /* divided by 2^32 */
164 }
unsigned long genrand_int32(void)
Definition: mt19937ar.cpp:102

◆ genrand_res53()

double genrand_res53 ( void  )

Definition at line 167 of file mt19937ar.cpp.

References genrand_int32().

168 {
169  unsigned long a=genrand_int32()>>5, b=genrand_int32()>>6;
170  return(a*67108864.0+b)*(1.0/9007199254740992.0);
171 }
unsigned long genrand_int32(void)
Definition: mt19937ar.cpp:102

◆ init_by_array()

void init_by_array ( unsigned long  init_key[],
int  key_length 
)

Definition at line 76 of file mt19937ar.cpp.

References init_genrand(), mt, and N.

Referenced by QuESTSeedRandom(), and QuESTSeedRandomDefault().

77 {
78  int i, j, k;
79  init_genrand(19650218UL);
80  i=1; j=0;
81  k = (N>key_length ? N : key_length);
82  for (; k; k--) {
83  mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525UL))
84  + init_key[j] + j; /* non linear */
85  mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */
86  i++; j++;
87  if (i>=N) { mt[0] = mt[N-1]; i=1; }
88  if (j>=key_length) j=0;
89  }
90  for (k=N-1; k; k--) {
91  mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941UL))
92  - i; /* non linear */
93  mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */
94  i++;
95  if (i>=N) { mt[0] = mt[N-1]; i=1; }
96  }
97 
98  mt[0] = 0x80000000UL; /* MSB is 1; assuring non-zero initial array */
99 }
#define N
Definition: mt19937ar.cpp:47
static unsigned long mt[N]
Definition: mt19937ar.cpp:53
void init_genrand(unsigned long s)
Definition: mt19937ar.cpp:57

◆ init_genrand()

void init_genrand ( unsigned long  s)

Definition at line 57 of file mt19937ar.cpp.

References mt, mti, and N.

Referenced by genrand_int32(), and init_by_array().

58 {
59  mt[0]= s & 0xffffffffUL;
60  for (mti=1; mti<N; mti++) {
61  mt[mti] =
62  (1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti);
63  /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */
64  /* In the previous versions, MSBs of the seed affect */
65  /* only MSBs of the array mt[]. */
66  /* 2002/01/09 modified by Makoto Matsumoto */
67  mt[mti] &= 0xffffffffUL;
68  /* for >32 bit machines */
69  }
70 }
#define N
Definition: mt19937ar.cpp:47
static int mti
Definition: mt19937ar.cpp:54
static unsigned long mt[N]
Definition: mt19937ar.cpp:53