[Gmp-commit] /home/hgfiles/gmp: Simplified calculation of block size n.

[mercurial at gmplib.org]

details:   /home/hgfiles/gmp/rev/4587ae5f3743
changeset: 13124:4587ae5f3743
user:      Niels M?ller <nisse at lysator.liu.se>
date:      Fri Dec 18 12:07:10 2009 +0100
description:
Simplified calculation of block size n.

diffstat:

 ChangeLog                |  3 +++
 mpn/generic/toom63_mul.c |  2 +-
 2 files changed, 4 insertions(+), 1 deletions(-)

diffs (24 lines):

diff -r 1c1ae1bee11d -r 4587ae5f3743 ChangeLog
--- a/ChangeLog	Fri Dec 18 11:41:41 2009 +0100
+++ b/ChangeLog	Fri Dec 18 12:07:10 2009 +0100
@@ -1,5 +1,8 @@
 2009-12-18  Niels Möller  <nisse at lysator.liu.se>
 
+	* mpn/generic/toom63_mul.c (mpn_toom63_mul): Simplified
+	calculation of block size n.
+
 	* mpn/generic/toom_eval_pm2exp.c (mpn_toom_eval_pm2exp): Fixed
 	output asserts.
 
diff -r 1c1ae1bee11d -r 4587ae5f3743 mpn/generic/toom63_mul.c
--- a/mpn/generic/toom63_mul.c	Fri Dec 18 11:41:41 2009 +0100
+++ b/mpn/generic/toom63_mul.c	Fri Dec 18 12:07:10 2009 +0100
@@ -137,7 +137,7 @@
 #define b2  (bp + 2 * n)
 
   ASSERT (an >= bn);
-  n = 1 + (3 * an >= 6 * bn ? (an - 1) / (size_t) 6 : (bn - 1) / (size_t) 3);
+  n = 1 + (an >= 2 * bn ? (an - 1) / (size_t) 6 : (bn - 1) / (size_t) 3);
 
   s = an - 5 * n;
   t = bn - 2 * n;