Displaying 2 results from an estimated 2 matches for "marsaglia_multicarry".
2003 Jan 09
1
RNG.c: unif_rand, MARSAGLIA_MULTICARRY (PR#2437)
Full_Name: Richard Simard
Version:
OS: Linux
Submission from: (NULL) (132.204.25.139)
In the random number generator MARSAGLIA_MULTICARRY, your algorithm is
different than the one in the original reference that you give in the R manual:
Marsaglia in his post to the mailing list {\it sci.stat.math} on September 29,
1997.
The last line in the R program has a ^
while Marsaglia's algorithm has a +.
This gives very different random nu...
2007 Sep 23
0
initial scrambling of seed in do_setseed / RNG_Init
...= 0.0; /* zap Box-Muller history */
/* Initial scrambling */
- for(j = 0; j < 50; j++)
- seed = (69069 * seed + 1);
+ /* Equivalent to 50 iterations of seed = (69069 * seed + 1) */
+ seed = (1100682473U * seed + 2358491998U);
switch(kind) {
case WICHMANN_HILL:
case MARSAGLIA_MULTICARRY:
Let me add also the following remark. The help(Random.user) page says
Optionally, the user can supply the entry point 'user_unif_init',
which is called with an 'unsigned int' argument when 'RNGkind' (or
'set.seed') is called, and is intended to be u...