| 112 | } |
| 113 | |
| 114 | static unsigned int rnd_int( rnd *r ) |
| 115 | { |
| 116 | unsigned int y; |
| 117 | int pos = r->cur++; |
| 118 | if( pos >= NSEEDS ) { |
| 119 | int kk; |
| 120 | for(kk=0;kk<NSEEDS-MAX;kk++) |
| 121 | r->seeds[kk] = r->seeds[kk+MAX] ^ (r->seeds[kk] >> 1) ^ mag01[r->seeds[kk] % 2]; |
| 122 | for(;kk<NSEEDS;kk++) |
| 123 | r->seeds[kk] = r->seeds[kk+(MAX-NSEEDS)] ^ (r->seeds[kk] >> 1) ^ mag01[r->seeds[kk] % 2]; |
| 124 | r->cur = 1; |
| 125 | pos = 0; |
| 126 | } |
| 127 | y = r->seeds[pos]; |
| 128 | y ^= (y << 7) & 0x2b5b2500; |
| 129 | y ^= (y << 15) & 0xdb8b0000; |
| 130 | y ^= (y >> 16); |
| 131 | return y; |
| 132 | } |
| 133 | |
| 134 | static double rnd_float( rnd *r ) |
| 135 | { |
no outgoing calls
no test coverage detected