Anomaly in mpn_sqrtrem and mpn_rottrem

Torbjörn Granlund tg at gmplib.org
Thu Jun 11 08:12:27 UTC 2015


nisse at lysator.liu.se (Niels Möller) writes:

  Here's a sketch with some more details. I've tried to work out both
  sqrt(B^{n-1} A) and sqrt(B^{n-2}). To my surprise, they seem
  independent, not mutually recursive.
  
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)...

-- 
Torbjörn
Please encrypt, key id 0xC8601622


More information about the gmp-devel mailing list