git: 9front

ref: c2a5d7e3f28e86aeef1e4e37e2b88e46af90b62c
dir: /sys/src/libc/port/frand.c/

View raw version 
#include	<u.h>
#include	<libc.h>

#define	MASK	0x7fffffffL
#define	NORM	(1.0/(1.0+MASK))

double
frand(void)
{
	double x;

	do {
		x = lrand() * NORM;
		x = (x + lrand()) * NORM;
	} while(x >= 1);
	return x;
}