Macros MPFR_EXP_INVALID (invalid exponent value) and MPFR_EXP_CHECK

added. Code update to use MPFR_GET_EXP and MPFR_SET_EXP instead of
MPFR_EXP to allow more bug detection related to special values.
Macros MPFR_SET_NAN, MPFR_SET_INF, MPFR_SET_ZERO and MPFR_INIT set
the exponent of the number to MPFR_EXP_INVALID if MPFR_EXP_CHECK
is defined. Compile with -DMPFR_EXP_CHECK and make check to see
the potential problems; currently, 40 of 76 tests fail.


git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@2301 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 23 insertions, 22 deletions

diff --git a/asinh.c b/asinh.c
index 8ca8f41d7..93830d2d7 100644
--- a/asinh.c
+++ b/asinh.c
@@ -1,6 +1,6 @@
 /* mpfr_asinh -- inverse hyperbolic sine
 
-Copyright 2001, 2002 Free Software Foundation.
+Copyright 2001, 2002, 2003 Free Software Foundation.
 
 This file is part of the MPFR Library.
 
@@ -82,27 +82,28 @@ mpfr_asinh (mpfr_ptr y, mpfr_srcptr x, mp_rnd_t rnd_mode)
   mpfr_save_emin_emax ();
 
   /* First computation of asinh */
-  do {
-
-    /* reactualisation of the precision */
-    mpfr_set_prec (t, Nt);
-    mpfr_set_prec (te, Nt);
-    mpfr_set_prec (ti, Nt);
-
-    /* compute asinh */
-    mpfr_mul (te, x, x, GMP_RNDD);                     /* x^2 */
-    mpfr_add_ui (ti, te, 1, GMP_RNDD);                 /* x^2+1 */
-    mpfr_sqrt (t, ti, GMP_RNDN);                       /* sqrt(x^2+1) */
-    ((neg) ? mpfr_sub : mpfr_add) (t, t, x, GMP_RNDN); /* sqrt(x^2+1)+x */
-    mpfr_log (t, t, GMP_RNDN);                         /* ln(sqrt(x^2+1)+x)*/
-
-    /* estimation of the error -- see algorithms.ps */
-    err = Nt - (MAX(3 - MPFR_EXP(t), 0) + 1);
-
-    /* actualisation of the precision */
-    Nt += 10;
-
-  } while ((err < 0) || (!mpfr_can_round (t, err, GMP_RNDN, rnd_mode, Ny) || (MPFR_IS_ZERO(t))));
+  do
+    {
+      /* reactualisation of the precision */
+      mpfr_set_prec (t, Nt);
+      mpfr_set_prec (te, Nt);
+      mpfr_set_prec (ti, Nt);
+
+      /* compute asinh */
+      mpfr_mul (te, x, x, GMP_RNDD);                   /* x^2 */
+      mpfr_add_ui (ti, te, 1, GMP_RNDD);               /* x^2+1 */
+      mpfr_sqrt (t, ti, GMP_RNDN);                     /* sqrt(x^2+1) */
+      (neg ? mpfr_sub : mpfr_add) (t, t, x, GMP_RNDN); /* sqrt(x^2+1)+x */
+      mpfr_log (t, t, GMP_RNDN);                       /* ln(sqrt(x^2+1)+x)*/
+
+      /* error estimate -- see algorithms.ps */
+      err = Nt - (MAX(3 - MPFR_GET_EXP (t), 0) + 1);
+
+      /* actualisation of the precision */
+      Nt += 10;
+    }
+  while ((err < 0) || (!mpfr_can_round (t, err, GMP_RNDN, rnd_mode, Ny)
+                       || MPFR_IS_ZERO(t)));
 
   mpfr_restore_emin_emax ();