sqrt algorithm

Marco Bodrato bodrato at mail.dm.unipi.it
Wed Jul 29 22:00:50 UTC 2015


Ciao,

On Wed, July 29, 2015 9:24 pm, Niels Möller wrote:
> It could likely be made to work if the input is properly normalized,
> does the current code do that?

Right, current code does normalise in such a way that at most the highest
bit of the input is zero

> And for small sizes, one ought to sqrlo rather than mulmod_bnm1 to
> exploit the cancellation. So the tuning would be between sqrlo and
> sqrmod_bnm1. But it seems we only have mullo, not sqrlo?

Right, we currently don't have sqrlo...

m

-- 
http://bodrato.it/



More information about the gmp-devel mailing list