pub fn lu_is_exact(a: &[Vec<Rational>], f: &Lu) -> bool
Certify a factorization: is P·A == L·U exactly? (P·A is A with its rows permuted by perm.)
P·A == L·U
perm