Changeset 3941

Timestamp:

Nov 8, 2006, 4:21:45 PM (17 years ago)

Author:

duncan

Message:

change Interpolation_function (in interpolate.py) so it accepts geospatial data objects. See ticket#108

Location:

anuga_core/source/anuga/fit_interpolate

Files:

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

anuga_core/source/anuga/fit_interpolate/interpolate.py

@@ -459 +459 @@
             self.spatial = False
         else:    
-            vertex_coordinates = ensure_numeric(vertex_coordinates)
+            vertex_coordinates = ensure_absolute(vertex_coordinates)
 
             assert triangles is not None, 'Triangles array must be specified'

anuga_core/source/anuga/fit_interpolate/test_interpolate.py

@@ -1012 +1012 @@
         #Three timesteps
         time = [1.0, 5.0, 6.0]
-        
-        
+               
         #Setup mesh used to represent fitted function
         a = [0.0, 0.0]
@@ -1054 +1053 @@
             for id in range(len(interpolation_points)):
                 assert allclose(I(t, id), answer[id])
-
             t += 0.1    
-
 
         try:    
@@ -1064 +1061 @@
         else:
             raise 'Should raise exception'
+
             
-
-
     def test_interpolation_interface(self):
         # Test spatio-temporal interpolation