|
 |
d54aba |
diff -Naurd mpfr-3.1.6-a/PATCHES mpfr-3.1.6-b/PATCHES
|
|
|
d54aba |
--- mpfr-3.1.6-a/PATCHES 2017-12-15 12:37:44.074256548 +0000
|
|
|
d54aba |
+++ mpfr-3.1.6-b/PATCHES 2017-12-15 12:37:44.146256304 +0000
|
|
|
d54aba |
@@ -0,0 +1 @@
|
|
|
d54aba |
+root
|
|
|
d54aba |
diff -Naurd mpfr-3.1.6-a/VERSION mpfr-3.1.6-b/VERSION
|
|
|
d54aba |
--- mpfr-3.1.6-a/VERSION 2017-10-26 13:55:51.236013121 +0000
|
|
|
d54aba |
+++ mpfr-3.1.6-b/VERSION 2017-12-15 12:37:44.142256319 +0000
|
|
|
d54aba |
@@ -1 +1 @@
|
|
|
d54aba |
-3.1.6-p1
|
|
|
d54aba |
+3.1.6-p2
|
|
|
d54aba |
diff -Naurd mpfr-3.1.6-a/src/mpfr.h mpfr-3.1.6-b/src/mpfr.h
|
|
|
d54aba |
--- mpfr-3.1.6-a/src/mpfr.h 2017-10-26 13:55:51.232013138 +0000
|
|
|
d54aba |
+++ mpfr-3.1.6-b/src/mpfr.h 2017-12-15 12:37:44.142256319 +0000
|
|
|
d54aba |
@@ -27,7 +27,7 @@
|
|
|
d54aba |
#define MPFR_VERSION_MAJOR 3
|
|
|
d54aba |
#define MPFR_VERSION_MINOR 1
|
|
|
d54aba |
#define MPFR_VERSION_PATCHLEVEL 6
|
|
|
d54aba |
-#define MPFR_VERSION_STRING "3.1.6-p1"
|
|
|
d54aba |
+#define MPFR_VERSION_STRING "3.1.6-p2"
|
|
|
d54aba |
|
|
|
d54aba |
/* Macros dealing with MPFR VERSION */
|
|
|
d54aba |
#define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c))
|
|
|
d54aba |
diff -Naurd mpfr-3.1.6-a/src/root.c mpfr-3.1.6-b/src/root.c
|
|
|
d54aba |
--- mpfr-3.1.6-a/src/root.c 2017-01-01 01:39:09.000000000 +0000
|
|
|
d54aba |
+++ mpfr-3.1.6-b/src/root.c 2017-12-15 12:37:44.118256398 +0000
|
|
|
d54aba |
@@ -107,6 +107,11 @@
|
|
|
d54aba |
MPFR_RET_NAN;
|
|
|
d54aba |
}
|
|
|
d54aba |
|
|
|
d54aba |
+ /* Special case |x| = 1. Note that if x = -1, then k is odd
|
|
|
d54aba |
+ (NaN results have already been filtered), so that y = -1. */
|
|
|
d54aba |
+ if (mpfr_cmpabs (x, __gmpfr_one) == 0)
|
|
|
d54aba |
+ return mpfr_set (y, x, rnd_mode);
|
|
|
d54aba |
+
|
|
|
d54aba |
/* General case */
|
|
|
d54aba |
|
|
|
d54aba |
/* For large k, use exp(log(x)/k). The threshold of 100 seems to be quite
|
|
|
d54aba |
@@ -188,6 +193,14 @@
|
|
|
d54aba |
Assume all special cases have been eliminated before.
|
|
|
d54aba |
In the extended exponent range, overflows/underflows are not possible.
|
|
|
d54aba |
Assume x > 0, or x < 0 and k odd.
|
|
|
d54aba |
+ Also assume |x| <> 1 because log(1) = 0, which does not have an exponent
|
|
|
d54aba |
+ and would yield a failure in the error bound computation. A priori, this
|
|
|
d54aba |
+ constraint is quite artificial because if |x| is close enough to 1, then
|
|
|
d54aba |
+ the exponent of log|x| does not need to be used (in the code, err would
|
|
|
d54aba |
+ be 1 in such a domain). So this constraint |x| <> 1 could be avoided in
|
|
|
d54aba |
+ the code. However, this is an exact case easy to detect, so that such a
|
|
|
d54aba |
+ change would be useless. Values very close to 1 are not an issue, since
|
|
|
d54aba |
+ an underflow is not possible before the MPFR_GET_EXP.
|
|
|
d54aba |
*/
|
|
|
d54aba |
static int
|
|
|
d54aba |
mpfr_root_aux (mpfr_ptr y, mpfr_srcptr x, unsigned long k, mpfr_rnd_t rnd_mode)
|
|
|
d54aba |
@@ -219,7 +232,8 @@
|
|
|
d54aba |
mpfr_log (t, absx, MPFR_RNDN);
|
|
|
d54aba |
/* t = log|x| * (1 + theta) with |theta| <= 2^(-w) */
|
|
|
d54aba |
mpfr_div_ui (t, t, k, MPFR_RNDN);
|
|
|
d54aba |
- expt = MPFR_GET_EXP (t);
|
|
|
d54aba |
+ /* No possible underflow in mpfr_log and mpfr_div_ui. */
|
|
|
d54aba |
+ expt = MPFR_GET_EXP (t); /* assumes t <> 0 */
|
|
|
d54aba |
/* t = log|x|/k * (1 + theta) + eps with |theta| <= 2^(-w)
|
|
|
d54aba |
and |eps| <= 1/2 ulp(t), thus the total error is bounded
|
|
|
d54aba |
by 1.5 * 2^(expt - w) */
|
|
|
d54aba |
diff -Naurd mpfr-3.1.6-a/src/version.c mpfr-3.1.6-b/src/version.c
|
|
|
d54aba |
--- mpfr-3.1.6-a/src/version.c 2017-10-26 13:55:51.232013138 +0000
|
|
|
d54aba |
+++ mpfr-3.1.6-b/src/version.c 2017-12-15 12:37:44.142256319 +0000
|
|
|
d54aba |
@@ -25,5 +25,5 @@
|
|
|
d54aba |
const char *
|
|
|
d54aba |
mpfr_get_version (void)
|
|
|
d54aba |
{
|
|
|
d54aba |
- return "3.1.6-p1";
|
|
|
d54aba |
+ return "3.1.6-p2";
|
|
|
d54aba |
}
|
|
|
d54aba |
diff -Naurd mpfr-3.1.6-a/tests/troot.c mpfr-3.1.6-b/tests/troot.c
|
|
|
d54aba |
--- mpfr-3.1.6-a/tests/troot.c 2017-01-01 01:39:09.000000000 +0000
|
|
|
d54aba |
+++ mpfr-3.1.6-b/tests/troot.c 2017-12-15 12:37:44.118256398 +0000
|
|
|
d54aba |
@@ -405,6 +405,26 @@
|
|
|
d54aba |
mpfr_clears (x, y1, y2, (mpfr_ptr) 0);
|
|
|
d54aba |
}
|
|
|
d54aba |
|
|
|
d54aba |
+static void
|
|
|
d54aba |
+bug20171214 (void)
|
|
|
d54aba |
+{
|
|
|
d54aba |
+ mpfr_t x, y;
|
|
|
d54aba |
+ int inex;
|
|
|
d54aba |
+
|
|
|
d54aba |
+ mpfr_init2 (x, 805);
|
|
|
d54aba |
+ mpfr_init2 (y, 837);
|
|
|
d54aba |
+ mpfr_set_ui (x, 1, MPFR_RNDN);
|
|
|
d54aba |
+ inex = mpfr_root (y, x, 120, MPFR_RNDN);
|
|
|
d54aba |
+ MPFR_ASSERTN (inex == 0);
|
|
|
d54aba |
+ MPFR_ASSERTN (mpfr_cmp_ui (y, 1) == 0);
|
|
|
d54aba |
+ mpfr_set_si (x, -1, MPFR_RNDN);
|
|
|
d54aba |
+ inex = mpfr_root (y, x, 121, MPFR_RNDN);
|
|
|
d54aba |
+ MPFR_ASSERTN (inex == 0);
|
|
|
d54aba |
+ MPFR_ASSERTN (mpfr_cmp_si (y, -1) == 0);
|
|
|
d54aba |
+ mpfr_clear (x);
|
|
|
d54aba |
+ mpfr_clear (y);
|
|
|
d54aba |
+}
|
|
|
d54aba |
+
|
|
|
d54aba |
int
|
|
|
d54aba |
main (void)
|
|
|
d54aba |
{
|
|
|
d54aba |
@@ -415,6 +435,7 @@
|
|
|
d54aba |
|
|
|
d54aba |
tests_start_mpfr ();
|
|
|
d54aba |
|
|
|
d54aba |
+ bug20171214 ();
|
|
|
d54aba |
exact_powers (3, 1000);
|
|
|
d54aba |
special ();
|
|
|
d54aba |
bigint ();
|