Esempio n. 1
0
    print(xg.shape)







    fine_grid = row_stack([xg.flatten(), yg.flatten()])

    true_vals = fun(fine_grid)
    interp_vals = sp.interpolate(fine_grid)
#    interp_new_vals = sp.interpolate_new(fine_grid)
    interp_new_vals = interp_vals_3

    interp_sg_vals = sg.interpolate(fine_grid)


#    plot(fine_grid.flatten(), true_vals.flatten())
#    plot(fine_grid.flatten(), interp_vals.flatten())
#    show()


    X = xg
    Y = yg
    Z = (true_vals[0,:]).reshape(X.shape)
    Za = (interp_new_vals[0,:]).reshape(X.shape)
    Zb = (interp_sg_vals[0,:]).reshape(X.shape)
    Zc = (interp_vals[0,:]).reshape(X.shape)