Open fredrik-johansson opened 1 month ago
I will just state it here as well: For inverses, I believe for small $n$ the fastest method is via Newton iteration (and based off of GMP, I suppose this extends to all numbers). Hardcoded routines for inverses could be implemented for limb counts that are powers of two, just generalizing mpn_invert_limb
.
There is also the basecase algorithm used by mpfr_divhigh_n_basecase
and the variant described here: https://inria.hal.science/hal-04557431v1/document
There is also the basecase algorithm used by
mpfr_divhigh_n_basecase
and the variant described here: https://inria.hal.science/hal-04557431v1/document
I saw that one. I'm wondering how a fast mpn_invert
joined with Granlund-Möller 2n-by-n division algorithm would compare to Sukop's and Zimmermann's new algorithm.
Before I forget anything