Changeset 2684 for inundation/fit_interpolate/test_interpolate.py

Timestamp:

Apr 10, 2006, 3:15:05 PM (19 years ago)

Author:

duncan

Message:

Moving methods around, comments, a new test (to use when interpolate_function can handle discontinuous meshes)

File:

• inundation/fit_interpolate/test_interpolate.py (modified) (2 diffs)

inundation/fit_interpolate/test_interpolate.py

@@ -798 +798 @@
 
 
+    def fix_test_interpolation_precompute_points(self):
+        # looking at a discrete mesh
+        #
+    
+        #Three timesteps
+        time = [0.0, 60.0]    
+
+        #Setup mesh used to represent fitted function
+        points = [[ 15., -20.],
+                  [ 15.,  10.],
+                  [  0., -20.],
+                  [  0.,  10.],
+                  [  0., -20.],
+                  [ 15.,  10.]]
+        
+        triangles = [[0, 1, 2],
+                     [3, 4, 5]]
+
+        #New datapoints where interpolated values are sought
+        interpolation_points = [[ 1.,  0.], [0.,1.]]
+
+        #One quantity
+        Q = zeros( (2,6), Float )
+
+        #Linear in time and space
+        for i, t in enumerate(time):
+            Q[i, :] = t*linear_function(points)
+        #print "Q", Q
+
+
+        
+        interp = Interpolate(points, triangles)
+        f = array([linear_function(points),2*linear_function(points) ])
+        f = transpose(f)
+        #print "f",f
+        z = interp.interpolate(f, interpolation_points)
+        answer = [linear_function(interpolation_points),
+                  2*linear_function(interpolation_points) ]
+        answer = transpose(answer)
+        print "z",z
+        print "answer",answer
+        assert allclose(z, answer)
+
+
+        #Check interpolation of one quantity using interpolaton points)
+        I = Interpolation_function(time, Q,
+                                   vertex_coordinates = points,
+                                   triangles = triangles, 
+                                   interpolation_points = interpolation_points,
+                                   verbose = False)
+        print "I.precomputed_values", I.precomputed_values
+        
+        self.failUnless( I.precomputed_values['Attribute'][1] == 60.0,
+                        ' failed')
+        
     def test_interpolation_function_outside_point(self):
         # Test spatio-temporal interpolation
@@ -911 +966 @@
 
     suite = unittest.makeSuite(Test_Interpolate,'test')
-    #suite = unittest.makeSuite(Test_Interpolate,'test_points_outside_the_polygon')
+    #suite = unittest.makeSuite(Test_Interpolate,'test_interpolation_precompute_points')
     runner = unittest.TextTestRunner(verbosity=1)
     runner.run(suite)