Changeset 1898 for inundation/pyvolution/test_least_squares.py

Timestamp:

Oct 12, 2005, 9:09:20 AM (20 years ago)

Author:

duncan

Message:

comments, test added.

File:

• inundation/pyvolution/test_least_squares.py (modified) (3 diffs)

inundation/pyvolution/test_least_squares.py

@@ -573 +573 @@
         #Point on common vertex
         d6 = [0., 5.]
-
-
+        
         point_coords = [d0, d1, d2, d3, d4, d5, d6]
 
@@ -593 +592 @@
                    [5., 7.5, 7.5, -5.], # (2.5, 2.5)
                    [5., 6., 9., -2.],   # (4,1)
-                   [5., 10., 5., -10.]] # (0,5)
+                   [5., 10., 5., -10.]]  # (0,5)
 
         #print "***********"
@@ -599 +598 @@
         #print "answer",answer
         #print "***********"
+
+        #Should an error message be returned if points are outside
+        # of the mesh? Not currently.  
+
+        assert allclose(z, answer)
+
+
+    def test_interpolate_point_outside_of_mesh(self):
+        """Test linear interpolation of known values at vertices to
+        new points inside a triangle
+        """
+        a = [0.0, 0.0]
+        b = [0.0, 5.0]
+        c = [5.0, 0.0]
+        d = [5.0, 5.0]
+
+        vertices = [a, b, c, d]
+        triangles = [ [1,0,2], [2,3,1] ]   #bac, cdb
+
+        #Far away point
+        d7 = [-1., -1.]
+        
+        point_coords = [ d7]
+
+        interp = Interpolation(vertices, triangles, point_coords)
+
+        #Known values at vertices
+        #Functions are x+y, x+2y, 2x+y, x-y-5
+        f = [ [0., 0., 0., -5.],        # (0,0)
+              [5., 10., 5., -10.],      # (0,5)
+              [5., 5., 10.0, 0.],       # (5,0)
+              [10., 15., 15., -5.]]     # (5,5)
+
+        z = interp.interpolate(f)
+        answer = [ [0., 0., 0., 0.]] # (-1,-1)
+
+        #print "***********"
+        #print "z",z
+        #print "answer",answer
+        #print "***********"
+
+        #Should an error message be returned if points are outside
+        # of the mesh? Not currently.  
 
         assert allclose(z, answer)