Changeset 7855 for anuga_work/development/2010-projects/anuga_1d/base/limiters_python.py

Timestamp:

Jun 17, 2010, 5:34:13 PM (14 years ago)

Author:

steve

Message:

Changing for loop to numpy.where

File:

• anuga_work/development/2010-projects/anuga_1d/base/limiters_python.py (modified) (1 diff)

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

@@ -20 +20 @@
     from numpy import sign, abs, minimum, where
 
-    return where( (sign(a)*sign(b) > 0.0) & (sign(a)*sign(c)>0.0), 
-        sign(a)*minimum(abs(a),abs(b),abs(c)), 0.0 )
+    return where( (sign(a)*sign(b) > 0.0) & (sign(a)*sign(c)>0.0),
+        sign(a)*minimum(minimum(abs(a),abs(b)),abs(c)), 0.0 )
+