test_xpart_map.m

Go to the documentation of this file.

function ok = test_xpart_map
  % test for the XPartMap class
  %
  % This test checks especially the methods
  %   - XPartMap.refine() and
  %   - XPartMap.sub_xpart_map()

oned = XPartMap;

oned.refine(0);
ok = all(oned.intervals == [1 0 1 0 0 0.5 0 0.5 1]);
ok = ok & all(oned.leaf_ids  == [1 2]);
ok = ok & all(oned.leaf_enum == [1 2]);
oned.refine([1,2]);
ok = ok & all(oned.intervals == [1 0 1 3 0 0.5 5 0.5 1 0 0 0.25 0 0.25 0.5 0 0.5 0.75 0 0.75 1]);
ok = ok & all(oned.leaf_ids  == [3 5 4 6]);
ok = ok & all(oned.leaf_enum == [3 4 5 6]);
oned.refine([3,5]);
ok = ok & all(oned.intervals == [1 0 1 3 0 0.5 5 0.5 1 7 0 0.25 0 0.25 0.5 9 0.5 0.75 0 0.75 1 0 0 0.125 0 0.125 0.25 0 0.5 0.625 0 0.625 0.75]);
ok = ok & all(oned.leaf_ids  == [7 9 4 6 8 10]);
ok = ok & all(oned.leaf_enum == [7 8 4 9 10 6]);

[sxm, old_leaf_enum] = sub_xpart_map(oned, [0.125 0.4]');
ok = ok & all(sxm.intervals == [1 0.1250 0.4000 2 0.1250 0.4000 4 0.1250 0.2500 0 0.2500 0.4000 0 0.1250 0.2500]);
ok = ok & all(sxm.leaf_ids == [3 4]);
ok = ok & all(sxm.leaf_enum == [4 3]);
ok = ok & all(sxm.refine_levels == [0 1 2 2 3]);

ok = ok & all(old_leaf_enum == [2 3]);

end