diff -r 000000000000 -r 39b1c755f19b patches/mpfr/2.4.1/100-remainder-neg.patch
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/patches/mpfr/2.4.1/100-remainder-neg.patch	Sun Apr 26 10:47:00 2009 +0000
@@ -0,0 +1,123 @@
+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";
+ }