Low WAFOM Niederreiter-Xing Sequence

R implementation of Low WAFOM Niederreiter-Xing Sequence, based on Shinsuke Mori, "Suuchi Sekibun no tameno QMC Ten Shuugou no Sekkei, Tansaku, oyobi sono Yuukousei", Master's Thesis, 2017, and Ryuichi Ohori, "Efficient Quasi Monte Carlo Integration by Adjusting the Derivation-sensitivity Parameter of Walsh Figure of Merit", Master's Thesis, 2015.

Porting to R by Mutsuo Saito. The R version does not return coordinate value zero, but returns value very near to zero, 2^-64.

The development of this code is partially supported by JST CREST.

Reference

Sample

Get available dimension number of Low WAFOM Niederreiter-Xing Sequence.

library(LowWAFOMNX)
lowWAFOMNX.dimMinMax()

Get available F2 dimension number of Low WAFOM Niederreiter-Xing Sequence.

lowWAFOMNX.dimF2MinMax(10)

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

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

Get digital shifted points.

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


Try the LowWAFOMNX package in your browser

Any scripts or data that you put into this service are public.

LowWAFOMNX documentation built on May 2, 2019, 6:55 a.m.