Changeset 6481


Ignore:
Timestamp:
Mar 10, 2009, 1:50:08 PM (15 years ago)
Author:
rwilson
Message:

Changed 'max(a.flat)' to 'num.max(a)' (ditto min()).

Location:
branches/numpy/anuga
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/numpy/anuga/abstract_2d_finite_volumes/domain.py

    r6304 r6481  
    879879
    880880            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))
    883883            msg += '  Histogram:\n'
    884884
     
    892892                else:
    893893                    # Closed upper interval
    894                     hi = max(self.max_speed.flat)
     894                    hi = num.max(self.max_speed)
    895895                    msg += '    [%f, %f]: %d\n' % (lo, hi, count)
    896896
  • branches/numpy/anuga/abstract_2d_finite_volumes/general_mesh.py

    r6463 r6481  
    128128
    129129        msg = 'Vertex indices reference non-existing coordinate sets'
    130         assert max(self.triangles.flat) < self.nodes.shape[0], msg
     130        assert num.max(self.triangles) < self.nodes.shape[0], msg
    131131
    132132        # FIXME: Maybe move to statistics?
     
    530530        Y = C[:,1:6:2].copy()
    531531
    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)
    536536
    537537        return xmin, xmax, ymin, ymax
  • branches/numpy/anuga/shallow_water/data_manager.py

    r6410 r6481  
    22362236        print '  Extent:'
    22372237        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))
    22392239        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))
    22412241        if timestep is not None:
    22422242            print '    t [s] = %f, len(t) == %d' %(times, 1)
     
    23642364
    23652365    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))
    23682368
    23692369    #Assign NODATA_value to all points outside bounding polygon (from interpolation mesh)
     
    26202620        print '  Extent:'
    26212621        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))
    26232623        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))
    26252625        print '    t [s] in [%f, %f], len(t) == %d' \
    26262626              % (min(times), max(times), len(times))
     
    26732673    if verbose:
    26742674        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))
    26762676
    26772677    # Assign NODATA_value to all points outside bounding polygon
     
    31503150        print '  Extent (lat/lon):'
    31513151        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))
    31533153        print '    lon in [%f, %f], len(lon) == %d' \
    3154               % (min(longitudes.flat), max(longitudes.flat),
     3154              % (num.min(longitudes), num.max(longitudes),
    31553155                 len(longitudes.flat))
    31563156        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()
    31603160        name = 'Amplitudes (ha) [cm]'
    31613161        print '  %s in [%f, %f]' % (name, min(q), max(q))
    31623162
    3163         q = uspeed.flat
     3163        q = uspeed.flatten()
    31643164        name = 'Speeds (ua) [cm/s]'
    31653165        print '  %s in [%f, %f]' % (name, min(q), max(q))
    31663166
    3167         q = vspeed.flat
     3167        q = vspeed.flatten()
    31683168        name = 'Speeds (va) [cm/s]'
    31693169        print '  %s in [%f, %f]' % (name, min(q), max(q))
    31703170
    3171         q = elevations.flat
     3171        q = elevations.flatten()
    31723172        name = 'Elevations (e) [m]'
    31733173        print '  %s in [%f, %f]' % (name, min(q), max(q))
     
    32983298        print '  Extent:'
    32993299        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))
    33013301        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))
    33033303        print '    t [s] in [%f, %f], len(t) == %d' \
    33043304              % (min(times), max(times), len(times))
    33053305        print '  Quantities [SI units]:'
    33063306        for name in ['stage', 'xmomentum', 'ymomentum', 'elevation']:
    3307             q = outfile.variables[name][:].flat
     3307            q = outfile.variables[name][:].flatten()
    33083308            print '    %s in [%f, %f]' % (name, min(q), max(q))
    33093309
     
    34683468    fid.close()
    34693469
    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)]
    34713471
    34723472
     
    59375937            print 'Statistics:'
    59385938            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))
    59405940
    59415941    ##
     
    63196319            print 'Statistics:'
    63206320            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))
    63226322
    63236323    ##
Note: See TracChangeset for help on using the changeset viewer.