naginterfaces.library.rand.quasi_​normal

naginterfaces.library.rand.quasi_normal(xmean, std, n, comm)

quasi_normal generates a quasi-random sequence from a Normal (Gaussian) distribution. It must be preceded by a call to one of the initialization functions quasi_init() or quasi_init_scrambled().

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

https://support.nag.com/numeric/nl/nagdoc_30.3/flhtml/g05/g05yjf.html

Parameters

xmeanfloat, array-like, shape $\left(\text{idim}\right)$

Specifies, for each dimension, the mean of the Normal distribution.

stdfloat, array-like, shape $\left(\text{idim}\right)$

Specifies, for each dimension, the standard deviation of the Normal distribution.

nint

The number of quasi-random numbers required.

commdict, communication object, modified in place

Communication structure.

This argument must have been initialized by a prior call to quasi_init() or quasi_init_scrambled().

Returns

quasfloat, ndarray, shape $(n,\text{idim})$

Contains the $n$ quasi-random numbers of dimension idim

If $\text{sorder}=2$, $quas[i-1,j-1]$ holds the $j$th value for the $i$th dimension.

If $\text{sorder}=1$, $quas[i-1,j-1]$ holds the $i$th value for the $j$th dimension.

Raises

NagValueError

(errno $1$)

On entry, $comm$[‘iref’] has either not been initialized or has been corrupted.

(errno $2$)

On entry, $n=⟨value⟩$.

Constraint: $n\ge 0$.

(errno $3$)

On entry, $i=⟨value⟩$ and $std[i-1]=⟨value⟩$.

Constraint: $std[i-1]\ge 0.0$.

(errno $4$)

There have been too many calls to the generator.

Notes

quasi_normal generates a quasi-random sequence from a Normal distribution by first generating a uniform quasi-random sequence which is then transformed into a Normal sequence using the inverse of the Normal CDF. The type of uniform sequence used depends on the initialization function called and can include the low-discrepancy sequences proposed by Sobol, Faure or Niederreiter. If the initialization function quasi_init_scrambled() was used then the underlying uniform sequence is first scrambled prior to being transformed (see Notes for quasi_init_scrambled for details).

References

Bratley, P and Fox, B L, 1988, Algorithm 659: implementing Sobol’s quasirandom sequence generator, ACM Trans. Math. Software (14(1)), 88–100

Fox, B L, 1986, Algorithm 647: implementation and relative efficiency of quasirandom sequence generators, ACM Trans. Math. Software (12(4)), 362–376

Wichura, 1988, Algorithm AS 241: the percentage points of the Normal distribution, Appl. Statist. (37), 477–484

See also

naginterfaces.library.examples.rand.bb_inc_ex.main()