diff options
Diffstat (limited to 'src/sqrt.c')
-rw-r--r-- | src/sqrt.c | 11 |
1 files changed, 1 insertions, 10 deletions
diff --git a/src/sqrt.c b/src/sqrt.c index 3c56264ba..57ef54a35 100644 --- a/src/sqrt.c +++ b/src/sqrt.c @@ -133,16 +133,7 @@ mpfr_sqrt (mpfr_ptr r, mpfr_srcptr u, mpfr_rnd_t rnd_mode) /* sticky0 is non-zero iff the truncated part of the input is non-zero */ - /* for size above about 500 limbs, mpn_rootrem with NULL 2nd argument is - faster than mpn_sqrtrem, thus use it if available and if the user asked - to use GMP internal functions - (https://gmplib.org/list-archives/gmp-devel/2010-September/001654.html) */ -#if defined(WANT_GMP_INTERNALS) && defined(HAVE___GMPN_ROOTREM) - if (rrsize >= 500) - tsize = __gmpn_rootrem (rp, NULL, sp, rrsize, 2); - else -#endif - tsize = mpn_sqrtrem (rp, sp, sp, rrsize); + tsize = mpn_sqrtrem (rp, NULL, sp, rrsize); /* a return value of zero in mpn_sqrtrem indicates a perfect square */ sticky = sticky0 || tsize != 0; |