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/slowexp.c | |
| parent | 9a656848eaa2f9c96ce438eeb3c63e33933c0b2e (diff) | |
| download | glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.gz glibc-ca58f1dbeb62840dad345d6bfcca18c81db130a8.tar.bz2 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/slowexp.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/slowexp.c | 20 |
1 files changed, 10 insertions, 10 deletions
diff --git a/sysdeps/ieee754/dbl-64/slowexp.c b/sysdeps/ieee754/dbl-64/slowexp.c index 1c2779b0e9..3a3758b305 100644 --- a/sysdeps/ieee754/dbl-64/slowexp.c +++ b/sysdeps/ieee754/dbl-64/slowexp.c @@ -30,10 +30,10 @@ /**************************************************************************/ #include "mpa.h" -void mpexp(mp_no *x, mp_no *y, int p); +void __mpexp(mp_no *x, mp_no *y, int p); /*Converting from double precision to Multi-precision and calculating e^x */ -double slowexp(double x) { +double __slowexp(double x) { double w,z,res,eps=3.0e-26; #if 0 double y; @@ -45,20 +45,20 @@ double slowexp(double x) { mp_no mpx, mpy, mpz,mpw,mpeps,mpcor; p=6; - dbl_mp(x,&mpx,p); /* Convert a double precision number x */ + __dbl_mp(x,&mpx,p); /* Convert a double precision number x */ /* into a multiple precision number mpx with prec. p. */ - mpexp(&mpx, &mpy, p); /* Multi-Precision exponential function */ - dbl_mp(eps,&mpeps,p); - mul(&mpeps,&mpy,&mpcor,p); - add(&mpy,&mpcor,&mpw,p); - sub(&mpy,&mpcor,&mpz,p); + __mpexp(&mpx, &mpy, p); /* Multi-Precision exponential function */ + __dbl_mp(eps,&mpeps,p); + __mul(&mpeps,&mpy,&mpcor,p); + __add(&mpy,&mpcor,&mpw,p); + __sub(&mpy,&mpcor,&mpz,p); __mp_dbl(&mpw, &w, p); __mp_dbl(&mpz, &z, p); if (w == z) return w; else { /* if calculating is not exactly */ p = 32; - dbl_mp(x,&mpx,p); - mpexp(&mpx, &mpy, p); + __dbl_mp(x,&mpx,p); + __mpexp(&mpx, &mpy, p); __mp_dbl(&mpy, &res, p); return res; } |