diff -Naurd mpfr-2.4.1-a/PATCHES mpfr-2.4.1-b/PATCHES --- mpfr-2.4.1-a/PATCHES 2009-02-20 09:43:17.000000000 +0000 +++ mpfr-2.4.1-b/PATCHES 2009-02-27 16:56:29.000000000 +0000 @@ -0,0 +1 @@ +remainder-neg diff -Naurd mpfr-2.4.1-a/VERSION mpfr-2.4.1-b/VERSION --- mpfr-2.4.1-a/VERSION 2009-02-25 16:16:08.000000000 +0000 +++ mpfr-2.4.1-b/VERSION 2009-02-27 16:55:37.000000000 +0000 @@ -1 +1 @@ -2.4.1 +2.4.1-p1 diff -Naurd mpfr-2.4.1-a/mpfr.h mpfr-2.4.1-b/mpfr.h --- mpfr-2.4.1-a/mpfr.h 2009-02-25 16:16:08.000000000 +0000 +++ mpfr-2.4.1-b/mpfr.h 2009-02-27 16:55:38.000000000 +0000 @@ -27,7 +27,7 @@ #define MPFR_VERSION_MAJOR 2 #define MPFR_VERSION_MINOR 4 #define MPFR_VERSION_PATCHLEVEL 1 -#define MPFR_VERSION_STRING "2.4.1" +#define MPFR_VERSION_STRING "2.4.1-p1" /* Macros dealing with MPFR VERSION */ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c)) diff -Naurd mpfr-2.4.1-a/rem1.c mpfr-2.4.1-b/rem1.c --- mpfr-2.4.1-a/rem1.c 2009-02-20 09:43:17.000000000 +0000 +++ mpfr-2.4.1-b/rem1.c 2009-02-27 16:55:15.000000000 +0000 @@ -170,7 +170,12 @@ } if (mpz_cmp_ui (r, 0) == 0) - inex = mpfr_set_ui (rem, 0, GMP_RNDN); + { + inex = mpfr_set_ui (rem, 0, GMP_RNDN); + /* take into account sign of x */ + if (signx < 0) + mpfr_neg (rem, rem, GMP_RNDN); + } else { if (rnd_q == GMP_RNDN) @@ -190,6 +195,9 @@ *quo += 1; } } + /* take into account sign of x */ + if (signx < 0) + mpz_neg (r, r); inex = mpfr_set_z (rem, r, rnd); /* if ex > ey, rem should be multiplied by 2^ey, else by 2^ex */ MPFR_EXP (rem) += (ex > ey) ? ey : ex; @@ -198,13 +206,6 @@ if (quo) *quo *= sign; - /* take into account sign of x */ - if (signx < 0) - { - mpfr_neg (rem, rem, GMP_RNDN); - inex = -inex; - } - mpz_clear (mx); mpz_clear (my); mpz_clear (r); diff -Naurd mpfr-2.4.1-a/tests/tremquo.c mpfr-2.4.1-b/tests/tremquo.c --- mpfr-2.4.1-a/tests/tremquo.c 2009-02-20 09:43:15.000000000 +0000 +++ mpfr-2.4.1-b/tests/tremquo.c 2009-02-27 16:55:15.000000000 +0000 @@ -25,6 +25,36 @@ #include "mpfr-test.h" +static void +bug20090227 (void) +{ + mpfr_t x, y, r1, r2; + int inex1, inex2; + + mpfr_init2 (x, 118); + mpfr_init2 (y, 181); + mpfr_init2 (r1, 140); + mpfr_init2 (r2, 140); + mpfr_set_si (x, -1, GMP_RNDN); + mpfr_set_str_binary (y, "1.100100100001111110110101010001000100001011010001100001000110100110001001100011001100010100010111000000011011100000111001101000100101001000000100100111000001000100010100110011111010"); + inex1 = mpfr_remainder (r1, x, y, GMP_RNDU); + /* since the quotient is -1, r1 is the rounding of x+y */ + inex2 = mpfr_add (r2, x, y, GMP_RNDU); + if (mpfr_cmp (r1, r2)) + { + printf ("Error in mpfr_remainder (bug20090227)\n"); + printf ("Expected "); + mpfr_dump (r2); + printf ("Got "); + mpfr_dump (r1); + exit (1); + } + mpfr_clear (x); + mpfr_clear (y); + mpfr_clear (r1); + mpfr_clear (r2); +} + int main (int argc, char *argv[]) { @@ -50,6 +80,8 @@ tests_start_mpfr (); + bug20090227 (); + mpfr_init (x); mpfr_init (y); mpfr_init (r); diff -Naurd mpfr-2.4.1-a/version.c mpfr-2.4.1-b/version.c --- mpfr-2.4.1-a/version.c 2009-02-25 16:16:08.000000000 +0000 +++ mpfr-2.4.1-b/version.c 2009-02-27 16:55:38.000000000 +0000 @@ -25,5 +25,5 @@ const char * mpfr_get_version (void) { - return "2.4.1"; + return "2.4.1-p1"; } diff -Naurd mpfr-2.4.1-a/PATCHES mpfr-2.4.1-b/PATCHES --- mpfr-2.4.1-a/PATCHES 2009-03-04 13:15:05.000000000 +0000 +++ mpfr-2.4.1-b/PATCHES 2009-03-04 13:17:04.000000000 +0000 @@ -0,0 +1 @@ +assert diff -Naurd mpfr-2.4.1-a/VERSION mpfr-2.4.1-b/VERSION --- mpfr-2.4.1-a/VERSION 2009-02-27 16:55:37.000000000 +0000 +++ mpfr-2.4.1-b/VERSION 2009-03-04 13:16:46.000000000 +0000 @@ -1 +1 @@ -2.4.1-p1 +2.4.1-p2 diff -Naurd mpfr-2.4.1-a/mpfr-gmp.c mpfr-2.4.1-b/mpfr-gmp.c --- mpfr-2.4.1-a/mpfr-gmp.c 2009-02-20 09:43:17.000000000 +0000 +++ mpfr-2.4.1-b/mpfr-gmp.c 2009-03-04 13:16:39.000000000 +0000 @@ -301,7 +301,7 @@ if (linenum != -1) fprintf (stderr, "%d: ", linenum); } - fprintf (stderr, " assertion failed: %s\n", expr); + fprintf (stderr, "MPFR assertion failed: %s\n", expr); abort(); } diff -Naurd mpfr-2.4.1-a/mpfr.h mpfr-2.4.1-b/mpfr.h --- mpfr-2.4.1-a/mpfr.h 2009-02-27 16:55:38.000000000 +0000 +++ mpfr-2.4.1-b/mpfr.h 2009-03-04 13:16:46.000000000 +0000 @@ -27,7 +27,7 @@ #define MPFR_VERSION_MAJOR 2 #define MPFR_VERSION_MINOR 4 #define MPFR_VERSION_PATCHLEVEL 1 -#define MPFR_VERSION_STRING "2.4.1-p1" +#define MPFR_VERSION_STRING "2.4.1-p2" /* Macros dealing with MPFR VERSION */ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c)) diff -Naurd mpfr-2.4.1-a/version.c mpfr-2.4.1-b/version.c --- mpfr-2.4.1-a/version.c 2009-02-27 16:55:38.000000000 +0000 +++ mpfr-2.4.1-b/version.c 2009-03-04 13:16:46.000000000 +0000 @@ -25,5 +25,5 @@ const char * mpfr_get_version (void) { - return "2.4.1-p1"; + return "2.4.1-p2"; } diff -Naurd mpfr-2.4.1-a/tests/tfprintf.c mpfr-2.4.1-b/tests/tfprintf.c --- mpfr-2.4.1-a/tests/tfprintf.c 2009-02-20 09:43:15.000000000 +0000 +++ mpfr-2.4.1-b/tests/tfprintf.c 2009-03-10 01:12:45.000000000 +0000 @@ -195,7 +195,7 @@ check_vfprintf (fout, "a. %Pu, b. %c, c. %Zi%Zn", prec, ch, mpz, &mpz); check_length_with_cmp (6, mpz, 17, mpz_cmp_ui (mpz, 17), Zi); check_vfprintf (fout, "%% a. %#.0RNg, b. %Qx%Rn, c. %p", mpfr, mpq, &mpfr, - &i); + (void *) &i); check_length_with_cmp (7, mpfr, 16, mpfr_cmp_ui (mpfr, 16), Rg); #ifndef NPRINTF_T diff -Naurd mpfr-2.4.1-a/tests/tprintf.c mpfr-2.4.1-b/tests/tprintf.c --- mpfr-2.4.1-a/tests/tprintf.c 2009-02-20 09:43:15.000000000 +0000 +++ mpfr-2.4.1-b/tests/tprintf.c 2009-03-10 01:12:45.000000000 +0000 @@ -191,7 +191,8 @@ check_length (5, sz, 34, zu); check_vprintf ("a. %Pu, b. %c, c. %RUG, d. %Zi%Zn", prec, ch, mpfr, mpz, &mpz); check_length_with_cmp (6, mpz, 24, mpz_cmp_ui (mpz, 24), Zi); - check_vprintf ("%% a. %#.0RNg, b. %Qx%Rn c. %p", mpfr, mpq, &mpfr, &i); + check_vprintf ("%% a. %#.0RNg, b. %Qx%Rn c. %p", + mpfr, mpq, &mpfr, (void *) &i); check_length_with_cmp (7, mpfr, 16, mpfr_cmp_ui (mpfr, 16), Rg); #ifndef NPRINTF_T diff -Naurd mpfr-2.4.1-a/PATCHES mpfr-2.4.1-b/PATCHES --- mpfr-2.4.1-a/PATCHES 2009-03-10 01:44:39.000000000 +0000 +++ mpfr-2.4.1-b/PATCHES 2009-03-10 01:44:40.000000000 +0000 @@ -0,0 +1 @@ +vasprintf-mp_limb_t diff -Naurd mpfr-2.4.1-a/VERSION mpfr-2.4.1-b/VERSION --- mpfr-2.4.1-a/VERSION 2009-03-04 13:16:46.000000000 +0000 +++ mpfr-2.4.1-b/VERSION 2009-03-10 01:44:40.000000000 +0000 @@ -1 +1 @@ -2.4.1-p2 +2.4.1-p4 diff -Naurd mpfr-2.4.1-a/mpfr.h mpfr-2.4.1-b/mpfr.h --- mpfr-2.4.1-a/mpfr.h 2009-03-04 13:16:46.000000000 +0000 +++ mpfr-2.4.1-b/mpfr.h 2009-03-10 01:44:40.000000000 +0000 @@ -27,7 +27,7 @@ #define MPFR_VERSION_MAJOR 2 #define MPFR_VERSION_MINOR 4 #define MPFR_VERSION_PATCHLEVEL 1 -#define MPFR_VERSION_STRING "2.4.1-p2" +#define MPFR_VERSION_STRING "2.4.1-p4" /* Macros dealing with MPFR VERSION */ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c)) diff -Naurd mpfr-2.4.1-a/vasprintf.c mpfr-2.4.1-b/vasprintf.c --- mpfr-2.4.1-a/vasprintf.c 2009-02-20 09:43:17.000000000 +0000 +++ mpfr-2.4.1-b/vasprintf.c 2009-03-10 01:44:40.000000000 +0000 @@ -398,7 +398,7 @@ (void) va_arg ((ap), mpq_srcptr); \ break; \ case MP_LIMB_ARG: \ - (void) va_arg ((ap), mp_ptr); \ + (void) va_arg ((ap), mp_limb_t); \ break; \ case MP_LIMB_ARRAY_ARG: \ (void) va_arg ((ap), mp_ptr); \ diff -Naurd mpfr-2.4.1-a/version.c mpfr-2.4.1-b/version.c --- mpfr-2.4.1-a/version.c 2009-03-04 13:16:46.000000000 +0000 +++ mpfr-2.4.1-b/version.c 2009-03-10 01:44:40.000000000 +0000 @@ -25,5 +25,5 @@ const char * mpfr_get_version (void) { - return "2.4.1-p2"; + return "2.4.1-p4"; } diff -Naurd mpfr-2.4.1-a/PATCHES mpfr-2.4.1-b/PATCHES --- mpfr-2.4.1-a/PATCHES 2009-03-13 02:50:47.000000000 +0000 +++ mpfr-2.4.1-b/PATCHES 2009-03-13 02:50:48.000000000 +0000 @@ -0,0 +1 @@ +zeta_ui-shift diff -Naurd mpfr-2.4.1-a/VERSION mpfr-2.4.1-b/VERSION --- mpfr-2.4.1-a/VERSION 2009-03-10 01:44:40.000000000 +0000 +++ mpfr-2.4.1-b/VERSION 2009-03-13 02:50:48.000000000 +0000 @@ -1 +1 @@ -2.4.1-p4 +2.4.1-p5 diff -Naurd mpfr-2.4.1-a/mpfr.h mpfr-2.4.1-b/mpfr.h --- mpfr-2.4.1-a/mpfr.h 2009-03-10 01:44:40.000000000 +0000 +++ mpfr-2.4.1-b/mpfr.h 2009-03-13 02:50:48.000000000 +0000 @@ -27,7 +27,7 @@ #define MPFR_VERSION_MAJOR 2 #define MPFR_VERSION_MINOR 4 #define MPFR_VERSION_PATCHLEVEL 1 -#define MPFR_VERSION_STRING "2.4.1-p4" +#define MPFR_VERSION_STRING "2.4.1-p5" /* Macros dealing with MPFR VERSION */ #define MPFR_VERSION_NUM(a,b,c) (((a) << 16L) | ((b) << 8) | (c)) diff -Naurd mpfr-2.4.1-a/version.c mpfr-2.4.1-b/version.c --- mpfr-2.4.1-a/version.c 2009-03-10 01:44:40.000000000 +0000 +++ mpfr-2.4.1-b/version.c 2009-03-13 02:50:48.000000000 +0000 @@ -25,5 +25,5 @@ const char * mpfr_get_version (void) { - return "2.4.1-p4"; + return "2.4.1-p5"; } diff -Naurd mpfr-2.4.1-a/zeta_ui.c mpfr-2.4.1-b/zeta_ui.c --- mpfr-2.4.1-a/zeta_ui.c 2009-02-20 09:43:17.000000000 +0000 +++ mpfr-2.4.1-b/zeta_ui.c 2009-03-13 02:50:48.000000000 +0000 @@ -177,7 +177,9 @@ mpz_mul_ui (t, t, 2 * k - 1); } mpz_div_2exp (t, t, 1); - if (n < 1UL << (BITS_PER_MP_LIMB / 2)) + /* Warning: the test below assumes that an unsigned long + has no padding bits. */ + if (n < 1UL << ((sizeof(unsigned long) * CHAR_BIT) / 2)) /* (n - k + 1) * (n + k - 1) < n^2 */ mpz_divexact_ui (t, t, (n - k + 1) * (n + k - 1)); else