aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
Diffstat
-rw-r--r--ChangeLog14
-rw-r--r--sysdeps/ieee754/dbl-64/mpa.c33
-rw-r--r--sysdeps/powerpc/powerpc32/power4/fpu/mpa.c47
-rw-r--r--sysdeps/powerpc/powerpc64/power4/fpu/mpa.c47
4 files changed, 14 insertions, 127 deletions
diff --git a/ChangeLog b/ChangeLog
index 2f70639..8c1295f 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 0db1386..b0571f9 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 b52fa60..f167969 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 b52fa60..f167969 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,
generated by cgit v1.1 at 2024-12-28 11:14:28 +0000