[270] | 1 | /*
|
---|
| 2 | ** mt19937ar.h - MT Random functions
|
---|
| 3 | **
|
---|
[331] | 4 | ** Copyright (C) 1997 - 2016, Makoto Matsumoto and Takuji Nishimura,
|
---|
| 5 | ** All rights reserved.
|
---|
| 6 | **
|
---|
| 7 | ** Permission is hereby granted, free of charge, to any person obtaining
|
---|
| 8 | ** a copy of this software and associated documentation files (the
|
---|
| 9 | ** "Software"), to deal in the Software without restriction, including
|
---|
| 10 | ** without limitation the rights to use, copy, modify, merge, publish,
|
---|
| 11 | ** distribute, sublicense, and/or sell copies of the Software, and to
|
---|
| 12 | ** permit persons to whom the Software is furnished to do so, subject to
|
---|
| 13 | ** the following conditions:
|
---|
| 14 | **
|
---|
| 15 | ** The above copyright notice and this permission notice shall be
|
---|
| 16 | ** included in all copies or substantial portions of the Software.
|
---|
| 17 | **
|
---|
| 18 | ** THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
|
---|
| 19 | ** EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
|
---|
| 20 | ** MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
|
---|
| 21 | ** IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY
|
---|
| 22 | ** CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
|
---|
| 23 | ** TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
|
---|
| 24 | ** SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
|
---|
| 25 | **
|
---|
| 26 | ** [ MIT license: http://www.opensource.org/licenses/mit-license.php ]
|
---|
| 27 | **
|
---|
| 28 | ** Any feedback is very welcome.
|
---|
| 29 | ** http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html
|
---|
| 30 | ** email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space)
|
---|
| 31 | **
|
---|
| 32 | ** This version is modified by mruby developers. If you see any problem,
|
---|
| 33 | ** contact us first at https://github.com/mruby/mruby/issues
|
---|
[270] | 34 | */
|
---|
| 35 |
|
---|
| 36 | #define N 624
|
---|
| 37 |
|
---|
| 38 | typedef struct {
|
---|
| 39 | unsigned long mt[N];
|
---|
| 40 | int mti;
|
---|
| 41 | union {
|
---|
| 42 | unsigned long int_;
|
---|
| 43 | double double_;
|
---|
| 44 | } gen;
|
---|
| 45 |
|
---|
| 46 | mrb_int seed;
|
---|
| 47 | mrb_bool has_seed : 1;
|
---|
| 48 | } mt_state;
|
---|
| 49 |
|
---|
| 50 | void mrb_random_init_genrand(mt_state *, unsigned long);
|
---|
| 51 | unsigned long mrb_random_genrand_int32(mt_state *);
|
---|
| 52 | double mrb_random_genrand_real1(mt_state *t);
|
---|
| 53 |
|
---|
| 54 | /* initializes mt[N] with a seed */
|
---|
| 55 | void init_genrand(unsigned long s);
|
---|
| 56 |
|
---|
| 57 | /* initialize by an array with array-length */
|
---|
| 58 | /* init_key is the array for initializing keys */
|
---|
| 59 | /* key_length is its length */
|
---|
| 60 | /* slight change for C++, 2004/2/26 */
|
---|
| 61 | void init_by_array(unsigned long init_key[], int key_length);
|
---|
| 62 |
|
---|
| 63 | /* generates a random number on [0,0xffffffff]-interval */
|
---|
| 64 | unsigned long genrand_int32(void);
|
---|
| 65 |
|
---|
| 66 | /* generates a random number on [0,0x7fffffff]-interval */
|
---|
| 67 | long genrand_int31(void);
|
---|
| 68 |
|
---|
| 69 | /* These real versions are due to Isaku Wada, 2002/01/09 added */
|
---|
| 70 | /* generates a random number on [0,1]-real-interval */
|
---|
| 71 | double genrand_real1(void);
|
---|
| 72 |
|
---|
| 73 | /* generates a random number on [0,1)-real-interval */
|
---|
| 74 | double genrand_real2(void);
|
---|
| 75 |
|
---|
| 76 | /* generates a random number on (0,1)-real-interval */
|
---|
| 77 | double genrand_real3(void);
|
---|
| 78 |
|
---|
| 79 | /* generates a random number on [0,1) with 53-bit resolution*/
|
---|
| 80 | double genrand_res53(void);
|
---|