pub fn is_close(a: f64, b: f64, atol: f64, rtol: f64) -> bool
Check whether two values are within tolerance: |a - b| <= atol + rtol * |b|.
|a - b| <= atol + rtol * |b|