Function lgamma_1_to_4 

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