diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/s_atan.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/s_atan.c | 15 |
1 files changed, 11 insertions, 4 deletions
diff --git a/sysdeps/ieee754/dbl-64/s_atan.c b/sysdeps/ieee754/dbl-64/s_atan.c index 997e29b5e4..0b71c0d5c8 100644 --- a/sysdeps/ieee754/dbl-64/s_atan.c +++ b/sysdeps/ieee754/dbl-64/s_atan.c @@ -51,12 +51,19 @@ double atan(double x) { double cor,s1,ss1,s2,ss2,t1,t2,t3,t4,t5,t6,t7,t8,t9,t10,u,u2,u3, - v,vv,w,ww,y,yy,y1,y2,z,zz; - int i,ux,dx,p; + v,vv,w,ww,y,yy,z,zz; +#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]; @@ -210,7 +217,7 @@ for (i=0; i<M; i++) { dbl_mp(x,&mpx,p); __mpatan(&mpx,&mpy,p); dbl_mp(u9[i].d,&mpt1,p); mul(&mpy,&mpt1,&mperr,p); add(&mpy,&mperr,&mpy1,p); sub(&mpy,&mperr,&mpy2,p); - mp_dbl(&mpy1,&y1,p); mp_dbl(&mpy2,&y2,p); + __mp_dbl(&mpy1,&y1,p); __mp_dbl(&mpy2,&y2,p); if (y1==y2) return y1; } return y1; /*if unpossible to do exact computing */ |