Plot elementary box
lowdisc_plotelembox(b,d) lowdisc_plotelembox(b,d,u)
a matrix of doubles, integer value, the basis
a matrix of doubles, integer value, the log-b of number of divisions in direction 1 and 2
a npoints-by-2 matrix of doubles, the point set to plot. Default is to plot no point set.
Plot the elementary interval with volume 1/prod(b.^d). This interval has b(1)^d(1) divisions in direction 1 and b(2)^d(2) divisions in direction 2.
The values of b and d are expanded as need. For example, b can be a 1-by-1 matrix and d can be a 2-by-1 matrix. In this case, b is expanded to match the size of d.
If the point set u is provided, we plot it.
// Plot an elementary interval with volume 1/2^3 lowdisc_plotelembox(2,[2 1]) // Use a different basis for each direction. // This is useful for Halton sequence. scf(); lowdisc_plotelembox([2 3],[2 1]) // Plot all elementary intervals with volume 1/b^m=1/2^3 b = 2; m = 3; C = [ 0. 3. 1. 2. 2. 1. 3. 0. ]; n = size(C,"r"); for i = 1 : n scf(); lowdisc_plotelembox(b,C(i,:)); end // Plot the Halton sequence in 2 dimensions scf(); u=lowdisc_ldgen(2*3^2,2,"halton"); lowdisc_plotelembox([2 3],[1 2],u) | ![]() | ![]() |
"Random number generation and quasi-Monte Carlo methods", H. Niederreiter, CBMS-NSF Series in Applied Mathematics, No. 63, SIAM, Philadelphia, 1992.