Up \SpecFunc
unit uRand13;
{@abstract(13 Dimensional Low Discrepancy Sobol Sequences
This unit realizes the Lambda/Pi/Tau sequence by V.M. Sobol.
The Numerical Monte-Carlo Methods, Moscow, Mir, 1979.
The present version was written by E.B. Krissinel', (C) 1991.)
author: Eugene B. Krissinel http://www.fh.huji.ac.il/~krissinel/)
created: February 02, 1991)
last modified: February 02, 2003) }
interface
uses
uMatTypes;
const
Rand13Dim = 13;
type
Rand13Point = array [1..Rand13Dim] of TFloat;
var
Rand13Seed : longint;
{ This is the current Random Generator's state variable;
It may be saved and restored;
The initial value should be assigned by Randomize13. }
{ Call for the random initial point installation.}
procedure Randomize13;
procedure Rand13(M: integer; var Q : Rand13Point );
{ M is the dimension of random-point space: 1 <= M <= Rand13Dim,
Q is the current random point vector,
it is placed in the first M elements of Q. }
implementation
end. { of the Rand13_ unit }
Up \SpecFunc