Changeset f1f95f2 in mainline for softfloat/generic/common.c

Timestamp:

2006-02-13T23:06:04Z (19 years ago)

Author:

Josef Cejka <malyzelenyhnus@…>

Branches:

lfn, master, serial, ticket/834-toolchain-update, topic/msim-upgrade, topic/simplify-dev-export

Children:

d9f51ccc

Parents:

d3ca210

Message:

Fixed bugs in 64bit float division.

File:

• softfloat/generic/common.c (modified) (1 diff)

softfloat/generic/common.c

@@ -48 +48 @@
         };
         
-        cfrac += (0x1 << (64 - FLOAT64_FRACTION_SIZE - 3)); 
-        
         if ((cexp < 0) || ( cexp == 0 && (!(cfrac & (FLOAT64_HIDDEN_BIT_MASK << (64 - FLOAT64_FRACTION_SIZE - 1)))))) {
                 /* FIXME: underflow */
                 result.parts.exp = 0;
-                if ((cexp + FLOAT64_FRACTION_SIZE) < 0) {
+                if ((cexp + FLOAT64_FRACTION_SIZE + 1) < 0) { /* +1 is place for rounding */
                         result.parts.fraction = 0;
                         return result;
                 }
-                //cfrac >>= 1;
+                
                 while (cexp < 0) {
                         cexp++;
                         cfrac >>= 1;
                 }
+        
+                cfrac += (0x1 << (64 - FLOAT64_FRACTION_SIZE - 3)); 
+                
+                if (!(cfrac & (FLOAT64_HIDDEN_BIT_MASK << (64 - FLOAT64_HIDDEN_BIT_MASK - 1)))) {
                         
-                result.parts.fraction = ((cfrac >>(64 - FLOAT64_FRACTION_SIZE - 2) ) & (~FLOAT64_HIDDEN_BIT_MASK)); 
-
-                return result;
+                        result.parts.fraction = ((cfrac >>(64 - FLOAT64_FRACTION_SIZE - 2) ) & (~FLOAT64_HIDDEN_BIT_MASK)); 
+                        return result;
+                }       
+        } else {
+                cfrac += (0x1 << (64 - FLOAT64_FRACTION_SIZE - 3)); 
         }