PI                 74 missing/lgamma_r.c         s = sin(PI * f);
PI                 88 missing/tgamma.c         return sign * PI / (sin(PI * f) * exp(loggamma(1 - x)));