Lines Matching refs:absx
258 double absx, r, y, z, sqrtpow;
284 absx = fabs(x);
287 if (absx < 1e-20) {
297 if (absx > 200.0) {
307 y = absx + lanczos_g_minus_half;
309 if (absx > lanczos_g_minus_half) {
314 double q = y - absx;
319 z = q - absx;
323 r = -pi / sinpi(absx) / absx * exp(y) / lanczos_sum(absx);
325 if (absx < 140.0) {
326 r /= pow(y, absx - 0.5);
329 sqrtpow = pow(y, absx / 2.0 - 0.25);
335 r = lanczos_sum(absx) / exp(y);
337 if (absx < 140.0) {
338 r *= pow(y, absx - 0.5);
341 sqrtpow = pow(y, absx / 2.0 - 0.25);
359 double r, absx;
380 absx = fabs(x);
382 if (absx < 1e-20)
383 return -log(absx);
389 r = log(lanczos_sum(absx)) - lanczos_g;
390 r += (absx - 0.5) * (log(absx + lanczos_g - 0.5) - 1);
393 r = logpi - log(fabs(sinpi(absx))) - log(absx) - r;
513 double absx, cf;
517 absx = fabs(x);
518 if (absx < ERF_SERIES_CUTOFF)
521 cf = m_erfc_contfrac(absx);
531 double absx, cf;
535 absx = fabs(x);
536 if (absx < ERF_SERIES_CUTOFF)
539 cf = m_erfc_contfrac(absx);
