Function residual_is_zero 

pub fn residual_is_zero(
    a: &[Vec<Rational>],
    x: &[Rational],
    b: &[Rational],
) -> bool

Is A·x == b exactly, component by component? The independent re-check behind a certified solve.