Changeset 4882
- Timestamp:
- Dec 11, 2007, 5:01:43 PM (17 years ago)
- Location:
- anuga_core/source/anuga
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/fit_interpolate/fit.py
r4861 r4882 262 262 AtA = self.AtA #Did this for speed, did ~nothing 263 263 self.point_count += point_coordinates.shape[0] 264 #print "_build_matrix_AtA_Atz - self.point_count", self.point_count 265 if verbose: print 'Getting indices inside mesh boundary' 266 #print 'point_coordinates.shape', point_coordinates.shape 267 #print 'self.mesh.get_boundary_polygon()',\ 268 # self.mesh.get_boundary_polygon() 264 265 #if verbose: print 'Getting indices inside mesh boundary' 269 266 270 267 inside_poly_indices, outside_poly_indices = \ … … 281 278 triangles = self.mesh.triangles #Did this for speed, did ~nothing 282 279 for d, i in enumerate(inside_poly_indices): 283 # For each data_coordinate point284 if verbose and d%((n+10)/10)==0: print 'Doing %d of %d' %(d, n)280 # For each data_coordinate point 281 # if verbose and d%((n+10)/10)==0: print 'Doing %d of %d' %(d, n) 285 282 x = point_coordinates[i] 286 283 element_found, sigma0, sigma1, sigma2, k = \ -
anuga_core/source/anuga/utilities/polygon.py
r4856 r4882 73 73 """ 74 74 75 if verbose: print 'Checking input to inside_polygon'75 #if verbose: print 'Checking input to inside_polygon' 76 76 77 77 try: … … 137 137 """ 138 138 139 if verbose: print 'Checking input to outside_polygon'139 #if verbose: print 'Checking input to outside_polygon' 140 140 try: 141 141 points = ensure_numeric(points, Float) … … 179 179 """ 180 180 181 if verbose: print 'Checking input to outside_polygon'181 #if verbose: print 'Checking input to outside_polygon' 182 182 try: 183 183 points = ensure_numeric(points, Float) … … 258 258 259 259 260 if verbose: print 'Checking input to separate_points_by_polygon'260 #if verbose: print 'Checking input to separate_points_by_polygon' 261 261 262 262 … … 302 302 303 303 304 305 if verbose: print 'Allocating array for indices'306 307 304 indices = zeros( M, Int ) 308 309 #if verbose: print 'Calling C-version of inside poly'310 305 311 306 from polygon_ext import separate_points_by_polygon as sep_points … … 430 425 """ 431 426 432 if verbose: print 'Checking input to poly_xy'427 #if verbose: print 'Checking input to poly_xy' 433 428 434 429 try:
Note: See TracChangeset
for help on using the changeset viewer.