Ignore:
Timestamp:
Apr 6, 2010, 8:23:54 PM (14 years ago)
Author:
hudson
Message:

Ticket 113 is complete, and all tests pass.
A centroid list is built by Interpolation_function as it calculates the interpolation matrix, and this is passed out as extra quantities which are output into the gauge.csv file.

File:
1 edited

Legend:

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

    r7673 r7675  
    104104
    105105        interp = Interpolate(points, vertices)
    106         A, _, _ = interp._build_interpolation_matrix_A(data)
     106        A, _, _, _ = interp._build_interpolation_matrix_A(data)
    107107        assert num.allclose(A.todense(), [[1./3, 1./3, 1./3]])
    108108
     
    285285                      0., 0. , 0., 0., 0., 0.]]
    286286
    287         A,_,_ = interp._build_interpolation_matrix_A(data)
     287        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    288288        assert num.allclose(A.todense(), answer)
    289289       
     
    294294                      0., 0. , 0., 0., 0., 0.]]
    295295       
    296         A,_,_ = interp._build_interpolation_matrix_A(data)       
     296        A,_,_,_ = interp._build_interpolation_matrix_A(data)       
    297297        assert num.allclose(A.todense(), answer)
    298298
     
    305305                      0., 0. , 0., 0., 0., 0.]]
    306306                     
    307         A,_,_ = interp._build_interpolation_matrix_A(data)       
     307        A,_,_,_ = interp._build_interpolation_matrix_A(data)       
    308308        assert num.allclose(A.todense(), answer)
    309309
     
    326326                   [0., 0., 1.]]
    327327                   
    328         A,_,_ = interp._build_interpolation_matrix_A(data)
     328        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    329329        assert num.allclose(A.todense(), answer)
    330330
     
    347347        interp = Interpolate(points, vertices)
    348348
    349         A,_,_ = interp._build_interpolation_matrix_A(data)
     349        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    350350        assert num.allclose(A.todense(), answer)
    351351
     
    368368        interp = Interpolate(points, vertices)
    369369
    370         A,_,_ = interp._build_interpolation_matrix_A(data)
     370        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    371371        assert num.allclose(A.todense(), answer)
    372372
     
    387387        #print "interp.get_A()", interp.get_A()
    388388       
    389         A,_,_ = interp._build_interpolation_matrix_A(data)
     389        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    390390        results = A.todense()
    391391        assert num.allclose(num.sum(results, axis=1), 1.0)
     
    410410        answer = [third, third, third]
    411411               
    412         A,_,_ = interp._build_interpolation_matrix_A(data, output_centroids=True)
     412        A,_,_,_ = interp._build_interpolation_matrix_A(data, output_centroids=True)
    413413        results = A.todense()
    414414        assert num.allclose(results, answer)           
     
    430430        interp = Interpolate(points, vertices)
    431431       
    432         A,_,_ = interp._build_interpolation_matrix_A(data)
     432        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    433433        results = A.todense()
    434434        assert num.allclose(num.sum(results, axis=1), [1,1,1,0])
     
    461461
    462462
    463         A,_,_ = interp._build_interpolation_matrix_A(data)
     463        A,_,_,_ = interp._build_interpolation_matrix_A(data)
    464464        A = A.todense()
    465465        for i in range(A.shape[0]):
Note: See TracChangeset for help on using the changeset viewer.