Changeset 7720 for anuga_core/source/anuga/geometry/quad.py
- Timestamp:
- May 12, 2010, 10:13:37 AM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/geometry/quad.py
r7719 r7720 84 84 85 85 # option 1 - try splitting 4 ways 86 subregion11, subregion12 = subregion1.split()87 subregion21, subregion22 = subregion2.split()88 regions = [subregion11, subregion12, subregion21, subregion22]89 for region in regions:90 if region.is_trivial_in(new_region):91 self.children = [Cell(x, parent=self) for x in regions]92 self._insert(new_leaf)93 return86 #subregion11, subregion12 = subregion1.split() 87 #subregion21, subregion22 = subregion2.split() 88 #regions = [subregion11, subregion12, subregion21, subregion22] 89 #for region in regions: 90 #if region.is_trivial_in(new_region): 91 #self.children = [Cell(x, parent=self) for x in regions] 92 #self._insert(new_leaf) 93 #return 94 94 95 95 # option 2 - try splitting 2 ways 96 #if subregion1.is_trivial_in(new_region):97 #self.children = [Cell(subregion1, self), Cell(subregion2, self)]98 #self.children[0]._insert(new_leaf)99 #return100 #elif subregion2.is_trivial_in(new_region):101 #self.children = [Cell(subregion1, self), Cell(subregion2, self)]102 #self.children[1]._insert(new_leaf)103 #return96 if subregion1.is_trivial_in(new_region): 97 self.children = [Cell(subregion1, self), Cell(subregion2, self)] 98 self.children[0]._insert(new_leaf) 99 return 100 elif subregion2.is_trivial_in(new_region): 101 self.children = [Cell(subregion1, self), Cell(subregion2, self)] 102 self.children[1]._insert(new_leaf) 103 return 104 104 105 105 # recursion ended without finding a fit, so attach it as a leaf
Note: See TracChangeset
for help on using the changeset viewer.