fixed bug in mpfr_exp_3 (wrong assert)

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@2921 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 4 insertions, 5 deletions

diff --git a/exp3.c b/exp3.c
index 39eaca77c..57a1f7ad9 100644
--- a/exp3.c
+++ b/exp3.c
@@ -118,11 +118,10 @@ mpfr_exp_rational (mpfr_ptr y, mpz_srcptr p, int r, int m)
 
   diff = mpz_sizeinbase(P[0], 2) - precy;
   expo -= diff;
-  /* since prec_i_have, which overestimates nbits(P[0]), can only exceed
-     precy by nbits(i), we assume diff <= 0 here, but we did not manage
-     to prove it, thus the assertion should be checked. */
-  MPFR_ASSERTN(diff <= 0);
-  mpz_mul_2exp (P[0], P[0], -diff);
+  if (diff > 0)
+    mpz_div_2exp (P[0], P[0], diff);
+  else
+    mpz_mul_2exp (P[0], P[0], -diff);
 
   mpz_tdiv_q (S[0], S[0], P[0]);
   mpfr_set_z (y, S[0], GMP_RNDD);