<< lowdisc_plotbmbox Support Functions lowdisc_proj2d >>

Low Discrepancy >> Low Discrepancy > Support Functions > lowdisc_plotelembox

lowdisc_plotelembox

Plot elementary box

Calling Sequence

lowdisc_plotelembox(b,d)
lowdisc_plotelembox(b,d,u)

Parameters

b :

a matrix of doubles, integer value, the basis

d:

a matrix of doubles, integer value, the log-b of number of divisions in direction 1 and 2

u :

a npoints-by-2 matrix of doubles, the point set to plot. Default is to plot no point set.

Description

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.

Examples

// 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

Bibliography

"Random number generation and quasi-Monte Carlo methods", H. Niederreiter, CBMS-NSF Series in Applied Mathematics, No. 63, SIAM, Philadelphia, 1992.

Authors


Report an issue
<< lowdisc_plotbmbox Support Functions lowdisc_proj2d >>