From bcb5ecae93d5aafdf871668bc0903f974fe4631c Mon Sep 17 00:00:00 2001
From: skimo <skimo>
Date: Wed, 17 Nov 2004 13:42:32 +0000
Subject: [PATCH] multiply with correct coefficient 1/d

---
 barvinok.cc | 8 +++-----
 1 file changed, 3 insertions(+), 5 deletions(-)

diff --git a/barvinok.cc b/barvinok.cc
index 723e239..91286f6 100644
--- a/barvinok.cc
+++ b/barvinok.cc
@@ -438,15 +438,13 @@ struct dpoly_r {
     }
     dpoly_r *div(dpoly& d) {
 	dpoly_r *rc = new dpoly_r(len, dim);
-	rc->denom = power(d.coeff[0], len+1);
-	ZZ cur_d = rc->denom;
+	rc->denom = power(d.coeff[0], len);
+	ZZ inv_d = rc->denom / d.coeff[0];
 	ZZ coeff;
 
 	for (int i = 0; i < len; ++i) {
-	    cur_d /= d.coeff[0];
-
 	    for (int k = 0; k < c[i].size(); ++k) {
-		coeff = c[i][k]->coeff * cur_d;
+		coeff = c[i][k]->coeff * inv_d;
 		rc->add_term(i, c[i][k]->powers, coeff);
 	    }
 
-- 
2.11.4.GIT