new approximate inverse code for GMP
Paul.Zimmermann at loria.fr
Thu Aug 31 16:14:38 CEST 2006
I've revised my note on fast approximate inverse for GMP.
The code is available too: the mpn_invert2 routine in file
invert.c could replace the mpn_inv routine from
except it returns a result which may be 1 bit too small.
It would be great if someone could check the proof of Lemma 2.
PS: don't look at sections 3 (approximate quotient) and 4 (quotient with
remainder). They are obsolete with the "New truncating division interface and
code" on <http://www.swox.se/gmp/devel/>.
More information about the gmp-devel