def METHOD_NAME(test_data):
    r"""Test find natural neighbors for a point interpolation function."""
    xp, yp, z = test_data
    tri = Delaunay(list(zip(xp, yp)))
    sim_gridx = [30]
    sim_gridy = [30]
    members, tri_info = find_natural_neighbors(tri,
                                               list(zip(sim_gridx, sim_gridy)))
    val = natural_neighbor_point(xp, yp, z, (sim_gridx[0], sim_gridy[0]),
                                 tri, members[0], tri_info)
    truth = 1.009
    assert_almost_equal(truth, val, 3)