Lines Matching refs:absx
233 double absx, r, y, z, sqrtpow; in m_tgamma() local
258 absx = fabs(x); in m_tgamma()
261 if (absx < 1e-20) { in m_tgamma()
271 if (absx > 200.0) { in m_tgamma()
281 y = absx + lanczos_g_minus_half; in m_tgamma()
283 if (absx > lanczos_g_minus_half) { in m_tgamma()
288 double q = y - absx; in m_tgamma()
293 z = q - absx; in m_tgamma()
297 r = -pi / sinpi(absx) / absx * exp(y) / lanczos_sum(absx); in m_tgamma()
299 if (absx < 140.0) { in m_tgamma()
300 r /= pow(y, absx - 0.5); in m_tgamma()
303 sqrtpow = pow(y, absx / 2.0 - 0.25); in m_tgamma()
309 r = lanczos_sum(absx) / exp(y); in m_tgamma()
311 if (absx < 140.0) { in m_tgamma()
312 r *= pow(y, absx - 0.5); in m_tgamma()
315 sqrtpow = pow(y, absx / 2.0 - 0.25); in m_tgamma()
333 double r, absx; in m_lgamma() local
354 absx = fabs(x); in m_lgamma()
356 if (absx < 1e-20) in m_lgamma()
357 return -log(absx); in m_lgamma()
369 r = log(pi) - log(fabs(sinpi(absx))) - log(absx) - in m_lgamma()
370 (log(lanczos_sum(absx)) - lanczos_g + in m_lgamma()
371 (absx-0.5)*(log(absx+lanczos_g-0.5)-1)); in m_lgamma()
492 double absx, cf; in m_erf() local
496 absx = fabs(x); in m_erf()
497 if (absx < ERF_SERIES_CUTOFF) in m_erf()
500 cf = m_erfc_contfrac(absx); in m_erf()
510 double absx, cf; in m_erfc() local
514 absx = fabs(x); in m_erfc()
515 if (absx < ERF_SERIES_CUTOFF) in m_erfc()
518 cf = m_erfc_contfrac(absx); in m_erfc()