2013-03-26 Siddhesh Poyarekar <siddhesh@redhat.com>
+ * sysdeps/ieee754/dbl-64/branred.c (__branred): Remove
+ commented code.
+ * sysdeps/ieee754/dbl-64/dosincos.c (__dubsin): Likewise.
+ (__dubcos): Likewise.
+ * sysdeps/ieee754/dbl-64/e_asin.c (__ieee754_asin): Likewise.
+ (__ieee754_acos): Likewise.
+ * sysdeps/ieee754/dbl-64/e_atan2.c (__ieee754_atan2): Likewise.
+ * sysdeps/ieee754/dbl-64/e_exp.c (__ieee754_exp): Likewise.
+ (__exp1): Likewise.
+ * sysdeps/ieee754/dbl-64/e_log.c (__ieee754_log): Likewise.
+ * sysdeps/ieee754/dbl-64/e_pow.c (__ieee754_pow): Likewise.
+ (log1): Likewise.
+ (my_log2): Likewise.
+ (checkint): Likewise.
+ * sysdeps/ieee754/dbl-64/e_remainder.c
+ (__ieee754_remainder): Likewise.
+ * sysdeps/ieee754/dbl-64/s_atan.c (atan): Likewise.
+ * sysdeps/ieee754/dbl-64/s_sin.c (__sin): Likewise.
+ (bsloww): Likewise.
+ * sysdeps/ieee754/dbl-64/s_tan.c (tan): Likewise.
+
* sysdeps/ieee754/dbl-64/mpa-arch.h: New file.
* sysdeps/ieee754/dbl-64/mpa.c (norm): Use MANTISSA_T and
MANTISSA_STORE_T to store computations on mantissa. Use
__branred(double x, double *a, double *aa)
{
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;
#ifndef DLA_FMS
double p,hx,tx,hy,ty,q;
#endif
-#if 0
- double xx,y,yy,z,zz;
-#endif
mynumber u;
int4 k;
#ifndef DLA_FMS
double p,hx,tx,hy,ty,q;
#endif
-#if 0
- double xx,y,yy,z,zz;
-#endif
mynumber u;
int4 k;
u.x=x+big.x;
double x1,x2,xx,s1,s2,res1,p,t,res,r,cor,cc,y,c,z,w[2];
mynumber u,v;
int4 k,m,n;
-#if 0
- int4 nn;
-#endif
u.x = x;
m = u.i[HIGH_HALF];
__ieee754_acos(double x)
{
double x1,x2,xx,s1,s2,res1,p,t,res,r,cor,cc,y,c,z,w[2],eps;
-#if 0
- double fc;
-#endif
mynumber u,v;
int4 k,m,n;
-#if 0
- int4 nn;
-#endif
u.x = x;
m = u.i[HIGH_HALF];
k = 0x7fffffff&m;
__ieee754_atan2(double y,double x) {
int i,de,ux,dx,uy,dy;
-#if 0
- int p;
-#endif
static const int pr[MM]={6,8,10,20,32};
double ax,ay,u,du,u9,ua,v,vv,dv,t1,t2,t3,t7,t8,
z,zz,cor,s1,ss1,s2,ss2;
#ifndef DLA_FMS
double t4,t5,t6;
#endif
-#if 0
- double z1,z2;
-#endif
number num;
-#if 0
- mp_no mperr,mpt1,mpx,mpy,mpz,mpz1,mpz2;
-#endif
static const int ep= 59768832, /* 57*16**5 */
em=-59768832; /* -57*16**5 */
__ieee754_exp(double x) {
double bexp, t, eps, del, base, y, al, bet, res, rem, cor;
mynumber junk1, junk2, binexp = {{0,0}};
-#if 0
- int4 k;
-#endif
int4 i,j,m,n,ex;
double retval;
__exp1(double x, double xx, double error) {
double bexp, t, eps, del, base, y, al, bet, res, rem, cor;
mynumber junk1, junk2, binexp = {{0,0}};
-#if 0
- int4 k;
-#endif
int4 i,j,m,n,ex;
junk1.x = x;
#define M 4
static const int pr[M]={8,10,18,32};
int i,j,n,ux,dx,p;
-#if 0
- int k;
-#endif
double dbl_n,u,p0,q,r0,w,nln2a,luai,lubi,lvaj,lvbj,
sij,ssij,ttij,A,B,B0,y,y1,y2,polI,polII,sa,sb,
t1,t2,t7,t8,t,ra,rb,ww,
SECTION
__ieee754_pow(double x, double y) {
double z,a,aa,error, t,a1,a2,y1,y2;
-#if 0
- double gor=1.0;
-#endif
mynumber u,v;
int k;
int4 qx,qy;
SECTION
log1(double x, double *delta, double *error) {
int i,j,m;
-#if 0
- int n;
-#endif
double uu,vv,eps,nx,e,e1,e2,t,t1,t2,res,add=0;
-#if 0
- double cor;
-#endif
mynumber u,v;
#ifdef BIG_ENDI
mynumber
SECTION
my_log2(double x, double *delta, double *error) {
int i,j,m;
-#if 0
- int n;
-#endif
double uu,vv,eps,nx,e,e1,e2,t,t1,t2,res,add=0;
-#if 0
- double cor;
-#endif
double ou1,ou2,lu1,lu2,ov,lv1,lv2,a,a1,a2;
double y,yy,z,zz,j1,j2,j7,j8;
#ifndef DLA_FMS
checkint(double x) {
union {int4 i[2]; double x;} u;
int k,m,n;
-#if 0
- int l;
-#endif
u.x = x;
m = u.i[HIGH_HALF]&0x7fffffff; /* no sign */
if (m >= 0x7ff00000) return 0; /* x is +/-inf or NaN */
double __ieee754_remainder(double x, double y)
{
double z,d,xx;
-#if 0
- double yy;
-#endif
int4 kx,ky,n,nn,n1,m1,l;
-#if 0
- int4 m;
-#endif
mynumber u,t,w={{0,0}},v={{0,0}},ww={{0,0}},r;
u.x=x;
t.x=y;
#ifndef DLA_FMS
double t4,t5,t6;
#endif
-#if 0
- double y1,y2;
-#endif
int i,ux,dx;
-#if 0
- int p;
-#endif
static const int pr[M]={6,8,10,32};
number num;
-#if 0
- mp_no mpt1,mpx,mpy,mpy1,mpy2,mperr;
-#endif
num.d = x; ux = num.i[HIGH_HALF]; dx = num.i[LOW_HALF];
SECTION
__sin(double x){
double xx,res,t,cor,y,s,c,sn,ssn,cs,ccs,xn,a,da,db,eps,xn1,xn2;
-#if 0
- double w[2];
-#endif
mynumber u,v;
int4 k,m,n;
-#if 0
- int4 nn;
-#endif
double retval = 0;
SET_RESTORE_ROUND_53BIT (FE_TONEAREST);
bsloww(double x,double dx, double orig,int n) {
static const double th2_36 = 206158430208.0; /* 1.5*2**37 */
double y,x1,x2,xx,r,t,res,cor,w[2];
-#if 0
- double a,da,xn;
- union {int4 i[2]; double x;} v;
-#endif
x1=(x+th2_36)-th2_36;
y = aa.x*x1*x1*x1;
r=x+y;
int p;
number num,v;
mp_no mpa,mpt1,mpt2;
-#if 0
- mp_no mpy;
-#endif
double retval;