Special case biased exponents < 1; Get boundary for Inf right

1 files changed, 13 insertions, 1 deletions

diff --git a/insert-dbl.c b/insert-dbl.c
index c63eaa557..02e8a6b65 100644
--- a/insert-dbl.c
+++ b/insert-dbl.c
@@ -45,11 +45,23 @@ __gmp_scale2 (d, exp)
     x.d = d;
     exp += x.s.exp;
     x.s.exp = exp;
-    if (exp >= 2048)
+    if (exp >= 2047)
       {
+	/* Return +-infinity */
 	x.s.exp = 2047;
 	x.s.manl = x.s.manh = 0;
       }
+    else if (exp < 1)
+      {
+	x.s.exp = 1;		/* smallest exponent (biased) */
+	/* Divide result by 2 until we have scaled it to the right IEEE
+	   denormalized number, but stop if it becomes zero.  */
+	while (exp < 1 && x.d != 0)
+	  {
+	    x.d *= 0.5;
+	    exp++;
+	  }
+      }
     return x.d;
   }
 #else