def METHOD_NAME():
    r"""Test find natural neighbors for a point function."""
    x = list(range(10))
    y = list(range(10))
    gx, gy = np.meshgrid(x, y)
    pts = np.vstack([gx.ravel(), gy.ravel()]).T
    tri = Delaunay(pts)
    tri_match = tri.find_simplex([4.5, 4.5])
    nn = find_nn_triangles_point(tri, tri_match, [4.5, 4.5])
    # Can't rely on simplex indices, so need to check point indices
    truth = {(45, 55, 44), (55, 54, 44)}
    assert {tuple(verts) for verts in tri.simplices[nn]} == truth