naginterfaces.library.rand.init_​nonrepeat

naginterfaces.library.rand.init_nonrepeat(genid, subid=0)

init_nonrepeat initializes the selected base generator to generate a non-repeatable sequence of variates. The base generator can then be used by the group of pseudorandom number functions (see init_leapfrog()–init_skipahead(), permute(), sample(), times_garch_asym1()–times_mv_varma(), matrix_orthog()–matrix_2waytable(), copula_students_t(), copula_normal(), multivar_students_t(), multivar_normal() and dist_uniform01()–int_uniform()) and the quasi-random scrambled sequence initialization function, quasi_init_scrambled().

For full information please refer to the NAG Library document for g05kg

https://www.nag.com/numeric/nl/nagdoc_29.3/flhtml/g05/g05kgf.html

Parameters

genidint

Must contain the type of base generator to use.

$genid=1$

NAG basic generator.

$genid=2$

Wichmann Hill I generator.

$genid=3$

Mersenne Twister.

$genid=4$

Wichmann Hill II generator.

$genid=5$

ACORN generator.

$genid=6$

L’Ecuyer MRG32k3a generator.

See the G05 Introduction for details of each of the base generators.

subidint, optional

If $genid=2$, $subid$ indicates which of the $273$ sub-generators to use. In this case, the $\left((\left|subid\right|+272)mod273\right)+1$ sub-generator is used.

If $genid=5$, $subid$ indicates the values of $k$ and $p$ to use, where $k$ is the order of the generator, and $p$ controls the size of the modulus, $M$, with $M=2^{(p\times 30)}$.

If $subid<1$, the default values of $k=10$ and $p=2$ are used, otherwise values for $k$ and $p$ are calculated from the formula, $subid=k+1000(p-1)$.

If $genid=6$ and $mod(subid,2)=0$ the range of the generator is set to $(0,1]$, otherwise the range is set to $(0,1)$; in this case the sequence is identical to the implementation of MRG32k3a in TestU01 (see L’Ecuyer and Simard (2002)) for identical seeds.

For all other values of $genid$, $subid$ is not referenced.

Returns

statecommdict, RNG communication object

RNG communication structure.

Raises

NagValueError

(errno $1$)

On entry, $genid=⟨value⟩$.

Constraint: $genid=1$, $2$, $3$, $4$, $5$ or $6$.

Notes

init_nonrepeat selects a base generator through the input value of the arguments $genid$ and $subid$, and then initializes it based on the values taken from the real-time clock, resulting in the same base generator yielding different sequences of random numbers each time the calling program is run. It should be noted that there is no guarantee of statistical properties between sequences, only within sequences.

A definition of some of the terms used in this description, along with details of the various base generators can be found in the G05 Introduction.

References

L’Ecuyer, P and Simard, R, 2002, TestU01: a software library in ANSI C for empirical testing of random number generators, Departement d’Informatique et de Recherche Operationnelle, Universite de Montreal, https://www.iro.umontreal.ca/~lecuyer

Maclaren, N M, 1989, The generation of multiple independent sequences of pseudorandom numbers, Appl. Statist. (38), 351–359

Matsumoto, M and Nishimura, T, 1998, Mersenne twister: a 623-dimensionally equidistributed uniform pseudorandom number generator, ACM Transactions on Modelling and Computer Simulations

Wichmann, B A and Hill, I D, 2006, Generating good pseudo-random numbers, Computational Statistics and Data Analysis (51), 1614–1622

Wikramaratna, R S, 1989, ACORN - a new method for generating sequences of uniformly distributed pseudo-random numbers, Journal of Computational Physics (83), 16–31