fixed possible overflow when EXP=1024

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

1 files changed, 5 insertions, 4 deletions

diff --git a/get_ld.c b/get_ld.c
index 9ca8b81ab..749043572 100644
--- a/get_ld.c
+++ b/get_ld.c
@@ -20,6 +20,7 @@ along with the MPFR Library; see the file COPYING.LIB.  If not, write to
 the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
 MA 02111-1307, USA. */
 
+#include <stdio.h>
 #include <float.h>
 
 #include "gmp.h"
@@ -64,10 +65,10 @@ mpfr_get_ld (mpfr_srcptr x, mp_rnd_t rnd_mode)
       mpfr_set (y, x, rnd_mode);
       negative = MPFR_SIGN(y) < 0;
       e = MPFR_EXP(y);
-      if (e > 1024)
+      if (e > 1023)
         {
-          sh = e - 1024;
-          MPFR_EXP(y) = 1024;
+          sh = e - 1023;
+          MPFR_EXP(y) = 1023;
         }
       else if (e < -1021)
         {
@@ -78,7 +79,7 @@ mpfr_get_ld (mpfr_srcptr x, mp_rnd_t rnd_mode)
         {
           sh = 0;
         }
-      /* now -1021 <= e - sh = EXP(y) <= 1024 */
+      /* now -1021 <= e - sh = EXP(y) <= 1023 */
       r = 0.0;
       mpfr_init2 (z, DBL_MANT_DIG);