Changeset 1653
- Timestamp:
- 03/03/10 12:19:27 (5 months ago)
- Files:
-
- 1 modified
-
trunk/yt/lagos/Clump.py (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/yt/lagos/Clump.py
r1146 r1653 33 33 self.data = data 34 34 self.field = field 35 self.min = self.data[field].min()36 self.max = self.data[field].max()35 self.min_val = self.data[field].min() 36 self.max_val = self.data[field].max() 37 37 self.cached_fields = cached_fields 38 38 … … 107 107 f_ptr.write("\n") 108 108 109 def find_children(self, min , max= None):109 def find_children(self, min_val, max_val = None): 110 110 if self.children is not None: 111 111 print "Wiping out existing children clumps." 112 112 self.children = [] 113 if max is None: max = self.max114 contour_info = identify_contours(self.data, self.field, min , max,113 if max_val is None: max_val = self.max_val 114 contour_info = identify_contours(self.data, self.field, min_val, max_val, 115 115 self.cached_fields) 116 116 for cid in contour_info: … … 143 143 def __reduce__(self): 144 144 return (_reconstruct_clump, 145 (self.parent, self.field, self.min , self.max,145 (self.parent, self.field, self.min_val, self.max_val, 146 146 self.function_value, self.children, self.data, self.clump_info, self.function)) 147 147 … … 154 154 if iterable(parent): parent = parent[1] 155 155 if children is None: children = [] 156 obj.parent, obj.field, obj.min , obj.max, obj.function_value, obj.children, obj.clump_info, obj.function = \156 obj.parent, obj.field, obj.min_val, obj.max_val, obj.function_value, obj.children, obj.clump_info, obj.function = \ 157 157 parent, field, mi, ma, function_value, children, clump_info, function 158 158 # Now we override, because the parent/child relationship seems a bit … … 163 163 return obj 164 164 165 def find_clumps(clump, min , max, d_clump):166 print "Finding clumps: min: %e, max: %e, step: %f" % (min , max, d_clump)167 if min >= max: return168 clump.find_children(min )165 def find_clumps(clump, min_val, max_val, d_clump): 166 print "Finding clumps: min: %e, max: %e, step: %f" % (min_val, max_val, d_clump) 167 if min_val >= max_val: return 168 clump.find_children(min_val) 169 169 170 170 if (len(clump.children) == 1): 171 find_clumps(clump, min *d_clump, max, d_clump)171 find_clumps(clump, min_val*d_clump, max_val, d_clump) 172 172 173 173 elif (len(clump.children) > 0): … … 175 175 print "Investigating %d children." % len(clump.children) 176 176 for child in clump.children: 177 find_clumps(child, min *d_clump, max, d_clump)177 find_clumps(child, min_val*d_clump, max_val, d_clump) 178 178 if ((child.children is not None) and (len(child.children) > 0)): 179 179 these_children.append(child)
