Ignore:
Timestamp:
Jun 17, 2010, 5:34:13 PM (13 years ago)
Author:
steve
Message:

Changing for loop to numpy.where

File:
1 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/development/2010-projects/anuga_1d/base/limiters_python.py

    r7852 r7855  
    2020    from numpy import sign, abs, minimum, where
    2121
    22     return where( (sign(a)*sign(b) > 0.0) & (sign(a)*sign(c)>0.0),
    23         sign(a)*minimum(abs(a),abs(b),abs(c)), 0.0 )
     22    return where( (sign(a)*sign(b) > 0.0) & (sign(a)*sign(c)>0.0),
     23        sign(a)*minimum(minimum(abs(a),abs(b)),abs(c)), 0.0 )
     24
    2425
    2526
Note: See TracChangeset for help on using the changeset viewer.