Changeset 6481
- Timestamp:
- Mar 10, 2009, 1:50:08 PM (15 years ago)
- Location:
- branches/numpy/anuga
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/numpy/anuga/abstract_2d_finite_volumes/domain.py
r6304 r6481 879 879 880 880 msg += '------------------------------------------------\n' 881 msg += ' Speeds in [%f, %f]\n' % ( min(self.max_speed.flat),882 max(self.max_speed.flat))881 msg += ' Speeds in [%f, %f]\n' % (num.min(self.max_speed), 882 num.max(self.max_speed)) 883 883 msg += ' Histogram:\n' 884 884 … … 892 892 else: 893 893 # Closed upper interval 894 hi = max(self.max_speed.flat)894 hi = num.max(self.max_speed) 895 895 msg += ' [%f, %f]: %d\n' % (lo, hi, count) 896 896 -
branches/numpy/anuga/abstract_2d_finite_volumes/general_mesh.py
r6463 r6481 128 128 129 129 msg = 'Vertex indices reference non-existing coordinate sets' 130 assert max(self.triangles.flat) < self.nodes.shape[0], msg130 assert num.max(self.triangles) < self.nodes.shape[0], msg 131 131 132 132 # FIXME: Maybe move to statistics? … … 530 530 Y = C[:,1:6:2].copy() 531 531 532 xmin = min(X.flat)533 xmax = max(X.flat)534 ymin = min(Y.flat)535 ymax = max(Y.flat)532 xmin = num.min(X) 533 xmax = num.max(X) 534 ymin = num.min(Y) 535 ymax = num.max(Y) 536 536 537 537 return xmin, xmax, ymin, ymax -
branches/numpy/anuga/shallow_water/data_manager.py
r6410 r6481 2236 2236 print ' Extent:' 2237 2237 print ' x [m] in [%f, %f], len(x) == %d'\ 2238 %( min(x.flat), max(x.flat), len(x.flat))2238 %(num.min(x), num.max(x), len(x.flat)) 2239 2239 print ' y [m] in [%f, %f], len(y) == %d'\ 2240 %( min(y.flat), max(y.flat), len(y.flat))2240 %(num.min(y), num.max(y), len(y.flat)) 2241 2241 if timestep is not None: 2242 2242 print ' t [s] = %f, len(t) == %d' %(times, 1) … … 2364 2364 2365 2365 if verbose: 2366 print 'Interpolated values are in [%f, %f]' %( min(grid_values.flat),2367 max(grid_values.flat))2366 print 'Interpolated values are in [%f, %f]' %(num.min(grid_values), 2367 num.max(grid_values)) 2368 2368 2369 2369 #Assign NODATA_value to all points outside bounding polygon (from interpolation mesh) … … 2620 2620 print ' Extent:' 2621 2621 print ' x [m] in [%f, %f], len(x) == %d' \ 2622 % ( min(x.flat), max(x.flat), len(x.flat))2622 % (num.min(x), num.max(x), len(x.flat)) 2623 2623 print ' y [m] in [%f, %f], len(y) == %d' \ 2624 % ( min(y.flat), max(y.flat), len(y.flat))2624 % (num.min(y), num.max(y), len(y.flat)) 2625 2625 print ' t [s] in [%f, %f], len(t) == %d' \ 2626 2626 % (min(times), max(times), len(times)) … … 2673 2673 if verbose: 2674 2674 print 'Interpolated values are in [%f, %f]' \ 2675 % ( min(interpolated_values.flat), max(interpolated_values.flat))2675 % (num.min(interpolated_values), num.max(interpolated_values)) 2676 2676 2677 2677 # Assign NODATA_value to all points outside bounding polygon … … 3150 3150 print ' Extent (lat/lon):' 3151 3151 print ' lat in [%f, %f], len(lat) == %d' \ 3152 % ( min(latitudes.flat), max(latitudes.flat), len(latitudes.flat))3152 % (num.min(latitudes), num.max(latitudes), len(latitudes.flat)) 3153 3153 print ' lon in [%f, %f], len(lon) == %d' \ 3154 % ( min(longitudes.flat), max(longitudes.flat),3154 % (num.min(longitudes), num.max(longitudes), 3155 3155 len(longitudes.flat)) 3156 3156 print ' t in [%f, %f], len(t) == %d' \ 3157 % ( min(times.flat), max(times.flat), len(times.flat))3158 3159 q = amplitudes.flat 3157 % (num.min(times), num.max(times), len(times.flat)) 3158 3159 q = amplitudes.flatten() 3160 3160 name = 'Amplitudes (ha) [cm]' 3161 3161 print ' %s in [%f, %f]' % (name, min(q), max(q)) 3162 3162 3163 q = uspeed.flat 3163 q = uspeed.flatten() 3164 3164 name = 'Speeds (ua) [cm/s]' 3165 3165 print ' %s in [%f, %f]' % (name, min(q), max(q)) 3166 3166 3167 q = vspeed.flat 3167 q = vspeed.flatten() 3168 3168 name = 'Speeds (va) [cm/s]' 3169 3169 print ' %s in [%f, %f]' % (name, min(q), max(q)) 3170 3170 3171 q = elevations.flat 3171 q = elevations.flatten() 3172 3172 name = 'Elevations (e) [m]' 3173 3173 print ' %s in [%f, %f]' % (name, min(q), max(q)) … … 3298 3298 print ' Extent:' 3299 3299 print ' x [m] in [%f, %f], len(x) == %d' \ 3300 % ( min(x.flat), max(x.flat), len(x.flat))3300 % (num.min(x), num.max(x), len(x.flat)) 3301 3301 print ' y [m] in [%f, %f], len(y) == %d' \ 3302 % ( min(y.flat), max(y.flat), len(y.flat))3302 % (num.min(y), num.max(y), len(y.flat)) 3303 3303 print ' t [s] in [%f, %f], len(t) == %d' \ 3304 3304 % (min(times), max(times), len(times)) 3305 3305 print ' Quantities [SI units]:' 3306 3306 for name in ['stage', 'xmomentum', 'ymomentum', 'elevation']: 3307 q = outfile.variables[name][:].flat 3307 q = outfile.variables[name][:].flatten() 3308 3308 print ' %s in [%f, %f]' % (name, min(q), max(q)) 3309 3309 … … 3468 3468 fid.close() 3469 3469 3470 return [min(x), max(x), min(y), max(y), min(stage.flat), max(stage.flat)]3470 return [min(x), max(x), min(y), max(y), num.min(stage), num.max(stage)] 3471 3471 3472 3472 … … 5937 5937 print 'Statistics:' 5938 5938 print ' t in [%f, %f], len(t) == %d' \ 5939 % ( min(times.flat), max(times.flat), len(times.flat))5939 % (num.min(times), num.max(times), len(times.flat)) 5940 5940 5941 5941 ## … … 6319 6319 print 'Statistics:' 6320 6320 print ' t in [%f, %f], len(t) == %d' \ 6321 % ( min(times.flat), max(times.flat), len(times.flat))6321 % (num.min(times), num.max(times), len(times.flat)) 6322 6322 6323 6323 ##
Note: See TracChangeset
for help on using the changeset viewer.