Lines Matching refs:absx
233 double absx, r, y, z, sqrtpow;
258 absx = fabs(x);
261 if (absx < 1e-20) {
271 if (absx > 200.0) {
281 y = absx + lanczos_g_minus_half;
283 if (absx > lanczos_g_minus_half) {
288 double q = y - absx;
293 z = q - absx;
297 r = -pi / sinpi(absx) / absx * exp(y) / lanczos_sum(absx);
299 if (absx < 140.0) {
300 r /= pow(y, absx - 0.5);
303 sqrtpow = pow(y, absx / 2.0 - 0.25);
309 r = lanczos_sum(absx) / exp(y);
311 if (absx < 140.0) {
312 r *= pow(y, absx - 0.5);
315 sqrtpow = pow(y, absx / 2.0 - 0.25);
333 double r, absx;
354 absx = fabs(x);
356 if (absx < 1e-20)
357 return -log(absx);
369 r = log(pi) - log(fabs(sinpi(absx))) - log(absx) -
370 (log(lanczos_sum(absx)) - lanczos_g +
371 (absx-0.5)*(log(absx+lanczos_g-0.5)-1));
493 double absx, cf;
497 absx = fabs(x);
498 if (absx < ERF_SERIES_CUTOFF)
501 cf = m_erfc_contfrac(absx);
511 double absx, cf;
515 absx = fabs(x);
516 if (absx < ERF_SERIES_CUTOFF)
519 cf = m_erfc_contfrac(absx);
