15 lines
305 B
C
15 lines
305 B
C
#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;
|
|
}
|