summaryrefslogtreecommitdiff
path: root/log10.c
diff options
context:
space:
mode:
authorvlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2004-01-14 17:18:34 +0000
committervlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2004-01-14 17:18:34 +0000
commitf9106da753829c3adca5c0cf8c38c4608cc34218 (patch)
tree2cbb6e0fc012305dc760d50ec8cb93f5b8e746ce /log10.c
parent0a40b14f2ae6f51268df552f999845f23de4a321 (diff)
downloadmpfr-f9106da753829c3adca5c0cf8c38c4608cc34218.tar.gz
Changed some error messages into assertions.
Removed some useless #include's. git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@2622 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'log10.c')
-rw-r--r--log10.c17
1 files changed, 8 insertions, 9 deletions
diff --git a/log10.c b/log10.c
index ede7ab463..1a6b0b91a 100644
--- a/log10.c
+++ b/log10.c
@@ -1,6 +1,6 @@
/* mpfr_log10 -- logarithm in base 10.
-Copyright 2001, 2002, 2003 Free Software Foundation, Inc.
+Copyright 2001, 2002, 2003, 2004 Free Software Foundation, Inc.
This file is part of the MPFR Library.
@@ -19,7 +19,6 @@ along with the MPFR Library; see the file COPYING.LIB. If not, write to
the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
MA 02111-1307, USA. */
-#include <stdio.h>
#include "gmp.h"
#include "gmp-impl.h"
#include "mpfr.h"
@@ -69,7 +68,7 @@ mpfr_log10 (mpfr_ptr r, mpfr_srcptr a, mp_rnd_t rnd_mode)
else
MPFR_ASSERTN(0);
}
-
+
/* If a is negative, the result is NaN */
if (MPFR_UNLIKELY( MPFR_IS_NEG(a) ))
{
@@ -97,10 +96,10 @@ mpfr_log10 (mpfr_ptr r, mpfr_srcptr a, mp_rnd_t rnd_mode)
/* Declaration of the size variable */
mp_prec_t Nx = MPFR_PREC(a); /* Precision of input variable */
mp_prec_t Ny = MPFR_PREC(r); /* Precision of output variable */
-
+
mp_prec_t Nt; /* Precision of the intermediary variable */
long int err; /* Precision of error */
-
+
/* compute the precision of intermediary variable */
Nt = MAX(Nx, Ny);
/* the optimal number of bits : see algorithms.ps */
@@ -110,14 +109,14 @@ mpfr_log10 (mpfr_ptr r, mpfr_srcptr a, mp_rnd_t rnd_mode)
mpfr_init (t);
mpfr_init (tt);
-
+
/* First computation of log10 */
do {
/* reactualisation of the precision */
mpfr_set_prec (t, Nt);
- mpfr_set_prec (tt, Nt);
-
+ mpfr_set_prec (tt, Nt);
+
/* compute log10 */
mpfr_set_ui (t, 10, GMP_RNDN); /* 10 */
mpfr_log (t, t, GMP_RNDD); /* log(10) */
@@ -140,7 +139,7 @@ mpfr_log10 (mpfr_ptr r, mpfr_srcptr a, mp_rnd_t rnd_mode)
/* actualisation of the precision */
Nt += 10;
} while ((err < 0) || !ok);
-
+
inexact = mpfr_set (r, t, rnd_mode);
mpfr_clear (t);