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/mpexp.c | |
| parent | 9a656848eaa2f9c96ce438eeb3c63e33933c0b2e (diff) | |
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/mpexp.c')
| -rw-r--r-- | sysdeps/ieee754/dbl-64/mpexp.c | 28 |
1 files changed, 14 insertions, 14 deletions
diff --git a/sysdeps/ieee754/dbl-64/mpexp.c b/sysdeps/ieee754/dbl-64/mpexp.c index 8b7c1dc4bb..2f0b826605 100644 --- a/sysdeps/ieee754/dbl-64/mpexp.c +++ b/sysdeps/ieee754/dbl-64/mpexp.c @@ -35,7 +35,7 @@ /* Multi-Precision exponential function subroutine (for p >= 4, */ /* 2**(-55) <= abs(x) <= 1024). */ -void mpexp(mp_no *x, mp_no *y, int p) { +void __mpexp(mp_no *x, mp_no *y, int p) { int i,j,k,m,m1,m2,n; double a,b; @@ -75,30 +75,30 @@ void mpexp(mp_no *x, mp_no *y, int p) { } /* Compute s=x*2**(-m). Put result in mps */ - dbl_mp(a,&mpt1,p); - mul(x,&mpt1,&mps,p); + __dbl_mp(a,&mpt1,p); + __mul(x,&mpt1,&mps,p); /* Evaluate the polynomial. Put result in mpt2 */ mpone.e=1; mpone.d[0]=ONE; mpone.d[1]=ONE; mpk.e = 1; mpk.d[0] = ONE; mpk.d[1]=nn[n].d; - dvd(&mps,&mpk,&mpt1,p); - add(&mpone,&mpt1,&mpak,p); + __dvd(&mps,&mpk,&mpt1,p); + __add(&mpone,&mpt1,&mpak,p); for (k=n-1; k>1; k--) { - mul(&mps,&mpak,&mpt1,p); + __mul(&mps,&mpak,&mpt1,p); mpk.d[1]=nn[k].d; - dvd(&mpt1,&mpk,&mpt2,p); - add(&mpone,&mpt2,&mpak,p); + __dvd(&mpt1,&mpk,&mpt2,p); + __add(&mpone,&mpt2,&mpak,p); } - mul(&mps,&mpak,&mpt1,p); - add(&mpone,&mpt1,&mpt2,p); + __mul(&mps,&mpak,&mpt1,p); + __add(&mpone,&mpt1,&mpt2,p); /* Raise polynomial value to the power of 2**m. Put result in y */ for (k=0,j=0; k<m; ) { - mul(&mpt2,&mpt2,&mpt1,p); k++; + __mul(&mpt2,&mpt2,&mpt1,p); k++; if (k==m) { j=1; break; } - mul(&mpt1,&mpt1,&mpt2,p); k++; + __mul(&mpt1,&mpt1,&mpt2,p); k++; } - if (j) cpy(&mpt1,y,p); - else cpy(&mpt2,y,p); + if (j) __cpy(&mpt1,y,p); + else __cpy(&mpt2,y,p); return; } |