def METHOD_NAME(self, point: np.ndarray[Any, np.dtype[np.float64]]) -> float:
    """Computes the winding number of a closed triangulated surface at given point.
    Parameters:
        point: ``shape=(1, 3)``
            The point being tested.
    Returns:
        The winding number generalized to R^3. Its absolute value ``|wn|`` is 0 for
        points outside the polyhedron, 1 for points inside non-convex polyhedron and
        ``|wn| > 1`` for points inside overlapping polyhedron.
    """
    R = self.vertices - point
    R_triangles = [R[indexes] for indexes in self.connectivity]
    angles = np.array(list(map(self.solid_angle, R_triangles)))
    wn = np.sum(angles) / (4.0 * np.pi)
    return wn