Function is_contraction 

fn is_contraction(operator: &DMatrix<Complex>) -> bool

Whether ||L U||_inf < 1 holds for the operator I - L U, which proves it invertible.

Each row sum of |re| + |im| over the entries of L U rounds upward, so the sums bound the infinity norm from above: rounding can make the test fail, never pass.