def from_ell_to_index(ell):
    """
    Returns the range of column values assuming  a matrix with columns ordered
    with the m multipole , m ranging from -ell to +ell
    """
    return ell ** 2, ell ** 2 + 2 * ell + 1