Ignore:
Timestamp:
Jul 1, 2011, 2:52:36 PM (14 years ago)
Author:
steve
Message:

Added in Sudi's most recent 1d codes

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r7930 r8189  
    1212    from numpy import abs, where
    1313
    14     phi = where((abs(a) < abs(b)) & (a*b > 0.0), a, 0.0)
    15     phi = where((abs(b) < abs(a)) & (a*b > 0.0), b, phi)
     14    phi = where((abs(a) < abs(b)) & (a*b >= 0.0), a, 0.0)
     15    phi = where((abs(b) < abs(a)) & (a*b >= 0.0), b, phi)
    1616
    1717    return phi
    … …  
    2020    from numpy import sign, abs, minimum, where
    2121
    22     return where( (sign(a)*sign(b) > 0.0) & (sign(a)*sign(c)>0.0),
     22    return where( (sign(a)*sign(b) >= 0.0) & (sign(a)*sign(c)>0.0),
    2323        sign(a)*minimum(minimum(abs(a),abs(b)),abs(c)), 0.0 )
    2424
Note: See TracChangeset for help on using the changeset viewer.