Ignore:
Timestamp:
Apr 5, 2006, 12:07:43 PM (19 years ago)
Author:
duncan
Message:

bug fix

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/fit_interpolate/test_interpolate.py

    r2659 r2660  
    524524                      2*linear_function(point_coords) ]
    525525            answer = transpose(answer)
     526            #print "z",z
     527            #print "answer",answer
     528            assert allclose(z, answer)
     529           
     530        f = array([linear_function(vertices),2*linear_function(vertices),
     531                   2*linear_function(vertices) - 100  ])
     532        f = transpose(f)
     533        #print "f",f
     534        for blocking_max in range(len(point_coords)+2):
     535        #if True:
     536         #   blocking_max = 5
     537            z = interp.interpolate(f, point_coords,
     538                                   start_blocking_len=blocking_max)
     539            answer = array([linear_function(point_coords),
     540                      2*linear_function(point_coords) ,
     541                      2*linear_function(point_coords)-100 ])
     542            z = transpose(z)
    526543            #print "z",z
    527544            #print "answer",answer
Note: See TracChangeset for help on using the changeset viewer.