From 01ae893012ce58112229ac1e3bd1050d5afec72b Mon Sep 17 00:00:00 2001 From: "Steven G. Kargl" Date: Thu, 7 Apr 2005 18:26:37 +0000 Subject: [PATCH] simplify.c (gfc_simplify_exponent): Fix exponent(tiny(x)) From-SVN: r97792 --- gcc/fortran/ChangeLog | 4 ++++ gcc/fortran/simplify.c | 7 +++++++ 2 files changed, 11 insertions(+) diff --git a/gcc/fortran/ChangeLog b/gcc/fortran/ChangeLog index d7cbdc0f16e..e114853a20a 100644 --- a/gcc/fortran/ChangeLog +++ b/gcc/fortran/ChangeLog @@ -1,3 +1,7 @@ +2005-04-07 Steven G. Kargl + + * simplify.c (gfc_simplify_exponent): Fix exponent(tiny(x)) + 2005-04-06 Steven G. Kargl * invoke.texi: Remove documentation of -std=f90 diff --git a/gcc/fortran/simplify.c b/gcc/fortran/simplify.c index c2117148839..add391f195e 100644 --- a/gcc/fortran/simplify.c +++ b/gcc/fortran/simplify.c @@ -967,6 +967,7 @@ gfc_simplify_exp (gfc_expr * x) gfc_expr * gfc_simplify_exponent (gfc_expr * x) { + int i; mpfr_t tmp; gfc_expr *result; @@ -991,6 +992,12 @@ gfc_simplify_exponent (gfc_expr * x) gfc_mpfr_to_mpz (result->value.integer, tmp); + /* The model number for tiny(x) is b**(emin - 1) where b is the base and emin + is the smallest exponent value. So, we need to add 1 if x is tiny(x). */ + i = gfc_validate_kind (x->ts.type, x->ts.kind, false); + if (mpfr_cmp (x->value.real, gfc_real_kinds[i].tiny) == 0) + mpz_add_ui (result->value.integer,result->value.integer, 1); + mpfr_clear (tmp); return range_check (result, "EXPONENT"); -- 2.30.2