[Top][All Lists]
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [INFO] Yet Another sqrt algorithm.
From: |
Torbjorn Granlund |
Subject: |
Re: [INFO] Yet Another sqrt algorithm. |
Date: |
12 Dec 2001 16:11:20 +0100 |
User-agent: |
Gnus/5.0807 (Gnus v5.8.7) Emacs/20.7 |
address@hidden (Paul Zimmermann) writes:
The current release is GMP 4.0. Please download
<ftp://ftp.gnu.org/gnu/gmp/gmp-4.0.tar.gz>. The algorithm
used currently was developed by Paul Zimmermann. I don't
know if you code could improve some part of that algorithm,
since my understanding of the algorithm is limited.
Paul, what do you think?
It could improve the one-limb code, which currently uses a 8-bit
table approximation, together with 2 or 3 Newton's loops, depending
on the word size. See function mpn_sqrtrem1 in mpn/sqrtrem.c from GMP 4.0.
Want to work on improving the gmp/mpn/generic/sqrtrem.c
basecase code? The proposed code could perhaps be faster
for up to a small number of limbs. The current basecase is
just for a single limb.
--
Torbjörn