pub fn ln_gamma(x: f32) -> f32
Natural logarithm of the Gamma function.
Only works for positive values.