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...