diff options
author | Ulrich Drepper <drepper@redhat.com> | 2001-03-13 02:01:34 +0000 |
---|---|---|
committer | Ulrich Drepper <drepper@redhat.com> | 2001-03-13 02:01:34 +0000 |
commit | ca58f1dbeb62840dad345d6bfcca18c81db130a8 (patch) | |
tree | eb1b9705fc324e0852875514dda109641e9399de /sysdeps/ieee754/dbl-64/mpatan.c | |
parent | 9a656848eaa2f9c96ce438eeb3c63e33933c0b2e (diff) | |
download | glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.gz glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.xz glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.zip |
Update.
2001-03-12 Ulrich Drepper <drepper@redhat.com> * sysdeps/ieee754/dbl-64/e_remainder.c: Fix handling of boundary conditions. * sysdeps/ieee754/dbl-64/e_pow.c: Fix handling of boundary conditions. * sysdeps/ieee754/dbl-64/s_sin.c (__sin): Handle Inf and NaN correctly. (__cos): Likewise. * sysdeps/ieee754/dbl-64/e_asin.c (__ieee754_asin): Handle NaN correctly. (__ieee754_acos): Likewise. redefinition. * sysdeps/ieee754/dbl-64/endian.h: Define also one of BIG_ENDI and LITTLE_ENDI. * sysdeps/ieee754/dbl-64/MathLib.h (Init_Lib): Use void as parameter list.
Diffstat (limited to 'sysdeps/ieee754/dbl-64/mpatan.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/mpatan.c | 38 |
1 files changed, 19 insertions, 19 deletions
diff --git a/sysdeps/ieee754/dbl-64/mpatan.c b/sysdeps/ieee754/dbl-64/mpatan.c index 416014bab7..0c8556db60 100644 --- a/sysdeps/ieee754/dbl-64/mpatan.c +++ b/sysdeps/ieee754/dbl-64/mpatan.c @@ -33,9 +33,9 @@ #include "endian.h" #include "mpa.h" -void mpsqrt(mp_no *, mp_no *, int); +void __mpsqrt(mp_no *, mp_no *, int); -void mpatan(mp_no *x, mp_no *y, int p) { +void __mpatan(mp_no *x, mp_no *y, int p) { #include "mpatan.h" int i,m,n; @@ -66,36 +66,36 @@ void mpatan(mp_no *x, mp_no *y, int p) { mptwo.d[1] = TWO; /* Reduce x m times */ - mul(x,x,&mpsm,p); - if (m==0) cpy(x,&mps,p); + __mul(x,x,&mpsm,p); + if (m==0) __cpy(x,&mps,p); else { for (i=0; i<m; i++) { - add(&mpone,&mpsm,&mpt1,p); - mpsqrt(&mpt1,&mpt2,p); - add(&mpt2,&mpt2,&mpt1,p); - add(&mptwo,&mpsm,&mpt2,p); - add(&mpt1,&mpt2,&mpt3,p); - dvd(&mpsm,&mpt3,&mpt1,p); - cpy(&mpt1,&mpsm,p); + __add(&mpone,&mpsm,&mpt1,p); + __mpsqrt(&mpt1,&mpt2,p); + __add(&mpt2,&mpt2,&mpt1,p); + __add(&mptwo,&mpsm,&mpt2,p); + __add(&mpt1,&mpt2,&mpt3,p); + __dvd(&mpsm,&mpt3,&mpt1,p); + __cpy(&mpt1,&mpsm,p); } - mpsqrt(&mpsm,&mps,p); mps.d[0] = X[0]; + __mpsqrt(&mpsm,&mps,p); mps.d[0] = X[0]; } /* Evaluate a truncated power series for Atan(s) */ n=np[p]; mptwoim1.d[1] = twonm1[p].d; - dvd(&mpsm,&mptwoim1,&mpt,p); + __dvd(&mpsm,&mptwoim1,&mpt,p); for (i=n-1; i>1; i--) { mptwoim1.d[1] -= TWO; - dvd(&mpsm,&mptwoim1,&mpt1,p); - mul(&mpsm,&mpt,&mpt2,p); - sub(&mpt1,&mpt2,&mpt,p); + __dvd(&mpsm,&mptwoim1,&mpt1,p); + __mul(&mpsm,&mpt,&mpt2,p); + __sub(&mpt1,&mpt2,&mpt,p); } - mul(&mps,&mpt,&mpt1,p); - sub(&mps,&mpt1,&mpt,p); + __mul(&mps,&mpt,&mpt1,p); + __sub(&mps,&mpt1,&mpt,p); /* Compute Atan(x) */ mptwoim1.d[1] = twom[m].d; - mul(&mptwoim1,&mpt,y,p); + __mul(&mptwoim1,&mpt,y,p); return; } |