Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 1 | /* Random objects */ |
| 2 | |
| 3 | /* ------------------------------------------------------------------ |
| 4 | The code in this module was based on a download from: |
Hiroki Noda | f15fa87 | 2017-02-15 18:04:43 +0900 | [diff] [blame] | 5 | http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/MT2002/emt19937ar.html |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 6 | |
| 7 | It was modified in 2002 by Raymond Hettinger as follows: |
| 8 | |
Benjamin Peterson | d8e5f2d | 2010-08-24 18:08:22 +0000 | [diff] [blame] | 9 | * the principal computational lines untouched. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 10 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 11 | * renamed genrand_res53() to random_random() and wrapped |
| 12 | in python calling/return code. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 13 | |
Victor Stinner | 9f5fe79 | 2020-04-17 19:05:35 +0200 | [diff] [blame] | 14 | * genrand_uint32() and the helper functions, init_genrand() |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 15 | and init_by_array(), were declared static, wrapped in |
| 16 | Python calling/return code. also, their global data |
| 17 | references were replaced with structure references. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 18 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 19 | * unused functions from the original were deleted. |
| 20 | new, original C python code was added to implement the |
| 21 | Random() interface. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 22 | |
| 23 | The following are the verbatim comments from the original code: |
| 24 | |
| 25 | A C-program for MT19937, with initialization improved 2002/1/26. |
| 26 | Coded by Takuji Nishimura and Makoto Matsumoto. |
| 27 | |
| 28 | Before using, initialize the state by using init_genrand(seed) |
| 29 | or init_by_array(init_key, key_length). |
| 30 | |
| 31 | Copyright (C) 1997 - 2002, Makoto Matsumoto and Takuji Nishimura, |
| 32 | All rights reserved. |
| 33 | |
| 34 | Redistribution and use in source and binary forms, with or without |
| 35 | modification, are permitted provided that the following conditions |
| 36 | are met: |
| 37 | |
| 38 | 1. Redistributions of source code must retain the above copyright |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 39 | notice, this list of conditions and the following disclaimer. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 40 | |
| 41 | 2. Redistributions in binary form must reproduce the above copyright |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 42 | notice, this list of conditions and the following disclaimer in the |
| 43 | documentation and/or other materials provided with the distribution. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 44 | |
| 45 | 3. The names of its contributors may not be used to endorse or promote |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 46 | products derived from this software without specific prior written |
| 47 | permission. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 48 | |
| 49 | THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS |
| 50 | "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT |
| 51 | LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR |
| 52 | A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR |
| 53 | CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
| 54 | EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
| 55 | PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
| 56 | PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF |
| 57 | LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING |
| 58 | NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS |
| 59 | SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
| 60 | |
| 61 | |
| 62 | Any feedback is very welcome. |
Hiroki Noda | f15fa87 | 2017-02-15 18:04:43 +0900 | [diff] [blame] | 63 | http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html |
| 64 | email: m-mat @ math.sci.hiroshima-u.ac.jp (remove space) |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 65 | */ |
| 66 | |
| 67 | /* ---------------------------------------------------------------*/ |
| 68 | |
| 69 | #include "Python.h" |
Victor Stinner | 9f2a920 | 2016-09-06 17:03:03 -0700 | [diff] [blame] | 70 | #ifdef HAVE_PROCESS_H |
Victor Stinner | 9f5fe79 | 2020-04-17 19:05:35 +0200 | [diff] [blame] | 71 | # include <process.h> // getpid() |
Victor Stinner | 9f2a920 | 2016-09-06 17:03:03 -0700 | [diff] [blame] | 72 | #endif |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 73 | |
| 74 | /* Period parameters -- These are all magic. Don't change. */ |
| 75 | #define N 624 |
| 76 | #define M 397 |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 77 | #define MATRIX_A 0x9908b0dfU /* constant vector a */ |
| 78 | #define UPPER_MASK 0x80000000U /* most significant w-r bits */ |
| 79 | #define LOWER_MASK 0x7fffffffU /* least significant r bits */ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 80 | |
| 81 | typedef struct { |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 82 | PyObject *Random_Type; |
| 83 | PyObject *Long___abs__; |
| 84 | } _randomstate; |
| 85 | |
Hai Shi | f707d94 | 2020-03-16 21:15:01 +0800 | [diff] [blame] | 86 | static inline _randomstate* |
| 87 | get_random_state(PyObject *module) |
| 88 | { |
| 89 | void *state = PyModule_GetState(module); |
| 90 | assert(state != NULL); |
| 91 | return (_randomstate *)state; |
| 92 | } |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 93 | |
| 94 | static struct PyModuleDef _randommodule; |
| 95 | |
Hai Shi | f707d94 | 2020-03-16 21:15:01 +0800 | [diff] [blame] | 96 | #define _randomstate_global get_random_state(PyState_FindModule(&_randommodule)) |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 97 | |
| 98 | typedef struct { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 99 | PyObject_HEAD |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 100 | int index; |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 101 | uint32_t state[N]; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 102 | } RandomObject; |
| 103 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 104 | |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 105 | #include "clinic/_randommodule.c.h" |
| 106 | |
| 107 | /*[clinic input] |
| 108 | module _random |
| 109 | class _random.Random "RandomObject *" "&Random_Type" |
| 110 | [clinic start generated code]*/ |
| 111 | /*[clinic end generated code: output=da39a3ee5e6b4b0d input=f79898ae7847c321]*/ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 112 | |
| 113 | /* Random methods */ |
| 114 | |
| 115 | |
| 116 | /* generates a random number on [0,0xffffffff]-interval */ |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 117 | static uint32_t |
Victor Stinner | 9f5fe79 | 2020-04-17 19:05:35 +0200 | [diff] [blame] | 118 | genrand_uint32(RandomObject *self) |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 119 | { |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 120 | uint32_t y; |
| 121 | static const uint32_t mag01[2] = {0x0U, MATRIX_A}; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 122 | /* mag01[x] = x * MATRIX_A for x=0,1 */ |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 123 | uint32_t *mt; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 124 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 125 | mt = self->state; |
| 126 | if (self->index >= N) { /* generate N words at one time */ |
| 127 | int kk; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 128 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 129 | for (kk=0;kk<N-M;kk++) { |
| 130 | y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK); |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 131 | mt[kk] = mt[kk+M] ^ (y >> 1) ^ mag01[y & 0x1U]; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 132 | } |
| 133 | for (;kk<N-1;kk++) { |
| 134 | y = (mt[kk]&UPPER_MASK)|(mt[kk+1]&LOWER_MASK); |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 135 | mt[kk] = mt[kk+(M-N)] ^ (y >> 1) ^ mag01[y & 0x1U]; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 136 | } |
| 137 | y = (mt[N-1]&UPPER_MASK)|(mt[0]&LOWER_MASK); |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 138 | mt[N-1] = mt[M-1] ^ (y >> 1) ^ mag01[y & 0x1U]; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 139 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 140 | self->index = 0; |
| 141 | } |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 142 | |
| 143 | y = mt[self->index++]; |
| 144 | y ^= (y >> 11); |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 145 | y ^= (y << 7) & 0x9d2c5680U; |
| 146 | y ^= (y << 15) & 0xefc60000U; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 147 | y ^= (y >> 18); |
| 148 | return y; |
| 149 | } |
| 150 | |
| 151 | /* random_random is the function named genrand_res53 in the original code; |
| 152 | * generates a random number on [0,1) with 53-bit resolution; note that |
| 153 | * 9007199254740992 == 2**53; I assume they're spelling "/2**53" as |
| 154 | * multiply-by-reciprocal in the (likely vain) hope that the compiler will |
| 155 | * optimize the division away at compile-time. 67108864 is 2**26. In |
| 156 | * effect, a contains 27 random bits shifted left 26, and b fills in the |
| 157 | * lower 26 bits of the 53-bit numerator. |
Berker Peksag | 0ac70c0 | 2016-04-29 16:54:10 +0300 | [diff] [blame] | 158 | * The original code credited Isaku Wada for this algorithm, 2002/01/09. |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 159 | */ |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 160 | |
| 161 | /*[clinic input] |
| 162 | _random.Random.random |
| 163 | |
| 164 | self: self(type="RandomObject *") |
| 165 | |
| 166 | random() -> x in the interval [0, 1). |
| 167 | [clinic start generated code]*/ |
| 168 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 169 | static PyObject * |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 170 | _random_Random_random_impl(RandomObject *self) |
| 171 | /*[clinic end generated code: output=117ff99ee53d755c input=afb2a59cbbb00349]*/ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 172 | { |
Victor Stinner | 9f5fe79 | 2020-04-17 19:05:35 +0200 | [diff] [blame] | 173 | uint32_t a=genrand_uint32(self)>>5, b=genrand_uint32(self)>>6; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 174 | return PyFloat_FromDouble((a*67108864.0+b)*(1.0/9007199254740992.0)); |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 175 | } |
| 176 | |
| 177 | /* initializes mt[N] with a seed */ |
| 178 | static void |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 179 | init_genrand(RandomObject *self, uint32_t s) |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 180 | { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 181 | int mti; |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 182 | uint32_t *mt; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 183 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 184 | mt = self->state; |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 185 | mt[0]= s; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 186 | for (mti=1; mti<N; mti++) { |
| 187 | mt[mti] = |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 188 | (1812433253U * (mt[mti-1] ^ (mt[mti-1] >> 30)) + mti); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 189 | /* See Knuth TAOCP Vol2. 3rd Ed. P.106 for multiplier. */ |
| 190 | /* In the previous versions, MSBs of the seed affect */ |
| 191 | /* only MSBs of the array mt[]. */ |
| 192 | /* 2002/01/09 modified by Makoto Matsumoto */ |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 193 | } |
| 194 | self->index = mti; |
| 195 | return; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 196 | } |
| 197 | |
| 198 | /* initialize by an array with array-length */ |
| 199 | /* init_key is the array for initializing keys */ |
| 200 | /* key_length is its length */ |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 201 | static void |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 202 | init_by_array(RandomObject *self, uint32_t init_key[], size_t key_length) |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 203 | { |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 204 | size_t i, j, k; /* was signed in the original code. RDH 12/16/2002 */ |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 205 | uint32_t *mt; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 206 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 207 | mt = self->state; |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 208 | init_genrand(self, 19650218U); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 209 | i=1; j=0; |
| 210 | k = (N>key_length ? N : key_length); |
| 211 | for (; k; k--) { |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 212 | mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1664525U)) |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 213 | + init_key[j] + (uint32_t)j; /* non linear */ |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 214 | i++; j++; |
| 215 | if (i>=N) { mt[0] = mt[N-1]; i=1; } |
| 216 | if (j>=key_length) j=0; |
| 217 | } |
| 218 | for (k=N-1; k; k--) { |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 219 | mt[i] = (mt[i] ^ ((mt[i-1] ^ (mt[i-1] >> 30)) * 1566083941U)) |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 220 | - (uint32_t)i; /* non linear */ |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 221 | i++; |
| 222 | if (i>=N) { mt[0] = mt[N-1]; i=1; } |
| 223 | } |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 224 | |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 225 | mt[0] = 0x80000000U; /* MSB is 1; assuring non-zero initial array */ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 226 | } |
| 227 | |
| 228 | /* |
| 229 | * The rest is Python-specific code, neither part of, nor derived from, the |
| 230 | * Twister download. |
| 231 | */ |
| 232 | |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 233 | static int |
| 234 | random_seed_urandom(RandomObject *self) |
| 235 | { |
Victor Stinner | 1a1bd2e | 2020-04-17 19:13:06 +0200 | [diff] [blame] | 236 | uint32_t key[N]; |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 237 | |
| 238 | if (_PyOS_URandomNonblock(key, sizeof(key)) < 0) { |
| 239 | return -1; |
| 240 | } |
| 241 | init_by_array(self, key, Py_ARRAY_LENGTH(key)); |
| 242 | return 0; |
| 243 | } |
| 244 | |
| 245 | static void |
| 246 | random_seed_time_pid(RandomObject *self) |
| 247 | { |
| 248 | _PyTime_t now; |
| 249 | uint32_t key[5]; |
| 250 | |
| 251 | now = _PyTime_GetSystemClock(); |
Victor Stinner | 1a1bd2e | 2020-04-17 19:13:06 +0200 | [diff] [blame] | 252 | key[0] = (uint32_t)(now & 0xffffffffU); |
| 253 | key[1] = (uint32_t)(now >> 32); |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 254 | |
Victor Stinner | 1a1bd2e | 2020-04-17 19:13:06 +0200 | [diff] [blame] | 255 | key[2] = (uint32_t)getpid(); |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 256 | |
| 257 | now = _PyTime_GetMonotonicClock(); |
Victor Stinner | 1a1bd2e | 2020-04-17 19:13:06 +0200 | [diff] [blame] | 258 | key[3] = (uint32_t)(now & 0xffffffffU); |
| 259 | key[4] = (uint32_t)(now >> 32); |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 260 | |
| 261 | init_by_array(self, key, Py_ARRAY_LENGTH(key)); |
| 262 | } |
| 263 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 264 | static PyObject * |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 265 | random_seed(RandomObject *self, PyObject *arg) |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 266 | { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 267 | PyObject *result = NULL; /* guilty until proved innocent */ |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 268 | PyObject *n = NULL; |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 269 | uint32_t *key = NULL; |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 270 | size_t bits, keyused; |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 271 | int res; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 272 | |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 273 | if (arg == NULL || arg == Py_None) { |
| 274 | if (random_seed_urandom(self) < 0) { |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 275 | PyErr_Clear(); |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 276 | |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 277 | /* Reading system entropy failed, fall back on the worst entropy: |
| 278 | use the current time and process identifier. */ |
| 279 | random_seed_time_pid(self); |
| 280 | } |
| 281 | Py_RETURN_NONE; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 282 | } |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 283 | |
Larry Hastings | d60cd42 | 2012-06-24 02:52:21 -0700 | [diff] [blame] | 284 | /* This algorithm relies on the number being unsigned. |
| 285 | * So: if the arg is a PyLong, use its absolute value. |
| 286 | * Otherwise use its hash value, cast to unsigned. |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 287 | */ |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 288 | if (PyLong_CheckExact(arg)) { |
| 289 | n = PyNumber_Absolute(arg); |
| 290 | } else if (PyLong_Check(arg)) { |
Oren Milman | d780b2d | 2017-09-28 10:50:01 +0300 | [diff] [blame] | 291 | /* Calling int.__abs__() prevents calling arg.__abs__(), which might |
| 292 | return an invalid value. See issue #31478. */ |
Victor Stinner | 00d7cd8 | 2020-03-10 15:15:14 +0100 | [diff] [blame] | 293 | n = PyObject_CallOneArg(_randomstate_global->Long___abs__, arg); |
Oren Milman | d780b2d | 2017-09-28 10:50:01 +0300 | [diff] [blame] | 294 | } |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 295 | else { |
Larry Hastings | d60cd42 | 2012-06-24 02:52:21 -0700 | [diff] [blame] | 296 | Py_hash_t hash = PyObject_Hash(arg); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 297 | if (hash == -1) |
| 298 | goto Done; |
Larry Hastings | d60cd42 | 2012-06-24 02:52:21 -0700 | [diff] [blame] | 299 | n = PyLong_FromSize_t((size_t)hash); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 300 | } |
| 301 | if (n == NULL) |
| 302 | goto Done; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 303 | |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 304 | /* Now split n into 32-bit chunks, from the right. */ |
| 305 | bits = _PyLong_NumBits(n); |
| 306 | if (bits == (size_t)-1 && PyErr_Occurred()) |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 307 | goto Done; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 308 | |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 309 | /* Figure out how many 32-bit chunks this gives us. */ |
| 310 | keyused = bits == 0 ? 1 : (bits - 1) / 32 + 1; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 311 | |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 312 | /* Convert seed to byte sequence. */ |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 313 | key = (uint32_t *)PyMem_Malloc((size_t)4 * keyused); |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 314 | if (key == NULL) { |
Victor Stinner | a4ced86 | 2013-07-15 20:00:36 +0200 | [diff] [blame] | 315 | PyErr_NoMemory(); |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 316 | goto Done; |
Victor Stinner | a4ced86 | 2013-07-15 20:00:36 +0200 | [diff] [blame] | 317 | } |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 318 | res = _PyLong_AsByteArray((PyLongObject *)n, |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 319 | (unsigned char *)key, keyused * 4, |
| 320 | PY_LITTLE_ENDIAN, |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 321 | 0); /* unsigned */ |
| 322 | if (res == -1) { |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 323 | goto Done; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 324 | } |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 325 | |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 326 | #if PY_BIG_ENDIAN |
| 327 | { |
| 328 | size_t i, j; |
| 329 | /* Reverse an array. */ |
Zachary Ware | c15ea4c | 2015-05-17 23:46:22 -0500 | [diff] [blame] | 330 | for (i = 0, j = keyused - 1; i < j; i++, j--) { |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 331 | uint32_t tmp = key[i]; |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 332 | key[i] = key[j]; |
| 333 | key[j] = tmp; |
| 334 | } |
Mark Dickinson | 4cd6017 | 2012-12-21 21:52:49 +0000 | [diff] [blame] | 335 | } |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 336 | #endif |
Victor Stinner | e66987e | 2016-09-06 16:33:52 -0700 | [diff] [blame] | 337 | init_by_array(self, key, keyused); |
| 338 | |
| 339 | Py_INCREF(Py_None); |
| 340 | result = Py_None; |
| 341 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 342 | Done: |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 343 | Py_XDECREF(n); |
| 344 | PyMem_Free(key); |
| 345 | return result; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 346 | } |
| 347 | |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 348 | /*[clinic input] |
| 349 | _random.Random.seed |
| 350 | |
| 351 | self: self(type="RandomObject *") |
| 352 | n: object = None |
| 353 | / |
| 354 | |
| 355 | seed([n]) -> None. |
| 356 | |
| 357 | Defaults to use urandom and falls back to a combination |
| 358 | of the current time and the process identifier. |
| 359 | [clinic start generated code]*/ |
| 360 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 361 | static PyObject * |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 362 | _random_Random_seed_impl(RandomObject *self, PyObject *n) |
| 363 | /*[clinic end generated code: output=0fad1e16ba883681 input=78d6ef0d52532a54]*/ |
| 364 | { |
| 365 | return random_seed(self, n); |
| 366 | } |
| 367 | |
| 368 | /*[clinic input] |
| 369 | _random.Random.getstate |
| 370 | |
| 371 | self: self(type="RandomObject *") |
| 372 | |
| 373 | getstate() -> tuple containing the current state. |
| 374 | [clinic start generated code]*/ |
| 375 | |
| 376 | static PyObject * |
| 377 | _random_Random_getstate_impl(RandomObject *self) |
| 378 | /*[clinic end generated code: output=bf6cef0c092c7180 input=b937a487928c0e89]*/ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 379 | { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 380 | PyObject *state; |
| 381 | PyObject *element; |
| 382 | int i; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 383 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 384 | state = PyTuple_New(N+1); |
| 385 | if (state == NULL) |
| 386 | return NULL; |
| 387 | for (i=0; i<N ; i++) { |
| 388 | element = PyLong_FromUnsignedLong(self->state[i]); |
| 389 | if (element == NULL) |
| 390 | goto Fail; |
| 391 | PyTuple_SET_ITEM(state, i, element); |
| 392 | } |
| 393 | element = PyLong_FromLong((long)(self->index)); |
| 394 | if (element == NULL) |
| 395 | goto Fail; |
| 396 | PyTuple_SET_ITEM(state, i, element); |
| 397 | return state; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 398 | |
| 399 | Fail: |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 400 | Py_DECREF(state); |
| 401 | return NULL; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 402 | } |
| 403 | |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 404 | |
| 405 | /*[clinic input] |
| 406 | _random.Random.setstate |
| 407 | |
| 408 | self: self(type="RandomObject *") |
| 409 | state: object |
| 410 | / |
| 411 | |
| 412 | setstate(state) -> None. Restores generator state. |
| 413 | [clinic start generated code]*/ |
| 414 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 415 | static PyObject * |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 416 | _random_Random_setstate(RandomObject *self, PyObject *state) |
| 417 | /*[clinic end generated code: output=fd1c3cd0037b6681 input=b3b4efbb1bc66af8]*/ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 418 | { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 419 | int i; |
| 420 | unsigned long element; |
| 421 | long index; |
bladebryan | 9616a82 | 2017-04-21 23:10:46 -0700 | [diff] [blame] | 422 | uint32_t new_state[N]; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 423 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 424 | if (!PyTuple_Check(state)) { |
| 425 | PyErr_SetString(PyExc_TypeError, |
| 426 | "state vector must be a tuple"); |
| 427 | return NULL; |
| 428 | } |
| 429 | if (PyTuple_Size(state) != N+1) { |
| 430 | PyErr_SetString(PyExc_ValueError, |
| 431 | "state vector is the wrong size"); |
| 432 | return NULL; |
| 433 | } |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 434 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 435 | for (i=0; i<N ; i++) { |
| 436 | element = PyLong_AsUnsignedLong(PyTuple_GET_ITEM(state, i)); |
| 437 | if (element == (unsigned long)-1 && PyErr_Occurred()) |
| 438 | return NULL; |
bladebryan | 9616a82 | 2017-04-21 23:10:46 -0700 | [diff] [blame] | 439 | new_state[i] = (uint32_t)element; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 440 | } |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 441 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 442 | index = PyLong_AsLong(PyTuple_GET_ITEM(state, i)); |
| 443 | if (index == -1 && PyErr_Occurred()) |
| 444 | return NULL; |
Serhiy Storchaka | 178f0b6 | 2015-07-24 09:02:53 +0300 | [diff] [blame] | 445 | if (index < 0 || index > N) { |
| 446 | PyErr_SetString(PyExc_ValueError, "invalid state"); |
| 447 | return NULL; |
| 448 | } |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 449 | self->index = (int)index; |
bladebryan | 9616a82 | 2017-04-21 23:10:46 -0700 | [diff] [blame] | 450 | for (i = 0; i < N; i++) |
| 451 | self->state[i] = new_state[i]; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 452 | |
Serhiy Storchaka | 228b12e | 2017-01-23 09:47:21 +0200 | [diff] [blame] | 453 | Py_RETURN_NONE; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 454 | } |
| 455 | |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 456 | /*[clinic input] |
| 457 | |
| 458 | _random.Random.getrandbits |
| 459 | |
| 460 | self: self(type="RandomObject *") |
| 461 | k: int |
| 462 | / |
| 463 | |
| 464 | getrandbits(k) -> x. Generates an int with k random bits. |
| 465 | [clinic start generated code]*/ |
| 466 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 467 | static PyObject * |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 468 | _random_Random_getrandbits_impl(RandomObject *self, int k) |
| 469 | /*[clinic end generated code: output=b402f82a2158887f input=8c0e6396dd176fc0]*/ |
Raymond Hettinger | 2f726e9 | 2003-10-05 09:09:15 +0000 | [diff] [blame] | 470 | { |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 471 | int i, words; |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 472 | uint32_t r; |
| 473 | uint32_t *wordarray; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 474 | PyObject *result; |
Raymond Hettinger | 2f726e9 | 2003-10-05 09:09:15 +0000 | [diff] [blame] | 475 | |
Antoine Pitrou | 75a3378 | 2020-04-17 19:32:14 +0200 | [diff] [blame] | 476 | if (k < 0) { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 477 | PyErr_SetString(PyExc_ValueError, |
Antoine Pitrou | 75a3378 | 2020-04-17 19:32:14 +0200 | [diff] [blame] | 478 | "number of bits must be non-negative"); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 479 | return NULL; |
| 480 | } |
Raymond Hettinger | 2f726e9 | 2003-10-05 09:09:15 +0000 | [diff] [blame] | 481 | |
Antoine Pitrou | 75a3378 | 2020-04-17 19:32:14 +0200 | [diff] [blame] | 482 | if (k == 0) |
| 483 | return PyLong_FromLong(0); |
| 484 | |
Serhiy Storchaka | d8a0bac | 2013-01-04 12:18:35 +0200 | [diff] [blame] | 485 | if (k <= 32) /* Fast path */ |
Victor Stinner | 9f5fe79 | 2020-04-17 19:05:35 +0200 | [diff] [blame] | 486 | return PyLong_FromUnsignedLong(genrand_uint32(self) >> (32 - k)); |
Serhiy Storchaka | d8a0bac | 2013-01-04 12:18:35 +0200 | [diff] [blame] | 487 | |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 488 | words = (k - 1) / 32 + 1; |
Benjamin Peterson | 9b3d770 | 2016-09-06 13:24:00 -0700 | [diff] [blame] | 489 | wordarray = (uint32_t *)PyMem_Malloc(words * 4); |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 490 | if (wordarray == NULL) { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 491 | PyErr_NoMemory(); |
| 492 | return NULL; |
| 493 | } |
Raymond Hettinger | 2f726e9 | 2003-10-05 09:09:15 +0000 | [diff] [blame] | 494 | |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 495 | /* Fill-out bits of long integer, by 32-bit words, from least significant |
| 496 | to most significant. */ |
| 497 | #if PY_LITTLE_ENDIAN |
| 498 | for (i = 0; i < words; i++, k -= 32) |
| 499 | #else |
| 500 | for (i = words - 1; i >= 0; i--, k -= 32) |
| 501 | #endif |
| 502 | { |
Victor Stinner | 9f5fe79 | 2020-04-17 19:05:35 +0200 | [diff] [blame] | 503 | r = genrand_uint32(self); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 504 | if (k < 32) |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 505 | r >>= (32 - k); /* Drop least significant bits */ |
| 506 | wordarray[i] = r; |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 507 | } |
Raymond Hettinger | 2f726e9 | 2003-10-05 09:09:15 +0000 | [diff] [blame] | 508 | |
Serhiy Storchaka | dce0405 | 2015-05-13 15:02:12 +0300 | [diff] [blame] | 509 | result = _PyLong_FromByteArray((unsigned char *)wordarray, words * 4, |
| 510 | PY_LITTLE_ENDIAN, 0 /* unsigned */); |
| 511 | PyMem_Free(wordarray); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 512 | return result; |
Raymond Hettinger | 2f726e9 | 2003-10-05 09:09:15 +0000 | [diff] [blame] | 513 | } |
| 514 | |
| 515 | static PyObject * |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 516 | random_new(PyTypeObject *type, PyObject *args, PyObject *kwds) |
| 517 | { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 518 | RandomObject *self; |
| 519 | PyObject *tmp; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 520 | |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 521 | if (type == (PyTypeObject*)_randomstate_global->Random_Type && |
| 522 | !_PyArg_NoKeywords("Random()", kwds)) { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 523 | return NULL; |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 524 | } |
Georg Brandl | 02c4287 | 2005-08-26 06:42:30 +0000 | [diff] [blame] | 525 | |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 526 | self = (RandomObject *)PyType_GenericAlloc(type, 0); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 527 | if (self == NULL) |
| 528 | return NULL; |
| 529 | tmp = random_seed(self, args); |
| 530 | if (tmp == NULL) { |
| 531 | Py_DECREF(self); |
| 532 | return NULL; |
| 533 | } |
| 534 | Py_DECREF(tmp); |
| 535 | return (PyObject *)self; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 536 | } |
| 537 | |
| 538 | static PyMethodDef random_methods[] = { |
Pablo Galindo | 561612d | 2019-05-24 22:09:23 +0100 | [diff] [blame] | 539 | _RANDOM_RANDOM_RANDOM_METHODDEF |
| 540 | _RANDOM_RANDOM_SEED_METHODDEF |
| 541 | _RANDOM_RANDOM_GETSTATE_METHODDEF |
| 542 | _RANDOM_RANDOM_SETSTATE_METHODDEF |
| 543 | _RANDOM_RANDOM_GETRANDBITS_METHODDEF |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 544 | {NULL, NULL} /* sentinel */ |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 545 | }; |
| 546 | |
| 547 | PyDoc_STRVAR(random_doc, |
| 548 | "Random() -> create a random number generator with its own internal state."); |
| 549 | |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 550 | static PyType_Slot Random_Type_slots[] = { |
Victor Stinner | 2f90261 | 2019-10-01 12:45:52 +0200 | [diff] [blame] | 551 | {Py_tp_doc, (void *)random_doc}, |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 552 | {Py_tp_methods, random_methods}, |
| 553 | {Py_tp_new, random_new}, |
| 554 | {Py_tp_free, PyObject_Free}, |
| 555 | {0, 0}, |
| 556 | }; |
| 557 | |
| 558 | static PyType_Spec Random_Type_spec = { |
| 559 | "_random.Random", |
| 560 | sizeof(RandomObject), |
| 561 | 0, |
| 562 | Py_TPFLAGS_DEFAULT | Py_TPFLAGS_BASETYPE, |
| 563 | Random_Type_slots |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 564 | }; |
| 565 | |
| 566 | PyDoc_STRVAR(module_doc, |
| 567 | "Module implements the Mersenne Twister random number generator."); |
| 568 | |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 569 | static int |
| 570 | _random_traverse(PyObject *module, visitproc visit, void *arg) |
| 571 | { |
Hai Shi | f707d94 | 2020-03-16 21:15:01 +0800 | [diff] [blame] | 572 | Py_VISIT(get_random_state(module)->Random_Type); |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 573 | return 0; |
| 574 | } |
| 575 | |
| 576 | static int |
| 577 | _random_clear(PyObject *module) |
| 578 | { |
Hai Shi | f707d94 | 2020-03-16 21:15:01 +0800 | [diff] [blame] | 579 | Py_CLEAR(get_random_state(module)->Random_Type); |
| 580 | Py_CLEAR(get_random_state(module)->Long___abs__); |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 581 | return 0; |
| 582 | } |
| 583 | |
| 584 | static void |
| 585 | _random_free(void *module) |
| 586 | { |
| 587 | _random_clear((PyObject *)module); |
| 588 | } |
Martin v. Löwis | 1a21451 | 2008-06-11 05:26:20 +0000 | [diff] [blame] | 589 | |
| 590 | static struct PyModuleDef _randommodule = { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 591 | PyModuleDef_HEAD_INIT, |
| 592 | "_random", |
| 593 | module_doc, |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 594 | sizeof(_randomstate), |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 595 | NULL, |
| 596 | NULL, |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 597 | _random_traverse, |
| 598 | _random_clear, |
| 599 | _random_free, |
Martin v. Löwis | 1a21451 | 2008-06-11 05:26:20 +0000 | [diff] [blame] | 600 | }; |
| 601 | |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 602 | PyMODINIT_FUNC |
Martin v. Löwis | 1a21451 | 2008-06-11 05:26:20 +0000 | [diff] [blame] | 603 | PyInit__random(void) |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 604 | { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 605 | PyObject *m; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 606 | |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 607 | PyObject *Random_Type = PyType_FromSpec(&Random_Type_spec); |
| 608 | if (Random_Type == NULL) { |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 609 | return NULL; |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 610 | } |
| 611 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 612 | m = PyModule_Create(&_randommodule); |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 613 | if (m == NULL) { |
| 614 | Py_DECREF(Random_Type); |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 615 | return NULL; |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 616 | } |
Hai Shi | f707d94 | 2020-03-16 21:15:01 +0800 | [diff] [blame] | 617 | get_random_state(m)->Random_Type = Random_Type; |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 618 | |
| 619 | Py_INCREF(Random_Type); |
| 620 | PyModule_AddObject(m, "Random", Random_Type); |
| 621 | |
| 622 | /* Look up and save int.__abs__, which is needed in random_seed(). */ |
| 623 | PyObject *longval = NULL, *longtype = NULL; |
| 624 | longval = PyLong_FromLong(0); |
| 625 | if (longval == NULL) goto fail; |
| 626 | |
| 627 | longtype = PyObject_Type(longval); |
| 628 | if (longtype == NULL) goto fail; |
| 629 | |
| 630 | PyObject *abs = PyObject_GetAttrString(longtype, "__abs__"); |
| 631 | if (abs == NULL) goto fail; |
| 632 | |
| 633 | Py_DECREF(longtype); |
| 634 | Py_DECREF(longval); |
Hai Shi | f707d94 | 2020-03-16 21:15:01 +0800 | [diff] [blame] | 635 | get_random_state(m)->Long___abs__ = abs; |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 636 | |
Antoine Pitrou | f95a1b3 | 2010-05-09 15:52:27 +0000 | [diff] [blame] | 637 | return m; |
Dino Viehland | 04f0bbf | 2019-09-13 11:12:27 +0100 | [diff] [blame] | 638 | |
| 639 | fail: |
| 640 | Py_XDECREF(longtype); |
| 641 | Py_XDECREF(longval); |
| 642 | Py_DECREF(m); |
| 643 | return NULL; |
Raymond Hettinger | 40f6217 | 2002-12-29 23:03:38 +0000 | [diff] [blame] | 644 | } |