69 #include <sys/types.h>
78 # if !defined(_WIN32_WINNT) || _WIN32_WINNT < 0x0400
80 # define _WIN32_WINNT 0x400
81 # undef __WINCRYPT_H__
91 #define MATRIX_A 0x9908b0dfU
92 #define UMASK 0x80000000U
93 #define LMASK 0x7fffffffU
94 #define MIXBITS(u,v) ( ((u) & UMASK) | ((v) & LMASK) )
95 #define TWIST(u,v) ((MIXBITS((u),(v)) >> 1) ^ ((v)&1U ? MATRIX_A : 0U))
106 #define genrand_initialized(mt) ((mt)->next != 0)
107 #define uninit_genrand(mt) ((mt)->next = 0)
114 mt->
state[0] = s & 0xffffffff
U;
115 for (j=1; j<
N; j++) {
116 mt->
state[j] = (1812433253U * (mt->
state[j-1] ^ (mt->
state[j-1] >> 30)) + j);
121 mt->
state[j] &= 0xffffffff;
137 k = (
N>key_length ?
N : key_length);
144 if (j>=key_length) j=0;
146 for (k=
N-1; k; k--) {
154 mt->
state[0] = 0x80000000
U;
160 unsigned int *
p = mt->
state;
166 for (j=
N-
M+1; --j; p++)
167 *p = p[
M] ^
TWIST(p[0], p[1]);
170 *p = p[
M-
N] ^
TWIST(p[0], p[1]);
187 y ^= (y << 7) & 0x9d2c5680;
188 y ^= (y << 15) & 0xefc60000;
200 return(a*67108864.0+b)*(1.0/9007199254740992.0);
225 #define DEFAULT_SEED_CNT 4
235 struct MT *mt = &r->
mt;
262 #define BDIGITS(x) (RBIGNUM_DIGITS(x))
263 #define BITSPERDIG (SIZEOF_BDIGITS*CHAR_BIT)
264 #define BIGRAD ((BDIGIT_DBL)1 << BITSPERDIG)
265 #define DIGSPERINT (SIZEOF_INT/SIZEOF_BDIGITS)
266 #define BIGUP(x) ((BDIGIT_DBL)(x) << BITSPERDIG)
267 #define BIGDN(x) RSHIFT((x),BITSPERDIG)
268 #define BIGLO(x) ((BDIGIT)((x) & (BIGRAD-1)))
269 #define BDIGMAX ((BDIGIT)-1)
271 #define roomof(n, m) (int)(((n)+(m)-1) / (m))
272 #define numberof(array) (int)(sizeof(array) / sizeof((array)[0]))
273 #define SIZEOF_INT32 (31/CHAR_BIT + 1)
300 #if CHAR_BIT * SIZEOF_LONG > 64
307 #if 64 % BITSPERDIG == 0
323 return ldexp(r, -53);
342 if (ptr != &default_rand)
399 switch (
TYPE(seed)) {
404 fixnum_seed = -fixnum_seed;
405 buf[0] = (
unsigned int)(fixnum_seed & 0xffffffff);
406 #if SIZEOF_LONG > SIZEOF_INT32
407 if ((
long)(int32_t)fixnum_seed != fixnum_seed) {
408 if ((buf[1] = (
unsigned int)(fixnum_seed >> 32)) != 0) ++len;
424 memset(buf, 0, len *
sizeof(*buf));
426 for (i = (
int)(blen-1); 0 <=
i; i--) {
428 #if SIZEOF_BDIGITS < SIZEOF_INT32
432 if (!len && buf[j]) len = j;
448 if (buf != buf0)
xfree(buf);
484 #define DEFAULT_SEED_LEN (DEFAULT_SEED_CNT * (int)sizeof(int))
486 #if defined(S_ISCHR) && !defined(DOSISH)
487 # define USE_DEV_URANDOM 1
489 # define USE_DEV_URANDOM 0
500 #elif defined(_WIN32)
507 if ((fd = open(
"/dev/urandom", O_RDONLY
516 if (
fstat(fd, &statbuf) == 0 &&
S_ISCHR(statbuf.st_mode)) {
523 #elif defined(_WIN32)
524 if (CryptAcquireContext(&prov,
NULL,
NULL, PROV_RSA_FULL, CRYPT_VERIFYCONTEXT)) {
526 CryptReleaseContext(prov, 0);
532 seed[1] ^= (
unsigned int)tv.
tv_sec;
533 #
if SIZEOF_TIME_T > SIZEOF_INT
536 seed[2] ^= getpid() ^ (n++ << 16);
537 seed[3] ^= (
unsigned int)(
VALUE)&seed;
538 #if SIZEOF_VOIDP > SIZEOF_INT
539 seed[2] ^= (
unsigned int)((
VALUE)&seed >> SIZEOF_INT *
CHAR_BIT);
559 #if SIZEOF_INT32 / SIZEOF_BDIGITS > 1
560 digits[len-2] <= 1 && digits[len-1] == 0
599 struct MT *mt = &rnd1->
mt;
614 unsigned int x = mt->
state[
i];
615 #if SIZEOF_BDIGITS < SIZEOF_INT32
677 struct MT *mt = &rnd->
mt;
698 mt->
state[0] = (
unsigned int)x;
699 #if SIZEOF_LONG / SIZEOF_INT >= 2
702 #if SIZEOF_LONG / SIZEOF_INT >= 3
705 #if SIZEOF_LONG / SIZEOF_INT >= 4
717 #if SIZEOF_BDIGITS < SIZEOF_INT
727 }
while (--len % DIGSPERINT);
740 # elif SIZEOF_BDIGITS < SIZEOF_INT
743 }
while (--len % DIGSPERINT);
753 mt->
left = (
unsigned int)x;
812 unsigned long val, mask;
814 if (!limit)
return 0;
819 if ((mask >> (i * 32)) & 0xffffffff) {
833 unsigned long mask, lim, rnd;
841 #if SIZEOF_BDIGITS == 2
842 # define BIG_GET32(big,i) \
843 (RBIGNUM_DIGITS(big)[(i)*2] | \
844 ((i)*2+1 < RBIGNUM_LEN(big) ? \
845 (RBIGNUM_DIGITS(big)[(i)*2+1] << 16) : \
847 # define BIG_SET32(big,i,d) \
848 ((RBIGNUM_DIGITS(big)[(i)*2] = (d) & 0xffff), \
849 ((i)*2+1 < RBIGNUM_LEN(big) ? \
850 (RBIGNUM_DIGITS(big)[(i)*2+1] = (d) >> 16) : \
854 # define BIG_GET32(big,i) (RBIGNUM_DIGITS(big)[(i)])
855 # define BIG_SET32(big,i,d) (RBIGNUM_DIGITS(big)[(i)] = (d))
860 for (i = len-1; 0 <=
i; i--) {
862 mask = mask ? 0xffffffff :
make_mask(lim);
897 #if SIZEOF_LONG * CHAR_BIT > 32
899 #elif defined HAVE_LONG_LONG
900 VALUE lim = ULL2NUM((LONG_LONG)0xffffffff+1);
916 if (d < 0.0 || d >= 1.0) {
974 if (endp) *endp = end;
990 if (!max)
return Qnil;
992 if (restrictive)
return Qnil;
1002 if (restrictive)
return Qnil;
1009 if (max == -1)
return Qnil;
1019 static inline double
1044 if ((max =
FIX2LONG(vmax) - excl) >= 0) {
1045 unsigned long r =
limited_rand(mt, (
unsigned long)max);
1080 return rb_float_new(+(+(+(r - 0.5) * max) * scale) + mid);
1084 else if (max == 0.0 && !excl) {
1145 else if (argc != 1) {
1232 if (argc == 0)
goto zero_arg;
1234 if (
NIL_P(vmax))
goto zero_arg;
1262 #define SIP_HASH_STREAMING 0
1263 #define sip_hash24 ruby_sip_hash24
1264 #if !defined _WIN32 && !defined BYTE_ORDER
1265 # ifdef WORDS_BIGENDIAN
1266 # define BYTE_ORDER BIG_ENDIAN
1268 # define BYTE_ORDER LITTLE_ENDIAN
1270 # ifndef LITTLE_ENDIAN
1271 # define LITTLE_ENDIAN 1234
1274 # define BIG_ENDIAN 4321
1301 struct MT *mt = &r->
mt;
1306 #if SIZEOF_ST_INDEX_T*CHAR_BIT > 4*8
1310 #if SIZEOF_ST_INDEX_T*CHAR_BIT > 8*8
1314 #if SIZEOF_ST_INDEX_T*CHAR_BIT > 12*8
1336 #ifdef HAVE_UINT64_T