diff options
author | Marco Bodrato <bodrato@mail.dm.unipi.it> | 2020-10-29 07:42:04 +0100 |
---|---|---|
committer | Marco Bodrato <bodrato@mail.dm.unipi.it> | 2020-10-29 07:42:04 +0100 |
commit | d0397f95bfcdfa2c3952b06295d69c849fc73558 (patch) | |
tree | ae0dc3c077bc81d95cf5794dbdf33da1a8cc61c0 /mpz/stronglucas.c | |
parent | 4dc90a38398bcd52993515631e6ffcca9b4f5697 (diff) | |
download | gmp-d0397f95bfcdfa2c3952b06295d69c849fc73558.tar.gz |
limb_apprsqrt: Slightly faster formula
Diffstat (limited to 'mpz/stronglucas.c')
-rw-r--r-- | mpz/stronglucas.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/mpz/stronglucas.c b/mpz/stronglucas.c index e9a71370e..cf9e7ff40 100644 --- a/mpz/stronglucas.c +++ b/mpz/stronglucas.c @@ -5,7 +5,7 @@ CERTAIN TO BE SUBJECT TO INCOMPATIBLE CHANGES OR DISAPPEAR COMPLETELY IN FUTURE GNU MP RELEASES. -Copyright 2018 Free Software Foundation, Inc. +Copyright 2018, 2020 Free Software Foundation, Inc. Contributed by Marco Bodrato. @@ -52,7 +52,7 @@ limb_apprsqrt (mp_limb_t x) ASSERT (x > 2); count_leading_zeros (s, x); s = (GMP_LIMB_BITS - s) >> 1; - return ((CNST_LIMB(1) << s) + (x >> s)) >> 1; + return ((CNST_LIMB(1) << (s - 1)) + (x >> 1 >> s)); } static int |