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

Timestamp:

2006-02-13T20:09:55Z (20 years ago)

Author:

Josef Cejka <malyzelenyhnus@…>

Branches:

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

Children:

f1f95f2

Parents:

e979fea

Message:

Fixed most problems with 64bit arithmetic (but division is still buggy).

File:

• softfloat/generic/common.c (modified) (3 diffs)

softfloat/generic/common.c

@@ -48 +48 @@
         };
         
-        cfrac >>= 1;
+        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) {
+                        result.parts.fraction = 0;
+                        return result;
+                }
+                //cfrac >>= 1;
+                while (cexp < 0) {
+                        cexp++;
+                        cfrac >>= 1;
+                }
+                        
+                result.parts.fraction = ((cfrac >>(64 - FLOAT64_FRACTION_SIZE - 2) ) & (~FLOAT64_HIDDEN_BIT_MASK)); 
+
+                return result;
+        }
+        
         ++cexp;
-        cfrac += (0x1 << (64 - FLOAT64_FRACTION_SIZE - 3)); 
 
         if (cfrac & (FLOAT64_HIDDEN_BIT_MASK << (64 - FLOAT64_FRACTION_SIZE - 1 ))) {
                 ++cexp;
                 cfrac >>= 1;
-                }       
+        }       
 
         /* check overflow */
@@ -65 +83 @@
         }
 
-        if (cexp < 0) {
-                /* FIXME: underflow */
-                result.parts.exp = 0;
-                if ((cexp + FLOAT64_FRACTION_SIZE) < 0) {
-                        result.parts.fraction = 0;
-                        return result;
-                }
-                cfrac >>= 1;
-                while (cexp < 0) {
-                        cexp ++;
-                        cfrac >>= 1;
-                }
-                return result;
-                
-        } else {
-                cexp ++; /*normalized*/
-                result.parts.exp = (__u32)cexp;
-        }
+        result.parts.exp = (__u32)cexp;
         
         result.parts.fraction = ((cfrac >>(64 - FLOAT64_FRACTION_SIZE - 2 ) ) & (~FLOAT64_HIDDEN_BIT_MASK)); 
@@ -88 +89 @@
         return result;  
 }
+