summaryrefslogblamecommitdiff
path: root/src/prng/__rand48_step.c
blob: 105f810ff017d51e8979b049bae0167ec612a497 (plain) (tree)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
                   
                   











                                                              
#include <stdlib.h>
#include <stdint.h>

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