Changeset 2660 for inundation/fit_interpolate/interpolate.py

Timestamp:

Apr 5, 2006, 12:07:43 PM (19 years ago)

Author:

duncan

Message:

bug fix

File:

• inundation/fit_interpolate/interpolate.py (modified) (3 diffs)

inundation/fit_interpolate/interpolate.py

@@ -229 +229 @@
                 start=0
                 # creating a dummy array to concatenate to.
-                z = zeros((0,2)) # 2 for points in 2D space
+                
+                f = ensure_numeric(f, Float)
+                z = zeros((0,f.shape[1]))
                 for end in range(start_blocking_len
                                  ,len(point_coordinates)
@@ -235 +237 @@
                     t = self.interpolate_block(f, point_coordinates[start:end],
                                                verbose=verbose)
+                    #print "t", t
+                    #print "z", z 
                     z = concatenate((z,t))
                     start = end
@@ -269 +273 @@
         #print "self.inside_poly_indices", self.inside_poly_indices
         #print "z", z
-        #print "z.size", z.size
         for i in self.outside_poly_indices: 
             z[i] = INF