def METHOD_NAME(points: np.ndarray) -> bool:
    """
    Test if all points are collinear
    Parameters
    ----------
    points
        The array of points
    Returns
    -------
    True if all points are collinear, else False
    """
    # it is expensive to compare all combinations but in reality we will find
    # nonlinear case fast and finish the loop, the probability when having
    # thousands of points, and they are collinear is low
    for x, y, z in itertools.combinations(points, 3):
        if not _collinear(x, y, z):
            return False
    return True