Lines Matching refs:absx
303 double absx, r, y, z, sqrtpow;
329 absx = fabs(x);
332 if (absx < 1e-20) {
342 if (absx > 200.0) {
352 y = absx + lanczos_g_minus_half;
354 if (absx > lanczos_g_minus_half) {
359 double q = y - absx;
364 z = q - absx;
368 r = -pi / m_sinpi(absx) / absx * exp(y) / lanczos_sum(absx);
370 if (absx < 140.0) {
371 r /= pow(y, absx - 0.5);
374 sqrtpow = pow(y, absx / 2.0 - 0.25);
380 r = lanczos_sum(absx) / exp(y);
382 if (absx < 140.0) {
383 r *= pow(y, absx - 0.5);
386 sqrtpow = pow(y, absx / 2.0 - 0.25);
405 double absx;
426 absx = fabs(x);
428 if (absx < 1e-20)
429 return -log(absx);
435 r = log(lanczos_sum(absx)) - lanczos_g;
436 r += (absx - 0.5) * (log(absx + lanczos_g - 0.5) - 1);
439 r = logpi - log(fabs(m_sinpi(absx))) - log(absx) - r;
566 double absx, cf;
570 absx = fabs(x);
571 if (absx < ERF_SERIES_CUTOFF)
574 cf = m_erfc_contfrac(absx);
588 double absx, cf;
592 absx = fabs(x);
593 if (absx < ERF_SERIES_CUTOFF)
596 cf = m_erfc_contfrac(absx);
648 double absx, absy, c, m, r;
654 absx = fabs(x);
656 m = fmod(absx, absy);
687 Here absx is exactly halfway between two multiples of absy,
690 absx = n * absy + m
698 0.5 * (absx - m) = (n/2) * absy
703 fmod(0.5 * (absx - m), absy) = |
709 Note that all steps in fmod(0.5 * (absx - m), absy)
714 r = m - 2.0 * fmod(0.5 * (absx - m), absy);