diff options
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; } |