| /** @file |
| |
| A C-program for MT19937-64 (2004/9/29 version). |
| Coded by Takuji Nishimura and Makoto Matsumoto. |
| |
| This is a 64-bit version of Mersenne Twister pseudorandom number |
| generator. |
| |
| Before using, initialize the state by using init_genrand64(seed) |
| or init_by_array64(init_key, key_length). |
| |
| Copyright (C) 2004, Makoto Matsumoto and Takuji Nishimura, |
| All rights reserved. |
| |
| Redistribution and use in source and binary forms, with or without |
| modification, are permitted provided that the following conditions |
| are met: |
| |
| 1. Redistributions of source code must retain the above copyright |
| notice, this list of conditions and the following disclaimer. |
| |
| 2. Redistributions in binary form must reproduce the above copyright |
| notice, this list of conditions and the following disclaimer in the |
| documentation and/or other materials provided with the distribution. |
| |
| 3. The names of its contributors may not be used to endorse or promote |
| products derived from this software without specific prior written |
| permission. |
| |
| THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
| "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
| LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR |
| A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR |
| CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
| EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
| PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
| PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF |
| LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING |
| NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS |
| SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| |
| References: |
| T. Nishimura, ``Tables of 64-bit Mersenne Twisters'' |
| ACM Transactions on Modeling and |
| Computer Simulation 10. (2000) 348--357. |
| M. Matsumoto and T. Nishimura, |
| ``Mersenne Twister: a 623-dimensionally equidistributed |
| uniform pseudorandom number generator'' |
| ACM Transactions on Modeling and |
| Computer Simulation 8. (Jan. 1998) 3--30. |
| |
| Any feedback is very welcome. |
| http://www.math.hiroshima-u.ac.jp/~m-mat/MT/emt.html |
| email: m-mat @ math.sci.hiroshima-u.ac.jp (remove spaces) |
| */ |
| |
| #include "tscore/ink_platform.h" |
| #include "tscore/ink_rand.h" |
| |
| #define NN 312 |
| #define MM 156 |
| #define MATRIX_A 0xB5026F5AA96619E9ULL |
| #define UM 0xFFFFFFFF80000000ULL /* Most significant 33 bits */ |
| #define LM 0x7FFFFFFFULL /* Least significant 31 bits */ |
| |
| static uint64_t mag01[2] = {0ULL, MATRIX_A}; |
| |
| InkRand::InkRand(uint64_t d) |
| { |
| seed(d); |
| } |
| |
| void |
| InkRand::seed(uint64_t seed) |
| { |
| mt[0] = seed; |
| for (mti = 1; mti < NN; mti++) { |
| mt[mti] = (6364136223846793005ULL * (mt[mti - 1] ^ (mt[mti - 1] >> 62)) + mti); |
| } |
| } |
| |
| uint64_t |
| InkRand::random() |
| { |
| int i; |
| uint64_t x; |
| |
| if (mti >= NN) { /* generate NN words at one time */ |
| for (i = 0; i < NN - MM; i++) { |
| x = (mt[i] & UM) | (mt[i + 1] & LM); |
| mt[i] = mt[i + MM] ^ (x >> 1) ^ mag01[static_cast<int>(x & 1ULL)]; |
| } |
| for (; i < NN - 1; i++) { |
| x = (mt[i] & UM) | (mt[i + 1] & LM); |
| mt[i] = mt[i + (MM - NN)] ^ (x >> 1) ^ mag01[static_cast<int>(x & 1ULL)]; |
| } |
| x = (mt[NN - 1] & UM) | (mt[0] & LM); |
| mt[NN - 1] = mt[MM - 1] ^ (x >> 1) ^ mag01[static_cast<int>(x & 1ULL)]; |
| |
| mti = 0; |
| } |
| |
| x = mt[mti++]; |
| |
| x ^= (x >> 29) & 0x5555555555555555ULL; |
| x ^= (x << 17) & 0x71D67FFFEDA60000ULL; |
| x ^= (x << 37) & 0xFFF7EEE000000000ULL; |
| x ^= (x >> 43); |
| |
| return x; |
| } |
| |
| double |
| InkRand::drandom() |
| { |
| return (random() >> 11) * (1.0 / 9007199254740991.0); |
| } |