Function ln_gamma_in 

fn ln_gamma_in<R: Scalar>(x: R) -> R

ln Γ(x) for x ≥ ½ by Lanczos’ approximation (g = 7, nine terms), written once against the mode: passively accurate to ~1e-15, and its derivatives are the approximant’s, digamma and trigamma to ~1e-13.