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, 4 insertions, 4 deletions

diff --git a/mul.c b/mul.c
index 66dcaf8a9..0262f6de6 100644
--- a/mul.c
+++ b/mul.c
@@ -1,6 +1,6 @@
 /* mpfr_mul -- multiply two floating-point numbers
 
-Copyright 1999, 2000, 2001, 2002 Free Software Foundation, Inc.
+Copyright 1999, 2000, 2001, 2002, 2003 Free Software Foundation, Inc.
 
 This file is part of the MPFR Library.
 
@@ -88,8 +88,8 @@ mpfr_mul (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mp_rnd_t rnd_mode)
       MPFR_RET(0); /* 0 * 0 is exact */
     }
 
-  bx = MPFR_EXP(b);
-  cx = MPFR_EXP(c);
+  bx = MPFR_GET_EXP (b);
+  cx = MPFR_GET_EXP (c);
   /* Note: the exponent of the exact result will be e = bx + cx + ec with
      ec in {-1,0,1} and the following assumes that e is representable. */
   MPFR_ASSERTN(MPFR_EMAX_MAX <= (MP_EXP_T_MAX >> 1));
@@ -155,7 +155,7 @@ mpfr_mul (mpfr_ptr a, mpfr_srcptr b, mpfr_srcptr c, mp_rnd_t rnd_mode)
         rnd_mode = GMP_RNDZ;
       return mpfr_set_underflow (a, rnd_mode, sign_product);
     }
-  MPFR_EXP(a) = ax;
+  MPFR_SET_EXP (a, ax);
 
   if (MPFR_SIGN(a) != sign_product)
     MPFR_CHANGE_SIGN(a);