From 9ca2b0db5726c53623f356754e5cc46985e7eaf8 Mon Sep 17 00:00:00 2001 From: Paul Brook Date: Wed, 25 Aug 2004 21:04:49 +0000 Subject: [PATCH] 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 --- gcc/fortran/ChangeLog | 5 +++++ gcc/fortran/arith.c | 9 ++++++--- 2 files changed, 11 insertions(+), 3 deletions(-) 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 + + PR fortran/17190 + * arith.c (gfc_mpfr_to_mpz): Workaround mpfr bug. + 2004-08-25 Paul Brook 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); } -- 2.30.2