fn lgamma_0p5_to_1( dx: f64, v_log: DoubleDouble, f_res: DoubleDouble, sum_parity: f64, ) -> DoubleDouble