diff options
| author | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2003-05-22 21:39:40 +0000 |
|---|---|---|
| committer | vlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4> | 2003-05-22 21:39:40 +0000 |
| commit | 44b4dd94bb98c8d9e7850ae401232bd1b2ea3028 (patch) | |
| tree | 9670f0ef8017d42ad2a2062dc08c63c022e450c8 /sqrt.c | |
| parent | 2f3cb289a102043a22bd32c5950db37199fb3fd2 (diff) | |
| download | mpfr-44b4dd94bb98c8d9e7850ae401232bd1b2ea3028.tar.gz | |
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
Diffstat (limited to 'sqrt.c')
| -rw-r--r-- | sqrt.c | 12 |
1 files changed, 6 insertions, 6 deletions
@@ -1,6 +1,6 @@ /* mpfr_sqrt -- square root of a floating-point number -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. @@ -92,7 +92,7 @@ mpfr_sqrt (mpfr_ptr r, mpfr_srcptr u, mp_rnd_t rnd_mode) /* Compare the mantissas */ - odd_exp_u = (unsigned int) MPFR_EXP(u) & 1; + odd_exp_u = (unsigned int) MPFR_GET_EXP (u) & 1; MPFR_ASSERTN(MPFR_PREC(r) <= MPFR_INTPREC_MAX - 3); rrsize = (MPFR_PREC(r) + 2 + odd_exp_u) / BITS_PER_MP_LIMB + 1; MPFR_ASSERTN(rrsize <= MP_SIZE_T_MAX/2); @@ -124,9 +124,9 @@ mpfr_sqrt (mpfr_ptr r, mpfr_srcptr u, mp_rnd_t rnd_mode) } } - MPFR_EXP(r) = MPFR_EXP(u) != MPFR_EMAX_MAX - ? (MPFR_EXP(u) + odd_exp_u) / 2 - : (MPFR_EMAX_MAX - 1) / 2 + 1; + MPFR_SET_EXP(r, MPFR_GET_EXP(u) != MPFR_EMAX_MAX + ? (MPFR_GET_EXP(u) + odd_exp_u) / 2 + : (MPFR_EMAX_MAX - 1) / 2 + 1); do { @@ -264,7 +264,7 @@ mpfr_sqrt (mpfr_ptr r, mpfr_srcptr u, mp_rnd_t rnd_mode) /* Is a shift necessary here? Isn't the result 1.0000...? */ mpn_rshift (rp, rp, rrsize, 1); rp[rrsize-1] |= MPFR_LIMB_HIGHBIT; - MPFR_EXP(r)++; + MPFR_SET_EXP (r, MPFR_GET_EXP (r) + 1); } fin: |