diff options
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_atan2.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/e_atan2.c | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_atan2.c b/sysdeps/ieee754/dbl-64/e_atan2.c index b77505db36..adf7a0d11b 100644 --- a/sysdeps/ieee754/dbl-64/e_atan2.c +++ b/sysdeps/ieee754/dbl-64/e_atan2.c @@ -374,9 +374,9 @@ static double normalized(double ax,double ay,double y, double z) { int p; mp_no mpx,mpy,mpz,mperr,mpz2,mpt1; p=6; - dbl_mp(ax,&mpx,p); dbl_mp(ay,&mpy,p); dvd(&mpy,&mpx,&mpz,p); - dbl_mp(ue.d,&mpt1,p); mul(&mpz,&mpt1,&mperr,p); - sub(&mpz,&mperr,&mpz2,p); __mp_dbl(&mpz2,&z,p); + __dbl_mp(ax,&mpx,p); __dbl_mp(ay,&mpy,p); __dvd(&mpy,&mpx,&mpz,p); + __dbl_mp(ue.d,&mpt1,p); __mul(&mpz,&mpt1,&mperr,p); + __sub(&mpz,&mperr,&mpz2,p); __mp_dbl(&mpz2,&z,p); return signArctan2(y,z); } /* Fix the sign and return after stage 1 or stage 2 */ @@ -392,10 +392,10 @@ static double atan2Mp(double x,double y,const int pr[]) mp_no mpx,mpy,mpz,mpz1,mpz2,mperr,mpt1; for (i=0; i<MM; i++) { p = pr[i]; - dbl_mp(x,&mpx,p); dbl_mp(y,&mpy,p); + __dbl_mp(x,&mpx,p); __dbl_mp(y,&mpy,p); __mpatan2(&mpy,&mpx,&mpz,p); - dbl_mp(ud[i].d,&mpt1,p); mul(&mpz,&mpt1,&mperr,p); - add(&mpz,&mperr,&mpz1,p); sub(&mpz,&mperr,&mpz2,p); + __dbl_mp(ud[i].d,&mpt1,p); __mul(&mpz,&mpt1,&mperr,p); + __add(&mpz,&mperr,&mpz1,p); __sub(&mpz,&mperr,&mpz2,p); __mp_dbl(&mpz1,&z1,p); __mp_dbl(&mpz2,&z2,p); if (z1==z2) return z1; } |