| #include <stdint.h> | |
| uint64_t __rand48_step(unsigned short *xi, unsigned short *lc) | |
| { | |
| uint64_t a, x; | |
| x = xi[0] | ((uint32_t)xi[1])<<16 | ((uint64_t)xi[2])<<32; | |
| a = lc[0] | ((uint32_t)lc[1])<<16 | ((uint64_t)lc[2])<<32; | |
| x = a*x + lc[3]; | |
| xi[0] = x; | |
| xi[1] = x>>16; | |
| xi[2] = x>>32; | |
| return x & 0xffffffffffffull; | |
| } |