Changeset 6183


Ignore:
Timestamp:
Jan 17, 2009, 2:45:35 PM (10 years ago)
Author:
ole
Message:

A little more testing of interpolate_polyline

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/fit_interpolate/test_interpolate.py

    r6182 r6183  
    19061906                       
    19071907
     1908        # Other and simpler numbers
     1909        f = num.array([1, 5, 13])       
     1910        vertex_coordinates = num.array([[0., 0., ],
     1911                                        [4., 4.],
     1912                                        [8., 8.]])       
     1913        point_coordinates = num.array([[0.1, 0.1],
     1914                                       [3.5, 3.5],
     1915                                       [4.0, 4.0],
     1916                                       [5.2, 5.2],
     1917                                       [7.0, 7.0],
     1918                                       [8.3, 8.3]])
     1919        gauge_neighbour_id = [1, 2, -1]
     1920                                               
     1921        z = interp.interpolate_polyline(f, vertex_coordinates, gauge_neighbour_id, point_coordinates)
     1922        z_ref = [1.1, 4.5, 5., 7.4, 11., 0.]
     1923        #print z
     1924        assert num.allclose(z, z_ref)               
     1925                       
    19081926#-------------------------------------------------------------
    19091927if __name__ == "__main__":
Note: See TracChangeset for help on using the changeset viewer.