This is a
simple random number generator providing deviates in the range 0 to 1, with period of
and to 6 or 7 digits accuracy. It is based upon Ahrens, Dieter & Grube
’s TOMS599 routines. Note
that there is no STATUS argument for efficiency.
RESULT = PDA_RAND( X
X = REAL (Given)
This is a dummy variable required by the Fortran
PDA_RAND = REAL
The pseudo-random deviate.
The initial seed MUST be set using routine PDA_RNSED (equivalent to NAG
G05CBF). If it has not, there is no guarantee that sensible values will be returned from
Ahrens, J.H., Dieter, U. & Grube, A., 1970,
numbers: a new proposal for the choice of multiplicators
", Computing, 6, pp.121–138.