ref: 1132d1b9df4c9ea2b857fa9778fe6762cd81ded8
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; }