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/mpatan2.c | |
| parent | 9a656848eaa2f9c96ce438eeb3c63e33933c0b2e (diff) | |
| download | glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.zip glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.gz glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.bz2 | |
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/mpatan2.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/mpatan2.c | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/sysdeps/ieee754/dbl-64/mpatan2.c b/sysdeps/ieee754/dbl-64/mpatan2.c index e41350e..2d1625b 100644 --- a/sysdeps/ieee754/dbl-64/mpatan2.c +++ b/sysdeps/ieee754/dbl-64/mpatan2.c @@ -37,12 +37,12 @@ #include "mpa.h" -void mpsqrt(mp_no *, mp_no *, int); -void mpatan(mp_no *, mp_no *, int); +void __mpsqrt(mp_no *, mp_no *, int); +void __mpatan(mp_no *, mp_no *, int); /* Multi-Precision Atan2(y,x) function subroutine, for p >= 4. */ /* y=0 is not permitted if x<=0. No error messages are given. */ -void mpatan2(mp_no *y, mp_no *x, mp_no *z, int p) { +void __mpatan2(mp_no *y, mp_no *x, mp_no *z, int p) { static const double ZERO = 0.0, ONE = 1.0; @@ -54,15 +54,15 @@ void mpatan2(mp_no *y, mp_no *x, mp_no *z, int p) { if (X[0] <= ZERO) { mpone.e = 1; mpone.d[0] = mpone.d[1] = ONE; - dvd(x,y,&mpt1,p); mul(&mpt1,&mpt1,&mpt2,p); + __dvd(x,y,&mpt1,p); __mul(&mpt1,&mpt1,&mpt2,p); if (mpt1.d[0] != ZERO) mpt1.d[0] = ONE; - add(&mpt2,&mpone,&mpt3,p); mpsqrt(&mpt3,&mpt2,p); - add(&mpt1,&mpt2,&mpt3,p); mpt3.d[0]=Y[0]; - mpatan(&mpt3,&mpt1,p); add(&mpt1,&mpt1,z,p); + __add(&mpt2,&mpone,&mpt3,p); __mpsqrt(&mpt3,&mpt2,p); + __add(&mpt1,&mpt2,&mpt3,p); mpt3.d[0]=Y[0]; + __mpatan(&mpt3,&mpt1,p); __add(&mpt1,&mpt1,z,p); } else - { dvd(y,x,&mpt1,p); - mpatan(&mpt1,z,p); + { __dvd(y,x,&mpt1,p); + __mpatan(&mpt1,z,p); } return; |