about summary refs log tree commit diff
path: root/sysdeps/ieee754/dbl-64/branred.c
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>2001-03-12 07:57:09 +0000
committerUlrich Drepper <drepper@redhat.com>2001-03-12 07:57:09 +0000
commit50944bca4bfeefc5de2ec205ad549eb669ed1008 (patch)
tree8b9a7141aef238c00893d27a39d3a34247374a82 /sysdeps/ieee754/dbl-64/branred.c
parent445028e348508c6bcf7cc41ffab4af1b984a11d2 (diff)
downloadglibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.tar.gz
glibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.tar.xz
glibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.zip
Fix warnings.
Diffstat (limited to 'sysdeps/ieee754/dbl-64/branred.c')
-rw-r--r--sysdeps/ieee754/dbl-64/branred.c53
1 files changed, 29 insertions, 24 deletions
diff --git a/sysdeps/ieee754/dbl-64/branred.c b/sysdeps/ieee754/dbl-64/branred.c
index 09b55602bf..2bceb947c0 100644
--- a/sysdeps/ieee754/dbl-64/branred.c
+++ b/sysdeps/ieee754/dbl-64/branred.c
@@ -5,9 +5,9 @@
  *
  * This program is free software; you can redistribute it and/or modify
  * it under the terms of the GNU Lesser General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or 
+ * the Free Software Foundation; either version 2 of the License, or
  * (at your option) any later version.
- * 
+ *
  * This program is distributed in the hope that it will be useful,
  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
@@ -15,14 +15,14 @@
  *
  * You should have received a copy of the GNU Lesser General Public License
  * along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.  
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
  */
 /*******************************************************************/
 /*                                                                 */
-/* MODULE_NAME: branred.c                                          */ 
+/* MODULE_NAME: branred.c                                          */
 /*                                                                 */
 /* FUNCTIONS:   branred                                            */
-/*                                                                 */ 
+/*                                                                 */
 /* FILES NEEDED: branred.h mydefs.h endian.h mpa.h                 */
 /*               mha.c                                             */
 /*                                                                 */
@@ -45,35 +45,41 @@
 /* x=n*pi/2+(a+aa), abs(a+aa)<pi/4, n=0,+-1,+-2,....               */
 /* Routine return integer (n mod 4)                                */
 /*******************************************************************/
-int branred(double x, double *a, double *aa) 
+int branred(double x, double *a, double *aa)
 {
-  int i,k,n;
-  mynumber  u,v,gor;
-  double r[6],s,t,sum,b,bb,sum1,sum2,b1,bb1,b2,bb2,x1,x2,t1,t2; 
-  
+  int i,k;
+#if 0
+  int n;
+#endif
+  mynumber  u,gor;
+#if 0
+  mynumber v;
+#endif
+  double r[6],s,t,sum,b,bb,sum1,sum2,b1,bb1,b2,bb2,x1,x2,t1,t2;
+
   x*=tm600.x;
   t=x*split;   /* split x to two numbers */
   x1=t-(t-x);
   x2=x-x1;
   sum=0;
   u.x = x1;
-  k = (u.i[HIGH_HALF]>>20)&2047; 
+  k = (u.i[HIGH_HALF]>>20)&2047;
   k = (k-450)/24;
   if (k<0)
     k=0;
   gor.x = t576.x;
   gor.i[HIGH_HALF] -= ((k*24)<<20);
-  for (i=0;i<6;i++) 
+  for (i=0;i<6;i++)
     { r[i] = x1*toverp[k+i]*gor.x; gor.x *= tm24.x; }
   for (i=0;i<3;i++) {
-    s=(r[i]+big.x)-big.x; 
+    s=(r[i]+big.x)-big.x;
     sum+=s;
     r[i]-=s;
   }
   t=0;
-  for (i=0;i<6;i++) 
+  for (i=0;i<6;i++)
     t+=r[5-i];
-  bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5]; 
+  bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5];
   s=(t+big.x)-big.x;
   sum+=s;
   t-=s;
@@ -85,7 +91,7 @@ int branred(double x, double *a, double *aa)
   bb1=bb;
   sum1=sum;
   sum=0;
-  
+
   u.x = x2;
   k = (u.i[HIGH_HALF]>>20)&2047;
   k = (k-450)/24;
@@ -96,14 +102,14 @@ int branred(double x, double *a, double *aa)
   for (i=0;i<6;i++)
     { r[i] = x2*toverp[k+i]*gor.x; gor.x *= tm24.x; }
   for (i=0;i<3;i++) {
-    s=(r[i]+big.x)-big.x; 
+    s=(r[i]+big.x)-big.x;
     sum+=s;
     r[i]-=s;
   }
   t=0;
   for (i=0;i<6;i++)
     t+=r[5-i];
-  bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5]; 
+  bb=(((((r[0]-t)+r[1])+r[2])+r[3])+r[4])+r[5];
   s=(t+big.x)-big.x;
  sum+=s;
  t-=s;
@@ -111,17 +117,17 @@ int branred(double x, double *a, double *aa)
  bb=(t-b)+bb;
  s=(sum+big1.x)-big1.x;
  sum-=s;
- 
+
  b2=b;
  bb2=bb;
  sum2=sum;
- 
+
  sum=sum1+sum2;
  b=b1+b2;
- bb = (ABS(b1)>ABS(b2))? (b1-b)+b2 : (b2-b)+b1;  
- if (b > 0.5) 
+ bb = (ABS(b1)>ABS(b2))? (b1-b)+b2 : (b2-b)+b1;
+ if (b > 0.5)
    {b-=1.0; sum+=1.0;}
- else if (b < -0.5) 
+ else if (b < -0.5)
    {b+=1.0; sum-=1.0;}
  s=b+(bb+bb1+bb2);
  t=((b-s)+bb)+(bb1+bb2);
@@ -136,4 +142,3 @@ int branred(double x, double *a, double *aa)
  *aa=t;
  return ((int) sum)&3; /* return quater of unit circle */
 }
-