Changeset 3018 for inundation/fit_interpolate

Timestamp:

May 30, 2006, 2:48:18 PM (19 years ago)

Author:

duncan

Message:

interpolate: Warning message if points are outside the mesh

Location:

inundation/fit_interpolate

Files:

• interpolate.py (modified) (1 diff)
• test_interpolate.py (modified) (3 diffs)

inundation/fit_interpolate/interpolate.py

@@ -235 +235 @@
         #print "self.outside_poly_indices",self.outside_poly_indices 
         #Build n x m interpolation matrix
+        if verbose and len(self.outside_poly_indices) >0:
+            print '\n WARNING: Points outside mesh boundary. \n'
         m = self.mesh.coordinates.shape[0] #Nbr of basis functions (1/vertex)
         n = point_coordinates.shape[0]     #Nbr of data points

inundation/fit_interpolate/test_interpolate.py

@@ -486 +486 @@
         #print "***********"
 
-        #Should an error message be returned if points are outside
-        # of the mesh? Not currently.  
-
         assert allclose(z, answer)
 
@@ -523 +520 @@
               [10., 15., 15., -5.]]     # (5,5)
 
-        z = interp.interpolate(f, point_coords)
+        z = interp.interpolate(f, point_coords) #, verbose=True)
         answer = array([ [INF, INF, INF, INF]]) # (-1,-1)
 
@@ -532 +529 @@
 
         #Should an error message be returned if points are outside
-        # of the mesh? Not currently.  
+        # of the mesh?
+        # A warning message is printed, if verbose is on.
 
         for i in range(4):