aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/dbl-64/branred.c
diff options
context:
space:
mode:
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 */
}
-