Function lgamma_around_2 

fn lgamma_around_2(x: f64, parity: f64, f_res: DoubleDouble) -> DoubleDouble

due to log(x) leading terms cancellation happens around 2, hence we’re using different approximation around LogGamma(2). Coefficients are ill-conditioned here so Minimal Newton form is used