Changeset 7720
- Timestamp:
- May 12, 2010, 10:13:37 AM (13 years ago)
- Location:
- anuga_core/source/anuga
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/fit_interpolate/test_search_functions.py
r7718 r7720 152 152 x = ensure_numeric([0.5, 0.5]) 153 153 154 triangles = root._trilist_from_data(root.search(x))155 156 154 found, sigma0, sigma1, sigma2, k = \ 157 root._search_triangles_of_vertices( triangles, x)155 root._search_triangles_of_vertices(root.search(x), x) 158 156 159 157 if k >= 0: -
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 -
anuga_core/source/anuga/geometry/test_geometry.py
r7716 r7720 86 86 cell = Cell(AABB(0,10, 0,5), None) 87 87 cell.insert([(AABB(1,3, 1, 3), 111), (AABB(8,9, 1, 2), 222)]) 88 88 89 assert len(cell.children) == 2 89 90 assert cell.parent == None -
anuga_core/source/anuga/geometry/test_meshquad.py
r7718 r7720 46 46 assert type(result) in [types.ListType,types.TupleType],\ 47 47 'should be a list' 48 index, parent = result[0]49 self.assertEqual( index, 1)48 49 self.assertEqual(result[0][0][0], 1) 50 50 51 51 … … 138 138 results = Q.search([4.5, 3]) 139 139 assert len(results) == 1 140 idx, _ = results[0] 141 self.assertEqual(idx, 2) 140 self.assertEqual(results[0][0][0], 2) 142 141 results = Q.search([5,4.]) 143 142 self.assertEqual(len(results),1) 144 idx, _ = results[0] 145 self.assertEqual(idx, 2) 143 self.assertEqual(results[0][0][0], 2) 146 144 147 145 def NOtest_num_visits(self):
Note: See TracChangeset
for help on using the changeset viewer.