## Sobol Sequence

R implementation of S. Joe and F. Y. Kuo, “Constructing Sobol sequences with better two-dimensional projections”, SIAM J. Sci. Comput. 30, 2635-2654 (2008). The implementation is based on the data file new-joe-kuo-6.21201 http://web.maths.unsw.edu.au/~fkuo/sobol/.

Get avairable dimension number of Sobol Sequence.

    library(SobolSequence)
    sobolSequence.dimMinMax()

    ## [1]     2 21201

Get avairable F2 dimension number of Sobol Sequence.

    sobolSequence.dimF2MinMax(10)

    ## [1] 10 31

And get points. Each row of returned matrix contains an s-dimensional point.

    s <- 4
    m <- 10
    c <- 2^m
    mat <- sobolSequence.points(dimR=s, dimF2=m, count=c)
    mat[1,]

    ## [1] 5.421011e-20 5.421011e-20 5.421011e-20 5.421011e-20

Get digital shifted points.

    s <- 4
    m <- 10
    c <- 2^m
    mat <- sobolSequence.points(dimR=s, dimF2=m, count=c, digitalShift=TRUE)
    mat[1,]

    ## [1] 0.71985296 0.57273998 0.55182147 0.06826342