WebSVN – HelenOS-historic – Diff – /uspace/trunk/softfloat/generic/add.c

 /** Add two Float32 numbers with same signs
  */
 float32 addFloat32(float32 a, float32 b)
+{
     int expdiff;
-    __u32 exp1, exp2,mant1, mant2;
+    __u32 exp1, exp2,frac1, frac2;
     expdiff = a.parts.exp - b.parts.exp;
     if (expdiff < 0) {
         if (isFloat32NaN(b)) {
-            //TODO: fix SigNaN
+            /* TODO: fix SigNaN */
             if (isFloat32SigNaN(b)) {
             };
             return b;
         };
         if (b.parts.exp == FLOAT32_MAX_EXPONENT) {
             return b;
+        }
-        mant1 = b.parts.mantisa;
+        frac1 = b.parts.fraction;
         exp1 = b.parts.exp;
-        mant2 = a.parts.mantisa;
+        frac2 = a.parts.fraction;
         exp2 = a.parts.exp;
         expdiff *= -1;
     } else {
-        if (isFloat32NaN(a)) {
+        if ((isFloat32NaN(a)) || (isFloat32NaN(b))) {
-            //TODO: fix SigNaN
+            /* TODO: fix SigNaN */
             if (isFloat32SigNaN(a) || isFloat32SigNaN(b)) {
             };
-            return a;
+            return (isFloat32NaN(a)?a:b);
         };
         if (a.parts.exp == FLOAT32_MAX_EXPONENT) {
             return a;
+        }
-        mant1 = a.parts.mantisa;
+        frac1 = a.parts.fraction;
         exp1 = a.parts.exp;
-        mant2 = b.parts.mantisa;
+        frac2 = b.parts.fraction;
         exp2 = b.parts.exp;
     };
     if (exp1 == 0) {
         /* both are denormalized */
-        mant1 += mant2;
+        frac1 += frac2;
-        if (mant1 & FLOAT32_HIDDEN_BIT_MASK ) {
+        if (frac1 & FLOAT32_HIDDEN_BIT_MASK ) {
             /* result is not denormalized */
             a.parts.exp = 1;
         };
-        a.parts.mantisa = mant1;
+        a.parts.fraction = frac1;
         return a;
     };
-    mant1 |= FLOAT32_HIDDEN_BIT_MASK; //add hidden bit
+    frac1 |= FLOAT32_HIDDEN_BIT_MASK; /* add hidden bit */
     if (exp2 == 0) {
         /* second operand is denormalized */
         --expdiff;
     } else {
         /* add hidden bit to second operand */
-        mant2 |= FLOAT32_HIDDEN_BIT_MASK;
+        frac2 |= FLOAT32_HIDDEN_BIT_MASK;
     };
     /* create some space for rounding */
-    mant1 <<= 6;
+    frac1 <<= 6;
-    mant2 <<= 6;
+    frac2 <<= 6;
+    if (expdiff < (FLOAT32_FRACTION_SIZE + 2) ) {
+        frac2 >>= expdiff;
+        frac1 += frac2;
+        };
-    if (expdiff > (FLOAT32_MANTISA_SIZE + 1) ) {
-         goto done;
-         };
-    mant2 >>= expdiff;
-    mant1 += mant2;
-done:
-    if (mant1 & (FLOAT32_HIDDEN_BIT_MASK << 7) ) {
+    if (frac1 & (FLOAT32_HIDDEN_BIT_MASK << 7) ) {
         ++exp1;
-        mant1 >>= 1;
+        frac1 >>= 1;
     };
-    /* rounding - if first bit after mantisa is set then round up */
+    /* rounding - if first bit after fraction is set then round up */
-    mant1 += (0x1 << 5);
+    frac1 += (0x1 << 5);
-    if (mant1 & (FLOAT32_HIDDEN_BIT_MASK << 7)) {
+    if (frac1 & (FLOAT32_HIDDEN_BIT_MASK << 7)) {
+        /* rounding overflow */
         ++exp1;
-        mant1 >>= 1;
+        frac1 >>= 1;
     };
+    if ((a.parts.exp == FLOAT32_MAX_EXPONENT ) || (a.parts.exp < exp1)) {
+            /* overflow - set infinity as result */
+            a.parts.exp = FLOAT32_MAX_EXPONENT;
+            a.parts.fraction = 0;
+            return a;
+            }
     a.parts.exp = exp1;
     /*Clear hidden bit and shift */
-    a.parts.mantisa = ((mant1 >> 6) & (~FLOAT32_HIDDEN_BIT_MASK)) ;
+    a.parts.fraction = ((frac1 >> 6) & (~FLOAT32_HIDDEN_BIT_MASK)) ;
     return a;
+}
 /** Add two Float64 numbers with same signs
  */
 float64 addFloat64(float64 a, float64 b)
+{
     int expdiff;
     __u32 exp1, exp2;
-    __u64 mant1, mant2;
+    __u64 frac1, frac2;
     expdiff = a.parts.exp - b.parts.exp;
     if (expdiff < 0) {
         if (isFloat64NaN(b)) {
-            //TODO: fix SigNaN
+            /* TODO: fix SigNaN */
             if (isFloat64SigNaN(b)) {
             };
             return b;
         };
         /* b is infinity and a not */
         if (b.parts.exp == FLOAT64_MAX_EXPONENT ) {
             return b;
+        }
-        mant1 = b.parts.mantisa;
+        frac1 = b.parts.fraction;
         exp1 = b.parts.exp;
-        mant2 = a.parts.mantisa;
+        frac2 = a.parts.fraction;
         exp2 = a.parts.exp;
         expdiff *= -1;
     } else {
         if (isFloat64NaN(a)) {
-            //TODO: fix SigNaN
+            /* TODO: fix SigNaN */
             if (isFloat64SigNaN(a) || isFloat64SigNaN(b)) {
             };
             return a;
         };
         /* a is infinity and b not */
         if (a.parts.exp == FLOAT64_MAX_EXPONENT ) {
             return a;
+        }
-        mant1 = a.parts.mantisa;
+        frac1 = a.parts.fraction;
         exp1 = a.parts.exp;
-        mant2 = b.parts.mantisa;
+        frac2 = b.parts.fraction;
         exp2 = b.parts.exp;
     };
     if (exp1 == 0) {
         /* both are denormalized */
-        mant1 += mant2;
+        frac1 += frac2;
-        if (mant1 & FLOAT64_HIDDEN_BIT_MASK) {
+        if (frac1 & FLOAT64_HIDDEN_BIT_MASK) {
             /* result is not denormalized */
             a.parts.exp = 1;
         };
-        a.parts.mantisa = mant1;
+        a.parts.fraction = frac1;
         return a;
     };
-    /* add hidden bit - mant1 is sure not denormalized */
+    /* add hidden bit - frac1 is sure not denormalized */
-    mant1 |= FLOAT64_HIDDEN_BIT_MASK;
+    frac1 |= FLOAT64_HIDDEN_BIT_MASK;
     /* second operand ... */
     if (exp2 == 0) {
         /* ... is denormalized */
         --expdiff;
     } else {
         /* is not denormalized */
-        mant2 |= FLOAT64_HIDDEN_BIT_MASK;
+        frac2 |= FLOAT64_HIDDEN_BIT_MASK;
     };
     /* create some space for rounding */
-    mant1 <<= 6;
+    frac1 <<= 6;
-    mant2 <<= 6;
+    frac2 <<= 6;
-    if (expdiff > (FLOAT64_MANTISA_SIZE + 1) ) {
+    if (expdiff < (FLOAT64_FRACTION_SIZE + 2) ) {
+        frac2 >>= expdiff;
-         goto done;
+        frac1 += frac2;
-         };
+        };
-    mant2 >>= expdiff;
-    mant1 += mant2;
-done:
-    if (mant1 & (FLOAT64_HIDDEN_BIT_MASK << 7) ) {
+    if (frac1 & (FLOAT64_HIDDEN_BIT_MASK << 7) ) {
         ++exp1;
-        mant1 >>= 1;
+        frac1 >>= 1;
     };
-    /* rounding - if first bit after mantisa is set then round up */
+    /* rounding - if first bit after fraction is set then round up */
-    mant1 += (0x1 << 5);
+    frac1 += (0x1 << 5);
-    if (mant1 & (FLOAT64_HIDDEN_BIT_MASK << 7)) {
+    if (frac1 & (FLOAT64_HIDDEN_BIT_MASK << 7)) {
+        /* rounding overflow */
         ++exp1;
-        mant1 >>= 1;
+        frac1 >>= 1;
     };
+    if ((a.parts.exp == FLOAT64_MAX_EXPONENT ) || (a.parts.exp < exp1)) {
+            /* overflow - set infinity as result */
+            a.parts.exp = FLOAT64_MAX_EXPONENT;
+            a.parts.fraction = 0;
+            return a;
+            }
     a.parts.exp = exp1;
     /*Clear hidden bit and shift */
-    a.parts.mantisa = ( (mant1 >> 6 ) & (~FLOAT64_HIDDEN_BIT_MASK));
+    a.parts.fraction = ( (frac1 >> 6 ) & (~FLOAT64_HIDDEN_BIT_MASK));
     return a;
+}

Subversion Repositories HelenOS-historic

(root)/uspace/trunk/softfloat/generic/add.c – Rev 734 → 804