Changeset 8969 for trunk/anuga_core/source/anuga/fit_interpolate/interpolate2d.py

Timestamp:

Sep 10, 2013, 8:12:49 AM (11 years ago)

Author:

steve

Message:

Got rid of comparison warnings with nans in test_interpolate2d.py

File:

• trunk/anuga_core/source/anuga/fit_interpolate/interpolate2d.py (modified) (2 diffs)

trunk/anuga_core/source/anuga/fit_interpolate/interpolate2d.py

@@ -142 +142 @@
 
     # Identify elements that are outside interpolation domain or NaN
+    oldset = numpy.seterr(invalid='ignore')  # Suppress comparison with nan warning
     outside = (xi < x[0]) + (eta < y[0]) + (xi > x[-1]) + (eta > y[-1])
     outside += numpy.isnan(xi) + numpy.isnan(eta)
-
+    numpy.seterr(**oldset)  # Restore warnings
+    
     # Restrict interpolation points to those that are inside the grid
     inside = -outside  # Invert boolean array to find elements inside
@@ -271 +273 @@
 
     if bounds_error:
-        xi0 = min(xi)
-        xi1 = max(xi)
-        eta0 = min(eta)
-        eta1 = max(eta)
+        xi0 = numpy.nanmin(xi)
+        xi1 = numpy.nanmax(xi)
+        eta0 = numpy.nanmin(eta)
+        eta1 = numpy.nanmax(eta)
 
         msg = ('Interpolation point xi=%f was less than the smallest '