fast inversion

2015-02-03 Thread paul zimmermann
Hi, I did update my notes and code about fast inversion (which were used as basis for Algorithm 3.5 ApproximateReciprocal in "Modern Computer Arithmetic"): http://www.loria.fr/~zimmerma/papers/invert.{c,pdf} The main change is that now the result is uniquely defined (instead of 2 ulps mar

Re: Adding support for R6 of MIPS architecture

2015-02-03 Thread Torbjörn Granlund
Marc Glisse writes: Apparently not, the motivation for the patch is that multu has disappeared... Then I see no other robust approach than making mpn/mipsnomultu_64 (or somesuch). Well, an analogous robust approach would be moving the dmultu code into mpn/mips64/dmultu and put the r6 code i