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/mplog.c | |
| parent | 9a656848eaa2f9c96ce438eeb3c63e33933c0b2e (diff) | |
| download | glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.gz | |
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/mplog.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/mplog.c | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/sysdeps/ieee754/dbl-64/mplog.c b/sysdeps/ieee754/dbl-64/mplog.c index 5957c43cc2..cca0c7e4ed 100644 --- a/sysdeps/ieee754/dbl-64/mplog.c +++ b/sysdeps/ieee754/dbl-64/mplog.c @@ -37,9 +37,9 @@ #include "endian.h" #include "mpa.h" -void mpexp(mp_no *, mp_no *, int); +void __mpexp(mp_no *, mp_no *, int); -void mplog(mp_no *x, mp_no *y, int p) { +void __mplog(mp_no *x, mp_no *y, int p) { #include "mplog.h" int i,m; #if 0 @@ -58,14 +58,14 @@ void mplog(mp_no *x, mp_no *y, int p) { /* Perform m newton iterations to solve for y: exp(y)-x=0. */ /* The iterations formula is: y(n+1)=y(n)+(x*exp(-y(n))-1). */ - cpy(y,&mpt1,p); + __cpy(y,&mpt1,p); for (i=0; i<m; i++) { mpt1.d[0]=-mpt1.d[0]; - mpexp(&mpt1,&mpt2,p); - mul(x,&mpt2,&mpt1,p); - sub(&mpt1,&mpone,&mpt2,p); - add(y,&mpt2,&mpt1,p); - cpy(&mpt1,y,p); + __mpexp(&mpt1,&mpt2,p); + __mul(x,&mpt2,&mpt1,p); + __sub(&mpt1,&mpone,&mpt2,p); + __add(y,&mpt2,&mpt1,p); + __cpy(&mpt1,y,p); } return; } |