Torbjorn Granlund <t...@gmplib.org> writes:
> Do you have a proposal for two separate functions?
I'll think about it...
> Perhaps broot.c should use Niels' function instead, I cannot recall
> the differences, if any.
My bsqrt uses an iteration converging to a^{-1/2}, and broot uses an
iteration converging to a^{1/n - 1}. Both division free.
They should use mulmod_bnm1 rather than mullo for larger sizes, but they
currently don't. I haven't done any serious benchmarking.
/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
