291     double absx, r, y, z, sqrtpow;  in m_tgamma()  local
317     absx = fabs(x);  in m_tgamma()
320     if (absx < 1e-20) {  in m_tgamma()
330     if (absx > 200.0) {  in m_tgamma()
340     y = absx + lanczos_g_minus_half;  in m_tgamma()
342     if (absx > lanczos_g_minus_half) {  in m_tgamma()
347         double q = y - absx;  in m_tgamma()
352         z = q - absx;  in m_tgamma()
356         r = -pi / m_sinpi(absx) / absx * exp(y) / lanczos_sum(absx);  in m_tgamma()
358         if (absx < 140.0) {  in m_tgamma()
359             r /= pow(y, absx - 0.5);  in m_tgamma()
362             sqrtpow = pow(y, absx / 2.0 - 0.25);  in m_tgamma()
368         r = lanczos_sum(absx) / exp(y);  in m_tgamma()
370         if (absx < 140.0) {  in m_tgamma()
371             r *= pow(y, absx - 0.5);  in m_tgamma()
374             sqrtpow = pow(y, absx / 2.0 - 0.25);  in m_tgamma()
393     double absx;  in m_lgamma()  local
414     absx = fabs(x);  in m_lgamma()
416     if (absx < 1e-20)  in m_lgamma()
417         return -log(absx);  in m_lgamma()
423     r = log(lanczos_sum(absx)) - lanczos_g;  in m_lgamma()
424     r += (absx - 0.5) * (log(absx + lanczos_g - 0.5) - 1);  in m_lgamma()
427         r = logpi - log(fabs(m_sinpi(absx))) - log(absx) - r;  in m_lgamma()
554     double absx, cf;  in m_erf()
558     absx = fabs(x);  in m_erf()
559     if (absx < ERF_SERIES_CUTOFF)  in m_erf()
562         cf = m_erfc_contfrac(absx);  in m_erf()
576     double absx, cf;  in m_erfc()
580     absx = fabs(x);  in m_erfc()
581     if (absx < ERF_SERIES_CUTOFF)  in m_erfc()
584         cf = m_erfc_contfrac(absx);  in m_erfc()
636         double absx, absy, c, m, r;  in m_remainder()  local
642         absx = fabs(x);  in m_remainder()
644         m = fmod(absx, absy);  in m_remainder()
702             r = m - 2.0 * fmod(0.5 * (absx - m), absy);  in m_remainder()