From: Paul Brook <paul@codesourcery.com>
Date: Wed, 25 Aug 2004 21:04:49 +0000 (+0000)
Subject: re PR fortran/17190 (MPFR semantics for mpfr_get_z_exp changed)
X-Git-Url: https://git.libre-soc.org/?a=commitdiff_plain;h=9ca2b0db5726c53623f356754e5cc46985e7eaf8;p=gcc.git

re PR fortran/17190 (MPFR semantics for mpfr_get_z_exp changed)

	PR fortran/17190
	* arith.c (gfc_mpfr_to_mpz): Workaround mpfr bug.

From-SVN: r86581
---

diff --git a/gcc/fortran/ChangeLog b/gcc/fortran/ChangeLog
index 85be1029682..dd6b48ea65a 100644
--- a/gcc/fortran/ChangeLog
+++ b/gcc/fortran/ChangeLog
@@ -1,3 +1,8 @@
+2004-08-25  Paul Brook  <paul@codesourcery.com>
+
+	PR fortran/17190
+	* arith.c (gfc_mpfr_to_mpz): Workaround mpfr bug.
+
 2004-08-25  Paul Brook  <paul@codesourcery.com>
 
 	PR fortran/17144
diff --git a/gcc/fortran/arith.c b/gcc/fortran/arith.c
index 03ee14c0999..5f558139401 100644
--- a/gcc/fortran/arith.c
+++ b/gcc/fortran/arith.c
@@ -106,17 +106,20 @@ int gfc_index_integer_kind;
    It's easily implemented with a few calls though.  */
 
 void
-gfc_mpfr_to_mpz(mpz_t z, mpfr_t x)
+gfc_mpfr_to_mpz (mpz_t z, mpfr_t x)
 {
   mp_exp_t e;
 
   e = mpfr_get_z_exp (z, x);
+  /* MPFR 2.0.1 (included with GMP 4.1) has a bug whereby mpfr_get_z_exp
+     may set the sign of z incorrectly.  Work around that here.  */
+  if (mpfr_sgn (x) != mpz_sgn (z))
+    mpz_neg (z, z);
+
   if (e > 0)
     mpz_mul_2exp (z, z, e);
   else
     mpz_tdiv_q_2exp (z, z, -e);
-  if (mpfr_sgn (x) < 0)
-    mpz_neg (z, z);
 }