pub fn psi_1piy(x: f64) -> f64
This routine computes the real part of the digamma function on the line 1+i y, \Re[\psi(1 + i y)].