From 0498709e07bf9dd0546a61220cd2a27fcd186ce0 Mon Sep 17 00:00:00 2001 From: nsz Date: Mon, 5 Mar 2012 04:11:38 +0100 Subject: [PATCH] rename __rem_pio2_slow to __rem_pio2_large --- src/math/__rem_pio2.c | 4 ++-- src/math/{__rem_pio2_slow.c => __rem_pio2_large.c} | 6 +++--- src/math/__rem_pio2f.c | 4 ++-- src/math/__rem_pio2l.h | 4 ++-- 4 files changed, 9 insertions(+), 9 deletions(-) rename src/math/{__rem_pio2_slow.c => __rem_pio2_large.c} (99%) diff --git a/src/math/__rem_pio2.c b/src/math/__rem_pio2.c index 4373e22..a7d779e 100644 --- a/src/math/__rem_pio2.c +++ b/src/math/__rem_pio2.c @@ -14,7 +14,7 @@ /* __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 */ #include "libm.h" @@ -164,7 +164,7 @@ medium: tx[2] = z; nx = 3; while (tx[nx-1] == zero) nx--; /* skip zero term */ - n = __rem_pio2_slow(tx,ty,e0,nx,1); + n = __rem_pio2_large(tx,ty,e0,nx,1); if (hx < 0) { y[0] = -ty[0]; y[1] = -ty[1]; diff --git a/src/math/__rem_pio2_slow.c b/src/math/__rem_pio2_large.c similarity index 99% rename from src/math/__rem_pio2_slow.c rename to src/math/__rem_pio2_large.c index 364f80e..35835f8 100644 --- a/src/math/__rem_pio2_slow.c +++ b/src/math/__rem_pio2_large.c @@ -10,10 +10,10 @@ * ==================================================== */ /* - * __rem_pio2_slow(x,y,e0,nx,prec) + * __rem_pio2_large(x,y,e0,nx,prec) * double x[],y[]; int e0,nx,prec; * - * __rem_pio2_slow return the last three digits of N with + * __rem_pio2_large return the last three digits of N with * y = x - N*pi/2 * so that |y| < pi/2. * @@ -276,7 +276,7 @@ one = 1.0, two24 = 1.67772160000000000000e+07, /* 0x41700000, 0x00000000 */ twon24 = 5.96046447753906250000e-08; /* 0x3E700000, 0x00000000 */ -int __rem_pio2_slow(double *x, double *y, int e0, int nx, int prec) +int __rem_pio2_large(double *x, double *y, int e0, int nx, int prec) { int32_t jz,jx,jv,jp,jk,carry,n,iq[20],i,j,k,m,q0,ih; double z,fw,f[20],fq[20],q[20]; diff --git a/src/math/__rem_pio2f.c b/src/math/__rem_pio2f.c index 7fa5029..965dc46 100644 --- a/src/math/__rem_pio2f.c +++ b/src/math/__rem_pio2f.c @@ -17,7 +17,7 @@ * * return the remainder of x rem pi/2 in *y * use double precision for everything except passing x - * use __rem_pio2_slow() for large x + * use __rem_pio2_large() for large x */ #include "libm.h" @@ -68,7 +68,7 @@ int __rem_pio2f(float x, double *y) e0 = (ix>>23) - 150; /* e0 = ilogb(|x|)-23; */ SET_FLOAT_WORD(z, ix - ((int32_t)(e0<<23))); tx[0] = z; - n = __rem_pio2_slow(tx,ty,e0,1,0); + n = __rem_pio2_large(tx,ty,e0,1,0); if (hx < 0) { *y = -ty[0]; return -n; diff --git a/src/math/__rem_pio2l.h b/src/math/__rem_pio2l.h index 3a747b1..130a8f0 100644 --- a/src/math/__rem_pio2l.h +++ b/src/math/__rem_pio2l.h @@ -17,7 +17,7 @@ /* 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) { -- 2.20.1