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)