From 77b558e2d393bac5f180f4bb77f054bf15a1acf8 Mon Sep 17 00:00:00 2001 From: "Steven G. Kargl" Date: Tue, 3 Oct 2006 21:28:05 +0000 Subject: [PATCH] arith.c (gfc_check_real_range): Use correct exponent range for subnormal numbers. 2006-10-03 Steven G. Kargl * arith.c (gfc_check_real_range): Use correct exponent range for subnormal numbers. From-SVN: r117414 --- gcc/fortran/ChangeLog | 5 +++++ 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 + + * arith.c (gfc_check_real_range): Use correct exponent range for + subnormal numbers. + 2006-10-03 Paul Thomas 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. */ -- 2.30.2