summaryrefslogtreecommitdiff
path: root/src/div.c
diff options
context:
space:
mode:
authorzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2017-12-11 17:02:05 +0000
committerzimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4>2017-12-11 17:02:05 +0000
commitc17961c02b2897ccb5346e3cc4a9b4a3e12197fb (patch)
tree5e7748752a1c7f4f63b258fdc16f20b082a939a1 /src/div.c
parentb29bec14c7145e6aa072cd4d2f570d9f834a9f83 (diff)
downloadmpfr-c17961c02b2897ccb5346e3cc4a9b4a3e12197fb.tar.gz
get rid of DEBUG macro (but keep corresponding commands in comment to help
understanding of the code and debugging) git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@11955 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'src/div.c')
-rw-r--r--src/div.c40
1 files changed, 0 insertions, 40 deletions
diff --git a/src/div.c b/src/div.c
index 86bbdd51c..3ef62b6d5 100644
--- a/src/div.c
+++ b/src/div.c
@@ -651,21 +651,6 @@ mpfr_div_2 (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
#endif /* !defined(MPFR_GENERIC_ABI) */
-#ifdef DEBUG2
-#define mpfr_mpn_print(ap,n) mpfr_mpn_print3 (ap,n,MPFR_LIMB_ZERO)
-static void
-mpfr_mpn_print3 (mpfr_limb_ptr ap, mp_size_t n, mp_limb_t cy)
-{
- mp_size_t i;
- for (i = 0; i < n; i++)
- printf ("+%lu*2^%lu", (unsigned long) ap[i], (unsigned long)
- (GMP_NUMB_BITS * i));
- if (cy)
- printf ("+2^%lu", (unsigned long) (GMP_NUMB_BITS * n));
- printf ("\n");
-}
-#endif
-
/* check if {ap, an} is zero */
static int
mpfr_mpn_cmpzero (mpfr_limb_ptr ap, mp_size_t an)
@@ -1052,9 +1037,6 @@ mpfr_div (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
else /* k=0: no more dividend limb */
extra_bit = mpfr_mpn_cmpzero (vp, l) == 0;
}
-#ifdef DEBUG
- printf ("extra_bit=%d\n", extra_bit);
-#endif
/* set exponent */
qexp = MPFR_GET_EXP (u) - MPFR_GET_EXP (v) + extra_bit;
@@ -1225,10 +1207,6 @@ mpfr_div (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
/* if Mulders' short division failed, we revert to division with remainder */
qh = mpn_divrem (qp, 0, ap + k, qqsize - k, bp, qsize - k);
/* warning: qh may be 1 if u1 == v1, but u < v */
-#ifdef DEBUG2
- printf ("q="); mpfr_mpn_print (qp, qsize);
- printf ("r="); mpfr_mpn_print (ap, qsize);
-#endif
k = qsize;
sticky_u = sticky_u || mpfr_mpn_cmpzero (ap, k);
@@ -1256,20 +1234,12 @@ mpfr_div (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
including the round bit, and 1 <= sh2 <= GMP_NUMB_BITS
is the number of bits in sticky3 */
inex = (sticky != MPFR_LIMB_ZERO) || (sticky3 != MPFR_LIMB_ZERO);
-#ifdef DEBUG
- printf ("sticky=%lu sticky3=%lu inex=%d\n",
- (unsigned long) sticky, (unsigned long) sticky3, inex);
-#endif
/* to round, we distinguish two cases:
(a) vsize <= qsize: we used the full divisor
(b) vsize > qsize: the divisor was truncated
*/
-#ifdef DEBUG
- printf ("vsize=%lu qsize=%lu\n",
- (unsigned long) vsize, (unsigned long) qsize);
-#endif
if (MPFR_LIKELY(vsize <= qsize)) /* use the full divisor */
{
if (MPFR_LIKELY(rnd_mode == MPFR_RNDN))
@@ -1302,10 +1272,6 @@ mpfr_div (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
{
round_bit = sticky3 & (MPFR_LIMB_ONE << (sh2 - 1));
sticky3 = sticky3 ^ round_bit;
-#ifdef DEBUG
- printf ("rb=%lu sb=%lu\n",
- (unsigned long) round_bit, (unsigned long) sticky3);
-#endif
}
if (sticky3 != MPFR_LIMB_ZERO && sticky3 != MPFR_LIMB_ONE)
{
@@ -1343,9 +1309,6 @@ mpfr_div (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
mpfr_mpn_cmp_aux (sp, k, up, usize - qqsize, extra_bit) :
mpfr_mpn_cmpzero (sp, k);
}
-#ifdef DEBUG
- printf ("cmp(q*v0,r+u0)=%d\n", cmp_s_r);
-#endif
/* now cmp_s_r > 0 if {sp, vsize} > {ap, qsize} + low(u)
cmp_s_r = 0 if {sp, vsize} = {ap, qsize} + low(u)
cmp_s_r < 0 if {sp, vsize} < {ap, qsize} + low(u) */
@@ -1403,9 +1366,6 @@ mpfr_div (mpfr_ptr q, mpfr_srcptr u, mpfr_srcptr v, mpfr_rnd_t rnd_mode)
cmp_s_r = 1; /* since in fact we subtracted
less than 1 */
}
-#ifdef DEBUG
- printf ("cmp(q*v0-(r+u0),v)=%d\n", cmp_s_r);
-#endif
if (cmp_s_r <= 0) /* q1-1 <= u/v < q1 */
{
if (sticky3 == MPFR_LIMB_ONE)