def test_coarsen(self):
# =====================================================================
# Simple Coarsening
# =====================================================================
# Define a new forest with two quadtrees
forest = Forest([Tree(4), Tree(4)])
# Coarsen, nothing should happen
forest.coarsen()
# Check that forest is as it was
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 2)
# Refine and coarsen again
forest.refine()
forest.coarsen()
# Check that forest is as it was
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 2)
# =====================================================================
# Coarsening Flag
# =====================================================================
# Refine and mark one grandchild
forest.refine()
# Check that forest now has 10 Tree
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 10)
forest.get_child(0).mark(1)
forest.coarsen(coarsening_flag=1)
# Tree 0 should not have children, while Tree 1 should have 4
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 6)
self.assertFalse(forest.get_child(0).has_children())
self.assertTrue(forest.get_child(1).has_children())
# Nothing is marked 1
self.assertFalse(any(
child.is_marked(1) for child in forest.traverse()))
# =====================================================================
# Coarsening with subforests
# =====================================================================
forest = Forest([Tree(2), Tree(2)])
forest.refine()
# Make a subforest 0, 00, 01, 1
forest.get_child(0).get_child(0).mark(2)
forest.root_subtrees(2)
count = 0
for node in forest.traverse(flag=2):
count += 1
self.assertEqual(count, 4)
# Coarsen the subforest
forest.coarsen(subforest_flag=2)
# Subforest now contains [0], [1]
count = 0
for dummy in forest.traverse(flag=2):
count += 1
self.assertEqual(count, 2)
# Forest still contains 4 nodes (0,1,00,01,10,11)
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 6)
# New subforest: 0, 00, 01, 1
forest.get_child(0).get_child(0).mark(2)
forest.root_subtrees(2)
# Count subforest nodes
count = 0
for dummy in forest.traverse(2):
count += 1
self.assertEqual(count, 4)
# Coarsening flag at a node not in the subforest
forest.get_child(1).get_child(1).mark(1)
forest.coarsen(subforest_flag=2, coarsening_flag=1)
# Count subforest nodes
count = 0
for dummy in forest.traverse(2):
count += 1
self.assertEqual(count, 4)
# Apply coarsening flag to a node in the subforest
forest.get_child(0).mark(1)
# Coarsen
forest.coarsen(subforest_flag=2, coarsening_flag=1)
# Now there should be 2 subnodes
count = 0
for dummy in forest.traverse(2):
count += 1
self.assertEqual(count, 2)
# Make sure the coarsening flag is deleted.
self.assertFalse(forest.get_child(0).is_marked(1))
# =====================================================================
# Coarsening with new_label
# =====================================================================
# TODO: TEST HERE.
# Subforest is: 0, 00, 01, 1
forest.get_child(0).get_child(0).mark(2)
forest.root_subtrees(2)
# Mark [0,0] with coarsening flag
forest.get_child(0).mark(1)
# Coarsen subforest and label with new_label
forest.coarsen(subforest_flag=2,
coarsening_flag=1,
new_label=3,
debug=True)
# Check that subforest still has the same nodes
count = 0
for dummy in forest.traverse(2):
count += 1
self.assertEqual(count, 4)
# Check that the new submesh has fewer
count = 0
for dummy in forest.traverse(3):
count += 1
self.assertEqual(count, 2)
#
# Now with no submesh
#
# Check that forest still has 6 nodes
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 6)
# Mark with coarsening flag
forest.get_child(0).mark(1)
# Coarsen
forest.coarsen(coarsening_flag=1, new_label=4)
# Check that forest still has 6 nodes
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 6)
# Check that subforest has 5 nodes
count = 0
for dummy in forest.traverse(4):
count += 1
self.assertEqual(count, 4)
def test_refine(self):
# =====================================================================
# Simple Refineement
# =====================================================================
# Define a new forest with two binary trees
forest = Forest([Tree(2), Tree(2)])
# Refine the forest indiscriminantly
forest.refine()
# Check wether the trees have been split
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 6)
# =====================================================================
# Refinement Label
# =====================================================================
# Mark second tree and refine only by its label
forest.get_child(1).mark(1)
forest.refine(refinement_flag=1)
# Nothing should have happened (because child is not a leaf)
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 6)
forest.get_child(1).get_child(1).mark(1)
forest.refine(refinement_flag=1)
count = 0
for dummy in forest.traverse():
count += 1
self.assertEqual(count, 8)
# =====================================================================
# Refinement of subforest
# =====================================================================
forest = Forest([Tree(2), Tree(2)])
forest.refine()
# Define subforest
forest.get_child(0).get_child(0).mark(1)
forest.root_subtrees(1)
# Check node count
count = 0
for dummy in forest.traverse(1):
count += 1
self.assertEqual(count, 4)
forest.refine(subforest_flag=1)
# Check node count
count = 0
for dummy in forest.traverse(1):
count += 1
self.assertEqual(count, 10)
forest.coarsen(subforest_flag=1)
count = 0
for dummy in forest.traverse(1):
count += 1
self.assertEqual(count, 4)
# Now try with a refinement flag
forest.get_child(1).mark(2)
forest.refine(subforest_flag=1, refinement_flag=2)
# Check node count
count = 0
for dummy in forest.traverse(1):
count += 1
self.assertEqual(count, 6)
# =====================================================================
# Refine with new label
# =====================================================================
forest = Forest([Tree(2), Tree(2)])
forest.refine()
# Define subforest
forest.get_child(0).get_child(0).mark(1)
forest.root_subtrees(1)
# Check node count
count = 0
for dummy in forest.traverse(1):
count += 1
self.assertEqual(count, 4)
forest.refine(subforest_flag=1, new_label=4)
# Node count of new label
count = 0
for dummy in forest.traverse(4):
count += 1
self.assertEqual(count, 10)
# Node count for original submesh
count = 0
for dummy in forest.traverse(1):
count += 1
self.assertEqual(count, 4)
# Refinement marker
forest.get_child(1).mark(3)
# Refine
forest.refine(subforest_flag=1, refinement_flag=3, new_label=5)
# Check node count
count = 0
for dummy in forest.traverse(5):
count += 1
self.assertEqual(count, 6)
def test_get_leaves(self):
#
# 1D
#
node = Tree(2)
forest = Forest([node])
leaves = forest.get_leaves()
# Only a ROOT node, it should be the only LEAF
self.assertEqual(leaves, [node], 'Cell should be its own leaf.')
#
# Split cell and L child - find leaves
#
node.split()
l_child = node.get_child(0)
l_child.split()
leaves = forest.get_leaves()
self.assertEqual(len(leaves), 3, 'Cell should have 3 leaves.')
#
# Depth first order
#
addresses_depth_first = [[0, 0, 0], [0, 0, 1], [0, 1]]
leaves = forest.get_leaves(mode='depth-first')
for i in range(len(leaves)):
leaf = leaves[i]
self.assertEqual(leaf.get_node_address(), addresses_depth_first[i],
'Incorrect order, depth first search.')
#
# Breadth first order
#
addresses_breadth_first = [[0, 1], [0, 0, 0], [0, 0, 1]]
leaves = node.get_leaves(mode='breadth-first')
for i in range(len(leaves)):
leaf = leaves[i]
self.assertEqual(leaf.get_node_address(),
addresses_breadth_first[i],
'Incorrect order, breadth first search.')
node.get_child(0).get_child(0).mark('1')
node.get_child(1).mark('1')
node.make_rooted_subtree('1')
leaves = node.get_leaves(subtree_flag='1')
self.assertEqual(len(leaves),2, \
'There should only be 2 flagged leaves')
#
# 2D
#
node = Tree(4)
forest = Forest([node])
#
# Split cell and SW child - find leaves
#
node.split()
sw_child = node.get_child(0)
sw_child.split()
leaves = node.get_leaves()
self.assertEqual(len(leaves), 7, 'Node should have 7 leaves.')
#
# Nested traversal
#
leaves = node.get_leaves()
self.assertEqual(leaves[0].get_node_address(),[0,1], \
'The first leaf in the nested enumeration should have address [1]')
leaves = node.get_leaves(mode='depth-first')
self.assertEqual(leaves[0].get_node_address(), [0,0,0], \
'First leaf in un-nested enumeration should be [0,0].')
#
# Merge SW child - find leaves
#
sw_child.delete_children()
leaves = node.get_leaves()
self.assertEqual(len(leaves), 4, 'Node should have 4 leaves.')
#
# Marked Leaves
#
node = Tree(4)
node.mark(1)
forest = Forest([node])
self.assertTrue(node in forest.get_leaves(flag=1), \
'Node should be a marked leaf node.')
self.assertTrue(node in forest.get_leaves(), \
'Node should be a marked leaf node.')
node.split()
sw_child = node.get_child(0)
sw_child.split()
sw_child.mark(1)
self.assertEqual(node.get_leaves(subtree_flag=1), \
[sw_child], 'SW child should be only marked leaf')
sw_child.remove()
self.assertEqual(forest.get_leaves(subforest_flag=1), \
[node], 'node should be only marked leaf')
#
# Nested traversal
#
node = Tree(4)
node.split()
forest = Forest([node])
for child in node.get_children():
child.split()
node.get_child(1).mark(1, recursive=True)
node.get_child(3).mark(1)
forest.root_subtrees(1)
leaves = forest.get_leaves(subforest_flag=1)
self.assertEqual(len(leaves), 7, 'This tree has 7 flagged LEAF nodes.')
self.assertEqual(leaves[0], node.get_child(0),
'The first leaf should be the NE child.')
self.assertEqual(leaves[3],
node.get_child(1).get_child(0),
'4th flagged leaf should be SE-NW grandchild.')