naginterfaces.library.rand.init_​leapfrog

naginterfaces.library.rand.init_leapfrog(n, k, statecomm)

init_leapfrog allows for the generation of multiple, independent, sequences of pseudorandom numbers using the leap-frog method.

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

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

Parameters

nint

$n$, the total number of sequences required.

kint

$k$, the number of the current sequence.

statecommdict, RNG communication object, modified in place

RNG communication structure.

This argument must have been initialized by a prior call to init_repeat() or init_nonrepeat().

Raises

NagValueError

(errno $1$)

On entry, $n=⟨value⟩$.

Constraint: $n\ge 1$.

(errno $2$)

On entry, $k=⟨value⟩$ and $n=⟨value⟩$.

Constraint: $0<k\le n$.

(errno $3$)

On entry, $statecomm$[‘state’] vector has been corrupted or not initialized.

(errno $4$)

On entry, cannot use leap-frog with the base generator defined by $statecomm$[‘state’].

Notes

init_leapfrog adjusts a base generator to allow multiple, independent, sequences of pseudorandom numbers to be generated via the leap-frog method (see the G05 Introduction for details).

If, prior to calling init_leapfrog the base generator defined by $statecomm$[‘state’] would produce random numbers $x_1,x_2,x_3,\dots$, then after calling init_leapfrog the generator will produce random numbers $x_k,x_{k+n},x_{k+2n},x_{k+3n},\dots$.

One of the initialization functions init_repeat() (for a repeatable sequence if computed sequentially) or init_nonrepeat() (for a non-repeatable sequence) must be called prior to the first call to init_leapfrog.

The leap-frog algorithm can be used in conjunction with the NAG basic generator, both the Wichmann–Hill I and Wichmann–Hill II generators, the Mersenne Twister and L’Ecuyer.

References

Knuth, D E, 1981, The Art of Computer Programming (Volume 2), (2nd Edition), Addison–Wesley