pub fn inv_digamma(x: f64) -> f64