290     double absx, r, y, z, sqrtpow;  in m_tgamma()  local
316     absx = fabs(x);  in m_tgamma()
319     if (absx < 1e-20) {  in m_tgamma()
329     if (absx > 200.0) {  in m_tgamma()
339     y = absx + lanczos_g_minus_half;  in m_tgamma()
341     if (absx > lanczos_g_minus_half) {  in m_tgamma()
346         double q = y - absx;  in m_tgamma()
351         z = q - absx;  in m_tgamma()
355         r = -pi / m_sinpi(absx) / absx * exp(y) / lanczos_sum(absx);  in m_tgamma()
357         if (absx < 140.0) {  in m_tgamma()
358             r /= pow(y, absx - 0.5);  in m_tgamma()
361             sqrtpow = pow(y, absx / 2.0 - 0.25);  in m_tgamma()
367         r = lanczos_sum(absx) / exp(y);  in m_tgamma()
369         if (absx < 140.0) {  in m_tgamma()
370             r *= pow(y, absx - 0.5);  in m_tgamma()
373             sqrtpow = pow(y, absx / 2.0 - 0.25);  in m_tgamma()
392     double absx;  in m_lgamma()  local
413     absx = fabs(x);  in m_lgamma()
415     if (absx < 1e-20)  in m_lgamma()
416         return -log(absx);  in m_lgamma()
422     r = log(lanczos_sum(absx)) - lanczos_g;  in m_lgamma()
423     r += (absx - 0.5) * (log(absx + lanczos_g - 0.5) - 1);  in m_lgamma()
426         r = logpi - log(fabs(m_sinpi(absx))) - log(absx) - r;  in m_lgamma()
553     double absx, cf;  in m_erf()
557     absx = fabs(x);  in m_erf()
558     if (absx < ERF_SERIES_CUTOFF)  in m_erf()
561         cf = m_erfc_contfrac(absx);  in m_erf()
575     double absx, cf;  in m_erfc()
579     absx = fabs(x);  in m_erfc()
580     if (absx < ERF_SERIES_CUTOFF)  in m_erfc()
583         cf = m_erfc_contfrac(absx);  in m_erfc()
635         double absx, absy, c, m, r;  in m_remainder()  local
641         absx = fabs(x);  in m_remainder()
643         m = fmod(absx, absy);  in m_remainder()
701             r = m - 2.0 * fmod(0.5 * (absx - m), absy);  in m_remainder()