Toom-Cook 4x3 and 5x2 multiplication: Toom3.5
Huseyin HISIL
huseyinhisil at iyte.edu.tr
Fri Nov 24 06:30:49 CET 2006
Dear All,
My apologies if I disturb your current concentration but
has anyone ever considered Karatsuba 3-way multiplication
in this fast multiplication context? Any idea? I guess
there should be upper and lower thresholds where it is
faster than both karatsuba and toomcook 3way. (More
multiplications and less overhead). Any comment is
appreciated.
You may want to check
http://arf.iyte.edu.tr/~hhisil/index/main/research/karatsuba_3way.txt
for a detailed explanation.
Regards,
Huseyin Hisil.
========================================
Faculty of Information Technology,
Queensland University of Technology
MS-712, 126 Margaret St
Brisbane Queensland 4001
AUSTRALIA
Phn: +61 7 3864 9557
Fax: +61 7 3221 2384
_--_|\
/. . . QUT
<| , , , |>
\ _.--._/
v
========================================
On Thu, 23 Nov 2006 10:28:21 +0100
"Marco Bodrato" <bodrato at mail.dm.unipi.it> wrote:
> Dear Paul,
>
> On Thu, 23 Nov 2006 09:20:29 +0100, Paul Zimmermann
>wrote
>> > Here are some timings: [...]
>>
>> These timings are great!
>
> On the web page you will also find graphs.
>
>> > - Toom3.5 will be faster than Toom4, for any size, and
>>this is obvious...
>>
>> Sorry it is not obvious to me (and maybe to other
>>readers of gmp-devel).
>> Please can you elaborate?
>
> Shortly. Toom4 now implements only 4x4. Toom3.5 may be
>used 4x3 or 5x2.
> If you compare 4x4 to 4x3 you will have exactly the same
>block sizes. In Toom4
> we have one more multiplication (which nullify, and can
>be avoided) and much
> more overhead (sums, divisions...). Comparing 5x2 to
>4x3, you have the
> same number of mul, sum, shift, div, but in 5x2 you have
>smaller operands.
> So, if they are all aplicable to some chosen operands,
>we have:
> time(5x2) < time(4x3) < time(4x4).
>Finaly, unbalanced Toom3.5, where applicable, will be
>faster than balanced Toom4.
>
>> Isn't it the contrary? mpz_mul does manage different
>>sizes, and I guess you
>> perform padding in mpz_tc3.
>
> Sorry it is not obvious to me :-) I do not know, I
>stealed your mpz_tc4 code :-)
> I guess you are right.
>
> Regards,
> Marco
>
> --
> http://bodrato.it/
>
> _______________________________________________
> gmp-devel mailing list
> gmp-devel at swox.com
> https://gmplib.org/mailman/listinfo/gmp-devel
More information about the gmp-devel
mailing list