diff options
Diffstat (limited to 'gcc')
-rw-r--r-- | gcc/fortran/ChangeLog | 5 | ||||
-rw-r--r-- | gcc/fortran/arith.c | 6 |
2 files changed, 9 insertions, 2 deletions
diff --git a/gcc/fortran/ChangeLog b/gcc/fortran/ChangeLog index 6e5584a01a4..f0455116236 100644 --- a/gcc/fortran/ChangeLog +++ b/gcc/fortran/ChangeLog @@ -1,3 +1,8 @@ +2006-10-03 Steven G. Kargl <kargl@gcc.gnu.org> + + * arith.c (gfc_check_real_range): Use correct exponent range for + subnormal numbers. + 2006-10-03 Paul Thomas <pault@gcc.gnu.org> PR fortran/29284 diff --git a/gcc/fortran/arith.c b/gcc/fortran/arith.c index b473fdb15ee..3541adc49a6 100644 --- a/gcc/fortran/arith.c +++ b/gcc/fortran/arith.c @@ -438,14 +438,16 @@ gfc_check_real_range (mpfr_t p, int kind) gfc_free (bin); #else mp_exp_t emin, emax; + int en; /* Save current values of emin and emax. */ emin = mpfr_get_emin (); emax = mpfr_get_emax (); /* Set emin and emax for the current model number. */ - mpfr_set_emin ((mp_exp_t) gfc_real_kinds[i].min_exponent - 1); - mpfr_set_emax ((mp_exp_t) gfc_real_kinds[i].max_exponent - 1); + en = gfc_real_kinds[i].min_exponent - gfc_real_kinds[i].digits + 1; + mpfr_set_emin ((mp_exp_t) en); + mpfr_set_emax ((mp_exp_t) gfc_real_kinds[i].max_exponent); mpfr_subnormalize (q, 0, GFC_RND_MODE); /* Reset emin and emax. */ |