From 981cb1a2a30b2fb2844a14e1b39481cda22cce2c Mon Sep 17 00:00:00 2001
From: zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>
Date: Fri, 24 Sep 2004 07:21:04 +0000
Subject: fixed double-rounding bug in get_ui

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@2986 280ebfd0-de03-0410-8827-d642c229c3f4
---
 get_ui.c | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)

(limited to 'get_ui.c')

diff --git a/get_ui.c b/get_ui.c
index 4ed644dae..b926c03d0 100644
--- a/get_ui.c
+++ b/get_ui.c
@@ -34,9 +34,32 @@ mpfr_get_ui (mpfr_srcptr f, mp_rnd_t rnd)
   if (!mpfr_fits_ulong_p (f, rnd) || MPFR_IS_ZERO(f))
     return (unsigned long) 0;
 
+  /* now 0 <= f and either f <= ULONG_MAX or
+     round(f, prec(ULONG_MAX), rnd) <= ULONG_MAX */
+
+  exp = MPFR_EXP(f);
+  if (exp <= 0) /* 0 < f < 1 */
+    {
+      /* return 0 for round to zero or -infinity,
+                1 for round to +infinity,
+                and for round to nearest:
+                   0 if f <= 1/2 (even rule)
+                   1 otherwise */
+      return (rnd == GMP_RNDZ || rnd == GMP_RNDD) ? 0
+        : ((rnd == GMP_RNDU) ? 1
+           : (mpfr_cmp_ui_2exp (f, 1, -1) <= 0 ? 0 : 1));
+    }
+
+  /* now exp > 0 */
+
   /* determine prec of unsigned long */
   for (s = ULONG_MAX, prec = 0; s != 0; s /= 2, prec ++);
 
+  /* since f < 2^EXP(f), we need at most min(prec,EXP(f)) bits */
+  prec = (exp < prec) ? exp : prec;
+  if (prec < MPFR_PREC_MIN)
+    prec = MPFR_PREC_MIN;
+
   /* first round to prec bits */
   mpfr_init2 (x, prec);
   mpfr_set (x, f, rnd);
-- 
cgit v1.2.1