aboutsummaryrefslogtreecommitdiff
path: root/sysdeps/ieee754/dbl-64/slowexp.c
diff options
context:
space:
mode:
authorUlrich Drepper <drepper@redhat.com>2001-03-12 07:57:09 +0000
committerUlrich Drepper <drepper@redhat.com>2001-03-12 07:57:09 +0000
commit50944bca4bfeefc5de2ec205ad549eb669ed1008 (patch)
tree8b9a7141aef238c00893d27a39d3a34247374a82 /sysdeps/ieee754/dbl-64/slowexp.c
parent445028e348508c6bcf7cc41ffab4af1b984a11d2 (diff)
downloadglibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.tar
glibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.tar.gz
glibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.tar.bz2
glibc-50944bca4bfeefc5de2ec205ad549eb669ed1008.zip
Fix warnings.
Diffstat (limited to 'sysdeps/ieee754/dbl-64/slowexp.c')
-rw-r--r--sysdeps/ieee754/dbl-64/slowexp.c27
1 files changed, 16 insertions, 11 deletions
diff --git a/sysdeps/ieee754/dbl-64/slowexp.c b/sysdeps/ieee754/dbl-64/slowexp.c
index d394dc3438..1c2779b0e9 100644
--- a/sysdeps/ieee754/dbl-64/slowexp.c
+++ b/sysdeps/ieee754/dbl-64/slowexp.c
@@ -5,9 +5,9 @@
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU Lesser General Public License as published by
- * the Free Software Foundation; either version 2 of the License, or
+ * the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
- *
+ *
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
@@ -15,12 +15,12 @@
*
* You should have received a copy of the GNU Lesser General Public License
* along with this program; if not, write to the Free Software
- * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
*/
/**************************************************************************/
/* MODULE_NAME:slowexp.c */
/* */
-/* FUNCTION:slowexp */
+/* FUNCTION:slowexp */
/* */
/* FILES NEEDED:mpa.h */
/* mpa.c mpexp.c */
@@ -34,10 +34,16 @@ 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 y,w,z,res,eps=3.0e-26;
- int orig,i,p;
+ double w,z,res,eps=3.0e-26;
+#if 0
+ double y;
+#endif
+ int p;
+#if 0
+ int orig,i;
+#endif
mp_no mpx, mpy, mpz,mpw,mpeps,mpcor;
-
+
p=6;
dbl_mp(x,&mpx,p); /* Convert a double precision number x */
/* into a multiple precision number mpx with prec. p. */
@@ -46,15 +52,14 @@ double slowexp(double x) {
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);
+ __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);
- mp_dbl(&mpy, &res, p);
+ __mp_dbl(&mpy, &res, p);
return res;
}
}
-