simplify.c (gfc_simplify_exponent): Fix exponent(tiny(x))

git-svn-id: svn+ssh://gcc.gnu.org/svn/gcc/trunk@97792 138bc75d-0d04-0410-961f-82ee72b054a4

diff --git a/gcc/fortran/ChangeLog b/gcc/fortran/ChangeLog
index d7cbdc0..e114853 100644 (file)
--- a/gcc/fortran/ChangeLog
+++ b/gcc/fortran/ChangeLog
@@ -1,3 +1,7 @@
+2005-04-07  Steven G. Kargl  <kargls@comcast.net>
+
+       * simplify.c (gfc_simplify_exponent): Fix exponent(tiny(x))
+
 2005-04-06  Steven G. Kargl  <kargls@comcast.net>
 
        * invoke.texi: Remove documentation of -std=f90 

diff --git a/gcc/fortran/simplify.c b/gcc/fortran/simplify.c
index c211714..add391f 100644 (file)
--- 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");