diff options
| -rw-r--r-- | ChangeLog | 14 | ||||
| -rw-r--r-- | sysdeps/ieee754/dbl-64/mpa.c | 33 | ||||
| -rw-r--r-- | sysdeps/powerpc/powerpc32/power4/fpu/mpa.c | 47 | ||||
| -rw-r--r-- | sysdeps/powerpc/powerpc64/power4/fpu/mpa.c | 47 |
4 files changed, 14 insertions, 127 deletions
diff --git a/ChangeLog b/ChangeLog index 2f706390a1..8c1295ff1c 100644 --- a/ChangeLog +++ b/ChangeLog @@ -1,5 +1,19 @@ 2013-01-04 Siddhesh Poyarekar <siddhesh@redhat.com> + * sysdeps/ieee754/dbl-64/mpa.c: Remove commented code. + * sysdeps/powerpc/powerpc32/power4/fpu/mpa.c (__cr): Remove. + (__cpymn): Likewise. + (norm): Remove commented code. + (denorm): Likewise. + (__mp_dbl): Likewise. + (__inv): Likewise. + * sysdeps/powerpc/powerpc64/power4/fpu/mpa.c (__cr): Remove. + (__cpymn): Likewise. + (norm): Remove commented code. + (denorm): Likewise. + (__mp_dbl): Likewise. + (__inv): Likewise. + * sysdeps/powerpc/powerpc32/power4/fpu/mpa.c (mpone): Define mp_no value for 1.0 and 2.0. (norm): Use RADIXI instead of radixi.d. diff --git a/sysdeps/ieee754/dbl-64/mpa.c b/sysdeps/ieee754/dbl-64/mpa.c index 0db13865fd..b0571f9f93 100644 --- a/sysdeps/ieee754/dbl-64/mpa.c +++ b/sysdeps/ieee754/dbl-64/mpa.c @@ -93,21 +93,6 @@ __acr(const mp_no *x, const mp_no *y, int p) { #endif -#if 0 -/* cr() compares the values of two multiple precision numbers */ -static int __cr(const mp_no *x, const mp_no *y, int p) { - int i; - - if (X[0] > Y[0]) i= 1; - else if (X[0] < Y[0]) i=-1; - else if (X[0] < ZERO ) i= __acr(y,x,p); - else i= __acr(x,y,p); - - return i; -} -#endif - - #ifndef NO___CPY /* Copy a multiple precision number. Set *y=*x. x=y is permissible. */ void __cpy(const mp_no *x, mp_no *y, int p) { @@ -117,24 +102,6 @@ void __cpy(const mp_no *x, mp_no *y, int p) { #endif -#if 0 -/* Copy a multiple precision number x of precision m into a */ -/* multiple precision number y of precision n. In case n>m, */ -/* the digits of y beyond the m'th are set to zero. In case */ -/* n<m, the digits of x beyond the n'th are ignored. */ -/* x=y is permissible. */ - -static void __cpymn(const mp_no *x, int m, mp_no *y, int n) { - - int i,k; - - EY = EX; k=MIN(m,n); - for (i=0; i <= k; i++) Y[i] = X[i]; - for ( ; i <= n; i++) Y[i] = ZERO; -} -#endif - - #ifndef NO___MP_DBL /* Convert a multiple precision number *x into a double precision */ /* number *y, normalized case (|x| >= 2**(-1022))) */ diff --git a/sysdeps/powerpc/powerpc32/power4/fpu/mpa.c b/sysdeps/powerpc/powerpc32/power4/fpu/mpa.c index b52fa60632..f167969ea3 100644 --- a/sysdeps/powerpc/powerpc32/power4/fpu/mpa.c +++ b/sysdeps/powerpc/powerpc32/power4/fpu/mpa.c @@ -23,9 +23,7 @@ /* FUNCTIONS: */ /* mcr */ /* acr */ -/* cr */ /* cpy */ -/* cpymn */ /* norm */ /* denorm */ /* mp_dbl */ @@ -86,19 +84,6 @@ int __acr(const mp_no *x, const mp_no *y, int p) { } -/* cr90 compares the values of two multiple precision numbers */ -int __cr(const mp_no *x, const mp_no *y, int p) { - int i; - - if (X[0] > Y[0]) i= 1; - else if (X[0] < Y[0]) i=-1; - else if (X[0] < ZERO ) i= __acr(y,x,p); - else i= __acr(x,y,p); - - return i; -} - - /* Copy a multiple precision number. Set *y=*x. x=y is permissible. */ void __cpy(const mp_no *x, mp_no *y, int p) { long i; @@ -110,34 +95,12 @@ void __cpy(const mp_no *x, mp_no *y, int p) { } -/* Copy a multiple precision number x of precision m into a */ -/* multiple precision number y of precision n. In case n>m, */ -/* the digits of y beyond the m'th are set to zero. In case */ -/* n<m, the digits of x beyond the n'th are ignored. */ -/* x=y is permissible. */ - -void __cpymn(const mp_no *x, int m, mp_no *y, int n) { - - long i,k; - long n2 = n; - long m2 = m; - - EY = EX; k=MIN(m2,n2); - for (i=0; i <= k; i++) Y[i] = X[i]; - for ( ; i <= n2; i++) Y[i] = ZERO; - - return; -} - /* Convert a multiple precision number *x into a double precision */ /* number *y, normalized case (|x| >= 2**(-1022))) */ static void norm(const mp_no *x, double *y, int p) { #define R RADIXI long i; -#if 0 - int k; -#endif double a,c,u,v,z[5]; if (p<5) { if (p==1) c = X[1]; @@ -191,9 +154,6 @@ static void denorm(const mp_no *x, double *y, int p) long i,k; long p2 = p; double c,u,z[5]; -#if 0 - double a,v; -#endif #define R RADIXI if (EX<-44 || (EX==-44 && X[1]<TWO5)) @@ -238,10 +198,6 @@ static void denorm(const mp_no *x, double *y, int p) /* The result is correctly rounded to the nearest/even. *x is left unchanged */ void __mp_dbl(const mp_no *x, double *y, int p) { -#if 0 - int i,k; - double a,c,u,v,z[5]; -#endif if (X[0] == ZERO) {*y = ZERO; return; } @@ -512,9 +468,6 @@ void __mul(const mp_no *x, const mp_no *y, mp_no *z, int p) { void __inv(const mp_no *x, mp_no *y, int p) { long i; -#if 0 - int l; -#endif double t; mp_no z,w; static const int np1[] = {0,0,0,0,1,2,2,2,2,3,3,3,3,3,3,3,3,3, diff --git a/sysdeps/powerpc/powerpc64/power4/fpu/mpa.c b/sysdeps/powerpc/powerpc64/power4/fpu/mpa.c index b52fa60632..f167969ea3 100644 --- a/sysdeps/powerpc/powerpc64/power4/fpu/mpa.c +++ b/sysdeps/powerpc/powerpc64/power4/fpu/mpa.c @@ -23,9 +23,7 @@ /* FUNCTIONS: */ /* mcr */ /* acr */ -/* cr */ /* cpy */ -/* cpymn */ /* norm */ /* denorm */ /* mp_dbl */ @@ -86,19 +84,6 @@ int __acr(const mp_no *x, const mp_no *y, int p) { } -/* cr90 compares the values of two multiple precision numbers */ -int __cr(const mp_no *x, const mp_no *y, int p) { - int i; - - if (X[0] > Y[0]) i= 1; - else if (X[0] < Y[0]) i=-1; - else if (X[0] < ZERO ) i= __acr(y,x,p); - else i= __acr(x,y,p); - - return i; -} - - /* Copy a multiple precision number. Set *y=*x. x=y is permissible. */ void __cpy(const mp_no *x, mp_no *y, int p) { long i; @@ -110,34 +95,12 @@ void __cpy(const mp_no *x, mp_no *y, int p) { } -/* Copy a multiple precision number x of precision m into a */ -/* multiple precision number y of precision n. In case n>m, */ -/* the digits of y beyond the m'th are set to zero. In case */ -/* n<m, the digits of x beyond the n'th are ignored. */ -/* x=y is permissible. */ - -void __cpymn(const mp_no *x, int m, mp_no *y, int n) { - - long i,k; - long n2 = n; - long m2 = m; - - EY = EX; k=MIN(m2,n2); - for (i=0; i <= k; i++) Y[i] = X[i]; - for ( ; i <= n2; i++) Y[i] = ZERO; - - return; -} - /* Convert a multiple precision number *x into a double precision */ /* number *y, normalized case (|x| >= 2**(-1022))) */ static void norm(const mp_no *x, double *y, int p) { #define R RADIXI long i; -#if 0 - int k; -#endif double a,c,u,v,z[5]; if (p<5) { if (p==1) c = X[1]; @@ -191,9 +154,6 @@ static void denorm(const mp_no *x, double *y, int p) long i,k; long p2 = p; double c,u,z[5]; -#if 0 - double a,v; -#endif #define R RADIXI if (EX<-44 || (EX==-44 && X[1]<TWO5)) @@ -238,10 +198,6 @@ static void denorm(const mp_no *x, double *y, int p) /* The result is correctly rounded to the nearest/even. *x is left unchanged */ void __mp_dbl(const mp_no *x, double *y, int p) { -#if 0 - int i,k; - double a,c,u,v,z[5]; -#endif if (X[0] == ZERO) {*y = ZERO; return; } @@ -512,9 +468,6 @@ void __mul(const mp_no *x, const mp_no *y, mp_no *z, int p) { void __inv(const mp_no *x, mp_no *y, int p) { long i; -#if 0 - int l; -#endif double t; mp_no z,w; static const int np1[] = {0,0,0,0,1,2,2,2,2,3,3,3,3,3,3,3,3,3, |