X-Git-Url: http://nsz.repo.hu/git/?p=libm;a=blobdiff_plain;f=src%2Fmath%2F__rem_pio2l.h;h=37f3bd28c19610c3ac9c272aa30c5e85dcafca1e;hp=3a747b18eda5e1d05f2c8503aa63001d2cddaa38;hb=59c2bd747bce9b1e35c893b2ace9bea57523a901;hpb=fc04675a09587031d5daebcb877b73976890b71d diff --git a/src/math/__rem_pio2l.h b/src/math/__rem_pio2l.h index 3a747b1..37f3bd2 100644 --- a/src/math/__rem_pio2l.h +++ b/src/math/__rem_pio2l.h @@ -13,11 +13,11 @@ * Optimized by Bruce D. Evans. */ #include "libm.h" -#if LD80 +#if LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384 /* ld80 version of __rem_pio2(x,y) * * return the remainder of x rem pi/2 in y[0]+y[1] - * use __rem_pio2_slow() + * use __rem_pio2_large() for large x */ #define BIAS (LDBL_MAX_EXP - 1) @@ -135,7 +135,7 @@ static inline int __rem_pio2l(long double x, long double *y) nx = 3; while (tx[nx-1] == zero) nx--; /* skip zero term */ - n = __rem_pio2_slow(tx,ty,e0,nx,2); + n = __rem_pio2_large(tx,ty,e0,nx,2); r = (long double)ty[0] + ty[1]; w = ty[1] - (r - ty[0]); if (expsign < 0) {