[src/rem1.c] Fixed mpfr_rem1 in reduced exponent range (the function

is actually simplified).

git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@7490 280ebfd0-de03-0410-8827-d642c229c3f4

1 files changed, 1 insertions, 3 deletions

diff --git a/src/rem1.c b/src/rem1.c
index 9964a02f1..3c931e9dc 100644
--- a/src/rem1.c
+++ b/src/rem1.c
@@ -196,9 +196,7 @@ mpfr_rem1 (mpfr_ptr rem, long *quo, mpfr_rnd_t rnd_q,
       /* 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;
+      inex = mpfr_set_z_2exp (rem, r, ex > ey ? ey : ex, rnd);
     }
 
   if (quo)