Torbjorn Granlund <t...@gmplib.org> writes: > Is your code different from perfpow.c's binv_sqroot?
I think so. I don't fully understand the binv_sqroot, but some differences: * As far as I see, it doesn't take any advantage of cancellation in the iteration. * It uses powlo in some way which I haven't figured out. * And they return different values, a^{1/2} and a^{-1/2}, respectively (but that's a minor issue; I'm not sure which value makes the most sense to return). * And the binv_sqrt function is a lot simpler than mine. 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 http://gmplib.org/mailman/listinfo/gmp-devel