tests/test-Halton.R

In randtoolbox: Toolbox for Pseudo and Quasi Random Number Generation and Random Generator Tests

library(randtoolbox)

#### outputs ####
print(halton(1, 5))

options(digits=15)
n <- 100
n <- 5
cbind( 1/get.primes(n), as.vector(halton(1, n) ) )

#### n argument ####
try(halton(-1))

#### dim argument ####
try(halton(1, 0))


#### init argument ####
halton(n)
randtoolbox:::.getrandtoolboxEnv(".halton.seed")
halton(n, init=TRUE)
randtoolbox:::.getrandtoolboxEnv(".halton.seed")
halton(n, init=FALSE)
try(halton(5, init="a"))

#### mixed argument ####

try(halton(1, mixed=3))
halton(n, mixed=TRUE)

#### usetime argument ####
#QMC with time machine start
halton(n, usetime=TRUE)
halton(n, usetime=TRUE)

#hybrid QMC with SFMT : test continuing the sequence (bug reported by Hiroyuki Kawakatsu)
setSeed(1234); 
halton(n, init=TRUE, mixed=TRUE, mexp=607)
halton(n, init=FALSE, mixed=TRUE)

setSeed(1234); 
halton(n, init=TRUE, mixed=TRUE, mexp=607)
halton(n, init=FALSE, mixed=TRUE) 

#### method argument ####
n <- 5
d <- 4
halton(n, d, method="C")
randtoolbox:::.getrandtoolboxEnv(".halton.seed")
halton(n, d, method="C", start=5)

#### normal argument ####
halton(n, normal=TRUE)
try(halton(3, normal=1))

#### mexp argument ####
halton(n, mexp=607)
try(halton(3, mexp=3))
try(halton(3, mexp="3"))


#### start argument ####
try(halton(3, start="3"))
halton(n, start=0)
halton(n, start=1)