[src/gamma.c] Added another assertion.

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

1 files changed, 6 insertions, 1 deletions

diff --git a/src/gamma.c b/src/gamma.c
index 397455df7..0cf171928 100644
--- a/src/gamma.c
+++ b/src/gamma.c
@@ -91,7 +91,12 @@ bits_fac (unsigned long n)
   mpfr_log2 (x, x, MPFR_RNDZ);
   r = mpfr_get_ui (x, MPFR_RNDU);  /* lower bound on ceil(x) */
   for (k = 2; k <= n; k *= 2)
-    r -= n / k;
+    {
+      /* Note: the approximation is accurate enough so that the
+         subtractions do not wrap. */
+      MPFR_ASSERTD (r >= n / k);
+      r -= n / k;
+    }
   mpfr_clear (x);
   mpfr_clear (y);
   MPFR_SAVE_EXPO_FREE (expo);