search for: congr_mm

Displaying 1 result from an estimated 1 matches for "congr_mm".

1997 Nov 24
0
R-alpha: random number generator -- S-plus's
...tosh <mcintosh@larch.bellcore.com> X + a correction of mine (B.R. wrote '69096' instead of '69069')) X X --> /u/maechler/tex/o.peopl/Ripley-S-RNG.tex */ X /*--- This code is perfectly tested out via S-code X X in --> /u/maechler/S/MISC/runif-tst.S */ X unsigned long congr_MM, taus_MM; X static double xuni() { X unsigned long n, lambda = 69069; X congr_MM *= lambda; X taus_MM ^= taus_MM >> 15; X taus_MM ^= taus_MM << 17; /* Orig (Splus <= 3.2): X n = taus_MM ^ congr_MM; X return ( (n>>1) / 2147483648.); X...