ref: 26f0aa33d91442fc7a56fc5f5c9b56b5f5f84ce2
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; }