|
SSJ V. 1.2.5. |
||||||||
PREV CLASS NEXT CLASS | FRAMES NO FRAMES | ||||||||
SUMMARY: NESTED | FIELD | CONSTR | METHOD | DETAIL: FIELD | CONSTR | METHOD |
java.lang.Objectumontreal.iro.lecuyer.hups.PointSet
umontreal.iro.lecuyer.hups.DigitalNet
umontreal.iro.lecuyer.hups.DigitalNetBase2
umontreal.iro.lecuyer.hups.DigitalSequenceBase2
umontreal.iro.lecuyer.hups.SobolSequence
public class SobolSequence
This class implements digital nets or digital sequences in base 2 formed by the first n = 2k points of a Sobol' sequence. Values of n up to 231 are allowed, in a maximum of 360 dimensions.
In Sobol's proposal, the generator matrices Cj are upper triangular matrices defined by a set of direction numbers
For each dimension j, the integers mj, c are defined by selecting a primitive polynomial over F2 of degree cj,
This list of primitive polynomials, as well as default choices for the direction numbers, are stored in precomputed tables. The ordered list of primitive polynomials was taken from Florent Chabaud's web site, at http://fchabaud.free.fr/. Each polynomial fj(z) is stored in the form of the integer 2cj + aj, 12cj-1 + ... + aj, cj, whose binary representation gives the polynomial coefficients.
For the set of direction numbers, there are several possibilities, based on different selection criteria. The original values proposed by Sobol', and implemented in the code of Bratley and Fox for j <= 40, were selected in terms of his properties A and A', which are equivalent to s-distribution with one and two bits of accuracy, respectively. The default direction numbers used here have been taken from Lemieux et al. For j <= 40, they are the same as in Bratley and Fox.
Constructor Summary | |
---|---|
SobolSequence(int n,
int dim)
Constructs a Sobol point set with at least n points in dim dimensions. |
|
SobolSequence(int k,
int w,
int dim)
Constructs a new digital net with n = 2k points and w = 31 output digits, in dim dimensions, formed by taking the first n points of the Sobol' sequence. |
Method Summary | |
---|---|
void |
extendSequence(int k)
Increases the number of points to n = 2k from now on. |
String |
toString()
Formats a string that contains the information about the point set. |
Methods inherited from class umontreal.iro.lecuyer.hups.DigitalSequenceBase2 |
---|
iteratorShift, iteratorShiftNoGray, toNet, toNetShiftCj |
Methods inherited from class umontreal.iro.lecuyer.hups.DigitalNet |
---|
eraseOriginalGeneratorMatrices, leftMatrixScrambleRest, resetGeneratorMatrices, unrandomize |
Methods inherited from class umontreal.iro.lecuyer.hups.PointSet |
---|
addRandomShift, addRandomShift, formatPoints, formatPoints, getDimension, getNumPoints, randomize, randomize, randomize, randomize |
Methods inherited from class java.lang.Object |
---|
equals, getClass, hashCode, notify, notifyAll, wait, wait, wait |
Constructor Detail |
---|
public SobolSequence(int k, int w, int dim)
k
- there will be 2^k pointsw
- number of output digitsdim
- dimension of the point setpublic SobolSequence(int n, int dim)
dim
- dimension of the point setn
- minimal number of pointsMethod Detail |
---|
public String toString()
PointSet
toString
in class DigitalNetBase2
public void extendSequence(int k)
DigitalSequenceBase2
extendSequence
in class DigitalSequenceBase2
k
- there will be 2^k points
|
SSJ V. 1.2.5. |
||||||||
PREV CLASS NEXT CLASS | FRAMES NO FRAMES | ||||||||
SUMMARY: NESTED | FIELD | CONSTR | METHOD | DETAIL: FIELD | CONSTR | METHOD |