Anomaly in mpn_sqrtrem and mpn_rottrem

Niels Möller nisse at
Thu Jun 11 09:00:19 UTC 2015

tg at (Torbjörn Granlund) writes:

> Before I try to understand the rest of your reasoning: What is B?  It's
> not the usual limb base, I presume, since then sqrt(B^{n-2}) =
> (B/2)^{n-2}, which I can compute completely without hard thinking
> (assuming the common case that B is a power of two)...

Sorry, there was an A missing. The function "sqrt_nm1" is intended to
compute a close approximation of sqrt(B^{n-2} A), with n input limbs A
and n-1 output limbs X. (Suggestions for better names are appreciated).

And B is the usual limb base.


Niels Möller. PGP-encrypted email is preferred. Keyid C0B98E26.
Internet email is subject to wholesale government surveillance.

More information about the gmp-devel mailing list