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
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