summaryrefslogtreecommitdiff
path: root/src/prng/drand48.c
blob: 08283e240b14fad45921469bc7aea39eadda3670 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
#include <stdlib.h>
#include <inttypes.h>
#include "rand48.h"

double erand48(unsigned short s[3])
{
	union {
		uint64_t u;
		double f;
	} x = { 0x3ff0000000000000ULL | __rand48_step(s, __seed48+3)<<4 };
	return x.f - 1.0;
}

double drand48(void)
{
	return erand48(__seed48);
}