t...@gmplib.org (Torbjörn Granlund) writes:
count_leading_zeros (ncnt, n0)
count_leading_zeros (dcnt, d0)
ni = ... extract high k bits from n0 ...
di = ... extract high \ell bits from d0 ...
q = qtab[ni<<\ell+di];
... adjust ...
Well, perhaps
q = qtab[ni<<\ell+di] << (ncnt - dcnt)
is closer.
if (n0 >= (d0 << 3))
q = n0 / d0;
else {
q = 0;
if (n0 > (d0 << 2))
q += 4
n0 -= d0 << 2
if (n0 > (d0 << 1))
q += 2
n0 -= d0 << 1
if (n0 > (d0 << 0))
q += 1
n0 -= d0
}
All '>' there should be '>='.
--
Torbjörn
Please encrypt, key id 0xC8601622
_______________________________________________
gmp-devel mailing list
gmp-devel@gmplib.org
https://gmplib.org/mailman/listinfo/gmp-devel
