t...@gmplib.org (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. Regards, /Niels -- Niels Möller. PGP-encrypted email is preferred. Keyid C0B98E26. Internet email is subject to wholesale government surveillance. _______________________________________________ gmp-devel mailing list gmp-devel@gmplib.org https://gmplib.org/mailman/listinfo/gmp-devel