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 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 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 b^m=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 | ![]() | ![]() |
"Random number generation and quasi-Monte Carlo methods", H. Niederreiter, CBMS-NSF Series in Applied Mathematics, No. 63, SIAM, Philadelphia, 1992.