Function sub_estimate 

pub fn sub_estimate<S: Simd>(
    simd: S,
    a: Quad<S::f64s>,
    b: Quad<S::f64s>,
) -> Quad<S::f64s>