paul zimmermann Paul.Zimmermann at
Wed Jul 20 13:23:24 UTC 2016


the attached file contains a new version (mpn_sqrtrem2_new) of mpn_sqrtrem2
for 64-bit computers. On my computer i5-4590 it runs in about 65 cycles
whereas the gmp-6.1.1 version runs in about 96 cycles.

It uses an auxiliary routine (mpn_rsqrtrem1) which computes a 41-bit
approximation of 1/sqrt(a).

Before I try to adapt those two routines to 32-bit limbs, I'd like someone
to review my code (comments are included).

-------------- next part --------------
A non-text attachment was scrubbed...
Name: sqrtrem.c
Type: text/x-csrc
Size: 27379 bytes
Desc: not available
URL: <>

More information about the gmp-devel mailing list