ref: 8f4053d1efa1ee75a7bea09e624f1134d71a289a
dir: /sys/src/libc/port/frand.c/
#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; }