lowdisc_plotelembox

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

Examples

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

Bibliography