Changeset e979fea in mainline for softfloat/generic/div.c

Timestamp:

2006-02-10T13:43:41Z (19 years ago)

Author:

Josef Cejka <malyzelenyhnus@…>

Branches:

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

Children:

d3ca210

Parents:

e6a40ac

Message:

Fixed some problems with 64 bit arithmetic but others still persisting.

File:

• softfloat/generic/div.c (modified) (2 diffs)

softfloat/generic/div.c

@@ -32 +32 @@
 #include<comparison.h>
 #include<mul.h>
+#include<common.h>
+
 
 float32 divFloat32(float32 a, float32 b) 
@@ -307 +309 @@
         }
         
-        /* pack and round */
-        
-        /* find first nonzero digit and shift result and detect possibly underflow */
-        while ((cexp > 0) && (cfrac) && (!(cfrac & (FLOAT64_HIDDEN_BIT_MASK << (64 - FLOAT64_FRACTION_SIZE - 1 ) )))) {
-                cexp--; 
-                cfrac <<= 1;
-                        /* TODO: fix underflow */
-        };
-        
-        
-        cfrac >>= 1;
-        ++cexp;
-        cfrac += (0x1 << (64 - FLOAT64_FRACTION_SIZE - 3)); 
-
-        if (cfrac & (FLOAT64_HIDDEN_BIT_MASK << (64 - FLOAT64_FRACTION_SIZE - 1 ))) {
-                ++cexp;
-                cfrac >>= 1;
-                }       
-
-        /* check overflow */
-        if (cexp >= FLOAT64_MAX_EXPONENT ) {
-                /* FIXME: overflow, return infinity */
-                result.parts.exp = FLOAT64_MAX_EXPONENT;
-                result.parts.fraction = 0;
-                return result;
-        }
-
-        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.fraction = ((cfrac >>(64 - FLOAT64_FRACTION_SIZE - 2 ) ) & (~FLOAT64_HIDDEN_BIT_MASK)); 
-        
-        return result;  
+        /* round and shift */
+        result = finishFloat64(cexp, cfrac, result.parts.sign);
+        return result;
+
 }