Forgotten functions in the graph adaptor extenders.
3 #include <lemon/bits/mingw32_rand.h>
6 static unsigned short _mingw_rand_state[10];
8 void _mingw_rand_next_state(unsigned short xsubi[3]) {
9 _mingw_rand_state[8] = _mingw_rand_state[9] = 0;
10 _mingw_rand_state[7] = _mingw_rand_state[6];
11 for (int i = 0; i < 3; ++i) {
12 unsigned long val = 0;
13 for (int j = 0; i + j < 3; ++j) {
14 val += (unsigned long)_mingw_rand_state[7 + i + j]
15 + (unsigned long)xsubi[i] * (unsigned long)_mingw_rand_state[3 + j];
16 _mingw_rand_state[7 + i + j] = val;
20 xsubi[0] = _mingw_rand_state[6];
21 xsubi[1] = _mingw_rand_state[7];
22 xsubi[2] = _mingw_rand_state[8];
25 long int lrand48(void) {
26 return nrand48(_mingw_rand_state);
29 long int nrand48(unsigned short xsubi[3]) {
30 _mingw_rand_next_state(xsubi);
31 return ((long)(xsubi[2] & ( (1 << 15) - 1) ) << 16) | (long)xsubi[1];
34 double drand48(void) {
35 return erand48(_mingw_rand_state);
38 double erand48(unsigned short xsubi[3]) {
39 return (double)nrand48(xsubi) / (1 << 31);
43 long int mrand48(void) {
44 return jrand48(_mingw_rand_state);
47 long int jrand48(unsigned short xsubi[3]) {
48 _mingw_rand_next_state(xsubi);
49 return ((long)xsubi[2] << 16) | (long)xsubi[1];
52 void srand48(long int seedval) {
53 _mingw_rand_state[0] = 0x330E;
54 _mingw_rand_state[1] = seedval & ( (1 << 16) - 1);
55 _mingw_rand_state[2] = seedval >> 16;
56 _mingw_rand_state[3] = 0xE66D;
57 _mingw_rand_state[4] = 0xDEEC;
58 _mingw_rand_state[5] = 0x0005;
59 _mingw_rand_state[6] = 0x000B;
62 unsigned short *seed48(unsigned short seed16v[3]) {
63 _mingw_rand_state[7] = _mingw_rand_state[0];
64 _mingw_rand_state[8] = _mingw_rand_state[1];
65 _mingw_rand_state[9] = _mingw_rand_state[2];
66 _mingw_rand_state[0] = seed16v[0];
67 _mingw_rand_state[1] = seed16v[1];
68 _mingw_rand_state[2] = seed16v[2];
69 return _mingw_rand_state + 7;
72 void lcong48(unsigned short param[7]) {
73 _mingw_rand_state[0] = param[0];
74 _mingw_rand_state[1] = param[1];
75 _mingw_rand_state[2] = param[2];
76 _mingw_rand_state[3] = param[3];
77 _mingw_rand_state[4] = param[4];
78 _mingw_rand_state[5] = param[5];
79 _mingw_rand_state[6] = param[6];