rename __rem_pio2_slow to __rem_pio2_large
authornsz <nsz@port70.net>
Mon, 5 Mar 2012 03:11:38 +0000 (04:11 +0100)
committernsz <nsz@port70.net>
Mon, 5 Mar 2012 03:11:38 +0000 (04:11 +0100)
src/math/__rem_pio2.c
src/math/__rem_pio2_large.c [moved from src/math/__rem_pio2_slow.c with 99% similarity]
src/math/__rem_pio2f.c
src/math/__rem_pio2l.h

index 4373e22..a7d779e 100644 (file)
@@ -14,7 +14,7 @@
 /* __rem_pio2(x,y)
  *
  * return the remainder of x rem pi/2 in y[0]+y[1]
 /* __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"
  */
 
 #include "libm.h"
@@ -164,7 +164,7 @@ medium:
        tx[2] = z;
        nx = 3;
        while (tx[nx-1] == zero) nx--;  /* skip zero term */
        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];
        if (hx < 0) {
                y[0] = -ty[0];
                y[1] = -ty[1];
similarity index 99%
rename from src/math/__rem_pio2_slow.c
rename to src/math/__rem_pio2_large.c
index 364f80e..35835f8 100644 (file)
  * ====================================================
  */
 /*
  * ====================================================
  */
 /*
- * __rem_pio2_slow(x,y,e0,nx,prec)
+ * __rem_pio2_large(x,y,e0,nx,prec)
  * double x[],y[]; int 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.
  *
  *              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 */
 
 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];
 {
        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];
index 7fa5029..965dc46 100644 (file)
@@ -17,7 +17,7 @@
  *
  * return the remainder of x rem pi/2 in *y
  * use double precision for everything except passing x
  *
  * 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"
  */
 
 #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;
        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;
        if (hx < 0) {
                *y = -ty[0];
                return -n;
index 3a747b1..130a8f0 100644 (file)
@@ -17,7 +17,7 @@
 /* ld80 version of __rem_pio2(x,y)
  *
  * return the remainder of x rem pi/2 in y[0]+y[1]
 /* 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)
  */
 
 #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 */
        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) {
        r = (long double)ty[0] + ty[1];
        w = ty[1] - (r - ty[0]);
        if (expsign < 0) {