diff options
Diffstat (limited to 'Build/source/libs/mpfr/mpfr-src/src/mpn_exp.c')
-rw-r--r-- | Build/source/libs/mpfr/mpfr-src/src/mpn_exp.c | 4 |
1 files changed, 3 insertions, 1 deletions
diff --git a/Build/source/libs/mpfr/mpfr-src/src/mpn_exp.c b/Build/source/libs/mpfr/mpfr-src/src/mpn_exp.c index 70417b870f1..de6e11fa088 100644 --- a/Build/source/libs/mpfr/mpfr-src/src/mpn_exp.c +++ b/Build/source/libs/mpfr/mpfr-src/src/mpn_exp.c @@ -91,7 +91,9 @@ mpfr_mpn_exp (mp_limb_t *a, mpfr_exp_t *exp_r, int b, mpfr_exp_t e, size_t n) n1 = (n * GMP_NUMB_BITS - bits) / GMP_NUMB_BITS; /* square of A : {c+2n1, 2(n-n1)} = {a+n1, n-n1}^2 */ - mpn_sqr_n (c + 2 * n1, a + n1, n - n1); + /* TODO: we should use a short square here, but this needs to redo + the error analysis */ + mpn_sqr (c + 2 * n1, a + n1, n - n1); /* set {c+n, 2n1-n} to 0 : {c, n} = {a, n}^2*K^n */ |