| OLD | NEW |
| 1 /* | 1 /* |
| 2 A C-program for MT19937, with initialization improved 2002/1/26. | 2 A C-program for MT19937, with initialization improved 2002/1/26. |
| 3 Coded by Takuji Nishimura and Makoto Matsumoto. | 3 Coded by Takuji Nishimura and Makoto Matsumoto. |
| 4 | 4 |
| 5 Before using, initialize the state by using init_genrand(seed) | 5 Before using, initialize the state by using init_genrand(seed) |
| 6 or init_by_array(init_key, key_length). | 6 or init_by_array(init_key, key_length). |
| 7 | 7 |
| 8 Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, | 8 Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, |
| 9 All rights reserved. | 9 All rights reserved. |
| 10 | 10 |
| (...skipping 39 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 50 #define UPPER_MASK 0x80000000UL /* most significant w-r bits */ | 50 #define UPPER_MASK 0x80000000UL /* most significant w-r bits */ |
| 51 #define LOWER_MASK 0x7fffffffUL /* least significant r bits */ | 51 #define LOWER_MASK 0x7fffffffUL /* least significant r bits */ |
| 52 | 52 |
| 53 MersenneTwister::MersenneTwister() : mt(N), mti(N+1) { | 53 MersenneTwister::MersenneTwister() : mt(N), mti(N+1) { |
| 54 } | 54 } |
| 55 | 55 |
| 56 MersenneTwister::~MersenneTwister() { | 56 MersenneTwister::~MersenneTwister() { |
| 57 } | 57 } |
| 58 | 58 |
| 59 /* initializes mt[N] with a seed */ | 59 /* initializes mt[N] with a seed */ |
| 60 void MersenneTwister::init_genrand(uint32 s) | 60 void MersenneTwister::init_genrand(uint32_t s) |
| 61 { | 61 { |
| 62 mt[0]= s & 0xffffffffUL; | 62 mt[0]= s & 0xffffffffUL; |
| 63 for (mti=1; mti<N; mti++) { | 63 for (mti=1; mti<N; mti++) { |
| 64 mt[mti] = | 64 mt[mti] = |
| 65 (1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti); | 65 (1812433253UL * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti); |
| 66 /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */ | 66 /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */ |
| 67 /* In the previous versions, MSBs of the seed affect */ | 67 /* In the previous versions, MSBs of the seed affect */ |
| 68 /* only MSBs of the array mt[]. */ | 68 /* only MSBs of the array mt[]. */ |
| 69 /* 2002/01/09 modified by Makoto Matsumoto */ | 69 /* 2002/01/09 modified by Makoto Matsumoto */ |
| 70 mt[mti] &= 0xffffffffUL; | 70 mt[mti] &= 0xffffffffUL; |
| 71 /* for >32 bit machines */ | 71 /* for >32 bit machines */ |
| 72 } | 72 } |
| 73 } | 73 } |
| 74 | 74 |
| 75 /* initialize by an array with array-length */ | 75 /* initialize by an array with array-length */ |
| 76 /* init_key is the array for initializing keys */ | 76 /* init_key is the array for initializing keys */ |
| 77 /* key_length is its length */ | 77 /* key_length is its length */ |
| 78 /* slight change for C++, 2004/2/26 */ | 78 /* slight change for C++, 2004/2/26 */ |
| 79 void MersenneTwister::init_by_array(uint32 init_key[], int key_length) | 79 void MersenneTwister::init_by_array(uint32_t init_key[], int key_length) |
| 80 { | 80 { |
| 81 int i, j, k; | 81 int i, j, k; |
| 82 init_genrand(19650218UL); | 82 init_genrand(19650218UL); |
| 83 i=1; j=0; | 83 i=1; j=0; |
| 84 k = (N>key_length ? N : key_length); | 84 k = (N>key_length ? N : key_length); |
| 85 for (; k; k--) { | 85 for (; k; k--) { |
| 86 mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525UL)) | 86 mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525UL)) |
| 87 + init_key[j] + j; /* non linear */ | 87 + init_key[j] + j; /* non linear */ |
| 88 mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */ | 88 mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */ |
| 89 i++; j++; | 89 i++; j++; |
| 90 if (i>=N) { mt[0] = mt[N-1]; i=1; } | 90 if (i>=N) { mt[0] = mt[N-1]; i=1; } |
| 91 if (j>=key_length) j=0; | 91 if (j>=key_length) j=0; |
| 92 } | 92 } |
| 93 for (k=N-1; k; k--) { | 93 for (k=N-1; k; k--) { |
| 94 mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941UL)) | 94 mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941UL)) |
| 95 - i; /* non linear */ | 95 - i; /* non linear */ |
| 96 mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */ | 96 mt[i] &= 0xffffffffUL; /* for WORDSIZE > 32 machines */ |
| 97 i++; | 97 i++; |
| 98 if (i>=N) { mt[0] = mt[N-1]; i=1; } | 98 if (i>=N) { mt[0] = mt[N-1]; i=1; } |
| 99 } | 99 } |
| 100 | 100 |
| 101 mt[0] = 0x80000000UL; /* MSB is 1; assuring non-zero initial array */ | 101 mt[0] = 0x80000000UL; /* MSB is 1; assuring non-zero initial array */ |
| 102 } | 102 } |
| 103 | 103 |
| 104 /* generates a random number on [0,0xffffffff]-interval */ | 104 /* generates a random number on [0,0xffffffff]-interval */ |
| 105 uint32 MersenneTwister::genrand_int32(void) | 105 uint32_t MersenneTwister::genrand_int32(void) |
| 106 { | 106 { |
| 107 uint32 y; | 107 uint32_t y; |
| 108 static uint32 mag01[2]={0x0UL, MATRIX_A}; | 108 static uint32_t mag01[2]={0x0UL, MATRIX_A}; |
| 109 /* mag01[x] = x * MATRIX_A for x=0,1 */ | 109 /* mag01[x] = x * MATRIX_A for x=0,1 */ |
| 110 | 110 |
| 111 if (mti >= N) { /* generate N words at one time */ | 111 if (mti >= N) { /* generate N words at one time */ |
| 112 int kk; | 112 int kk; |
| 113 | 113 |
| 114 if (mti == N+1) /* if init_genrand() has not been called, */ | 114 if (mti == N+1) /* if init_genrand() has not been called, */ |
| 115 init_genrand(5489UL); /* a default initial seed is used */ | 115 init_genrand(5489UL); /* a default initial seed is used */ |
| 116 | 116 |
| 117 for (kk=0;kk<N-M;kk++) { | 117 for (kk=0;kk<N-M;kk++) { |
| 118 y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK); | 118 y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK); |
| (...skipping 12 matching lines...) Expand all Loading... |
| 131 y = mt[mti++]; | 131 y = mt[mti++]; |
| 132 | 132 |
| 133 /* Tempering */ | 133 /* Tempering */ |
| 134 y ^= (y >> 11); | 134 y ^= (y >> 11); |
| 135 y ^= (y << 7) & 0x9d2c5680UL; | 135 y ^= (y << 7) & 0x9d2c5680UL; |
| 136 y ^= (y << 15) & 0xefc60000UL; | 136 y ^= (y << 15) & 0xefc60000UL; |
| 137 y ^= (y >> 18); | 137 y ^= (y >> 18); |
| 138 | 138 |
| 139 return y; | 139 return y; |
| 140 } | 140 } |
| OLD | NEW |