pub fn two_diff(a: f64, b: f64) -> (f64, f64)
Difference of a and b, plus the exact rounding error.
a
b
a - b == difference + error exactly, for any finite inputs.
a - b == difference + error