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