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)