ref: 0e54aa2a4df49f9a4fb7f1a275940332d3bee3c6
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; }