Wed, 12 Mar 2003 04:07:46 +0000
is a faster powm function eg at 750,000 decimal digits it's about 60% fas=
than mpz_powm and if you have the time even faster for larger ones.
Note: It's a large download 1.7M , as it includes some other things=20
Basically its a mpn based sliding base 2^k with barrett reduction .
For small values it is slower than the existing mpz_powm , all I need to =
add a basecase redc to sort this out , which is very easy.
I chose barrett reduction over redc reduction , because I have other uses=
this for which redc is no good (ie I need to test divisibility and/or siz=
residues), and powm is a good test case.