267     double absx, r, y, z, sqrtpow;  in m_tgamma()  local
293     absx = fabs(x);  in m_tgamma()
296     if (absx < 1e-20) {  in m_tgamma()
306     if (absx > 200.0) {  in m_tgamma()
316     y = absx + lanczos_g_minus_half;  in m_tgamma()
318     if (absx > lanczos_g_minus_half) {  in m_tgamma()
323         double q = y - absx;  in m_tgamma()
328         z = q - absx;  in m_tgamma()
332         r = -pi / m_sinpi(absx) / absx * exp(y) / lanczos_sum(absx);  in m_tgamma()
334         if (absx < 140.0) {  in m_tgamma()
335             r /= pow(y, absx - 0.5);  in m_tgamma()
338             sqrtpow = pow(y, absx / 2.0 - 0.25);  in m_tgamma()
344         r = lanczos_sum(absx) / exp(y);  in m_tgamma()
346         if (absx < 140.0) {  in m_tgamma()
347             r *= pow(y, absx - 0.5);  in m_tgamma()
350             sqrtpow = pow(y, absx / 2.0 - 0.25);  in m_tgamma()
369     double absx;  in m_lgamma()  local
390     absx = fabs(x);  in m_lgamma()
392     if (absx < 1e-20)  in m_lgamma()
393         return -log(absx);  in m_lgamma()
399     r = log(lanczos_sum(absx)) - lanczos_g;  in m_lgamma()
400     r += (absx - 0.5) * (log(absx + lanczos_g - 0.5) - 1);  in m_lgamma()
403         r = logpi - log(fabs(m_sinpi(absx))) - log(absx) - r;  in m_lgamma()
530     double absx, cf;  in m_erf()
534     absx = fabs(x);  in m_erf()
535     if (absx < ERF_SERIES_CUTOFF)  in m_erf()
538         cf = m_erfc_contfrac(absx);  in m_erf()
552     double absx, cf;  in m_erfc()
556     absx = fabs(x);  in m_erfc()
557     if (absx < ERF_SERIES_CUTOFF)  in m_erfc()
560         cf = m_erfc_contfrac(absx);  in m_erfc()
612         double absx, absy, c, m, r;  in m_remainder()  local
618         absx = fabs(x);  in m_remainder()
620         m = fmod(absx, absy);  in m_remainder()
678             r = m - 2.0 * fmod(0.5 * (absx - m), absy);  in m_remainder()