def find_neighbors(pindex, tri):
    """
    Adapted from @user2535797 at https://stackoverflow.com/questions/12374781/how-to-find-all-neighbors-of-a-given-point-\
    in-a-delaunay-triangulation-using-sci

    Parameters
    ----------
        pindex: integer
            a index of a single point
        tri:
            output from scipy.spatial.Delaunay() method

    Returns
    -------
        neighbors: numpy array
            neighbors for a data point

    """
    return tri.vertex_neighbor_vertices[1][tri.vertex_neighbor_vertices[0][pindex]:tri.vertex_neighbor_vertices[0][pindex+1]]