Changes in uspace/lib/softfloat/sub.c [a35b458:c0c38c7c] in mainline

File:

• uspace/lib/softfloat/sub.c (modified) (39 diffs)

uspace/lib/softfloat/sub.c

@@ -52 +52 @@
         uint32_t exp1, exp2, frac1, frac2;
         float32 result;
-
+        
         result.bin = 0;
-
+        
         expdiff = a.parts.exp - b.parts.exp;
         if ((expdiff < 0 ) || ((expdiff == 0) &&
@@ -62 +62 @@
                                 // TODO: fix SigNaN
                         }
-
+                        
                         return b;
                 }
-
+                
                 if (b.parts.exp == FLOAT32_MAX_EXPONENT) {
                         /* num -(+-inf) = -+inf */
@@ -71 +71 @@
                         return b;
                 }
-
+                
                 result.parts.sign = !a.parts.sign;
-
+                
                 frac1 = b.parts.fraction;
                 exp1 = b.parts.exp;
@@ -84 +84 @@
                                 // TODO: fix SigNaN
                         }
-
+                        
                         return a;
                 }
-
+                
                 if (a.parts.exp == FLOAT32_MAX_EXPONENT) {
                         if (b.parts.exp == FLOAT32_MAX_EXPONENT) {
@@ -95 +95 @@
                                 return result;
                         }
-
+                        
                         return a;
                 }
-
+                
                 result.parts.sign = a.parts.sign;
-
+                
                 frac1 = a.parts.fraction;
                 exp1 = a.parts.exp;
@@ -106 +106 @@
                 exp2 = b.parts.exp;
         }
-
+        
         if (exp1 == 0) {
                 /* both are denormalized */
@@ -114 +114 @@
                         return result;
                 }
-
+                
                 result.parts.exp = 0;
                 return result;
         }
-
+        
         /* add hidden bit */
         frac1 |= FLOAT32_HIDDEN_BIT_MASK;
-
+        
         if (exp2 == 0) {
                 /* denormalized */
@@ -129 +129 @@
                 frac2 |= FLOAT32_HIDDEN_BIT_MASK;
         }
-
+        
         /* create some space for rounding */
         frac1 <<= 6;
         frac2 <<= 6;
-
+        
         if (expdiff > FLOAT32_FRACTION_SIZE + 1)
                 goto done;
-
+        
         frac1 = frac1 - (frac2 >> expdiff);
-
+        
 done:
         /* TODO: find first nonzero digit and shift result and detect possibly underflow */
@@ -146 +146 @@
                 /* TODO: fix underflow - frac1 == 0 does not necessary means underflow... */
         }
-
+        
         /* rounding - if first bit after fraction is set then round up */
         frac1 += 0x20;
-
+        
         if (frac1 & (FLOAT32_HIDDEN_BIT_MASK << 7)) {
                 ++exp1;
                 frac1 >>= 1;
         }
-
+        
         /* Clear hidden bit and shift */
         result.parts.fraction = ((frac1 >> 6) & (~FLOAT32_HIDDEN_BIT_MASK));
         result.parts.exp = exp1;
-
+        
         return result;
 }
@@ -176 +176 @@
         uint64_t frac1, frac2;
         float64 result;
-
+        
         result.bin = 0;
-
+        
         expdiff = a.parts.exp - b.parts.exp;
         if ((expdiff < 0 ) ||
@@ -186 +186 @@
                                 // TODO: fix SigNaN
                         }
-
+                        
                         return b;
                 }
-
+                
                 if (b.parts.exp == FLOAT64_MAX_EXPONENT) {
                         /* num -(+-inf) = -+inf */
@@ -195 +195 @@
                         return b;
                 }
-
+                
                 result.parts.sign = !a.parts.sign;
-
+                
                 frac1 = b.parts.fraction;
                 exp1 = b.parts.exp;
@@ -208 +208 @@
                                 // TODO: fix SigNaN
                         }
-
+                        
                         return a;
                 }
-
+                
                 if (a.parts.exp == FLOAT64_MAX_EXPONENT) {
                         if (b.parts.exp == FLOAT64_MAX_EXPONENT) {
@@ -219 +219 @@
                                 return result;
                         }
-
+                        
                         return a;
                 }
-
+                
                 result.parts.sign = a.parts.sign;
-
+                
                 frac1 = a.parts.fraction;
                 exp1 = a.parts.exp;
@@ -230 +230 @@
                 exp2 = b.parts.exp;
         }
-
+        
         if (exp1 == 0) {
                 /* both are denormalized */
@@ -238 +238 @@
                         return result;
                 }
-
+                
                 result.parts.exp = 0;
                 return result;
         }
-
+        
         /* add hidden bit */
         frac1 |= FLOAT64_HIDDEN_BIT_MASK;
-
+        
         if (exp2 == 0) {
                 /* denormalized */
@@ -253 +253 @@
                 frac2 |= FLOAT64_HIDDEN_BIT_MASK;
         }
-
+        
         /* create some space for rounding */
         frac1 <<= 6;
         frac2 <<= 6;
-
+        
         if (expdiff > FLOAT64_FRACTION_SIZE + 1)
                 goto done;
-
+        
         frac1 = frac1 - (frac2 >> expdiff);
-
+        
 done:
         /* TODO: find first nonzero digit and shift result and detect possibly underflow */
@@ -270 +270 @@
                 /* TODO: fix underflow - frac1 == 0 does not necessary means underflow... */
         }
-
+        
         /* rounding - if first bit after fraction is set then round up */
         frac1 += 0x20;
-
+        
         if (frac1 & (FLOAT64_HIDDEN_BIT_MASK << 7)) {
                 ++exp1;
                 frac1 >>= 1;
         }
-
+        
         /* Clear hidden bit and shift */
         result.parts.fraction = ((frac1 >> 6) & (~FLOAT64_HIDDEN_BIT_MASK));
         result.parts.exp = exp1;
-
+        
         return result;
 }
@@ -300 +300 @@
         uint64_t frac1_hi, frac1_lo, frac2_hi, frac2_lo, tmp_hi, tmp_lo;
         float128 result;
-
+        
         result.bin.hi = 0;
         result.bin.lo = 0;
-
+        
         expdiff = a.parts.exp - b.parts.exp;
         if ((expdiff < 0 ) || ((expdiff == 0) &&
@@ -311 +311 @@
                                 // TODO: fix SigNaN
                         }
-
+                        
                         return b;
                 }
-
+                
                 if (b.parts.exp == FLOAT128_MAX_EXPONENT) {
                         /* num -(+-inf) = -+inf */
@@ -320 +320 @@
                         return b;
                 }
-
+                
                 result.parts.sign = !a.parts.sign;
-
+                
                 frac1_hi = b.parts.frac_hi;
                 frac1_lo = b.parts.frac_lo;
@@ -335 +335 @@
                                 // TODO: fix SigNaN
                         }
-
+                        
                         return a;
                 }
-
+                
                 if (a.parts.exp == FLOAT128_MAX_EXPONENT) {
                         if (b.parts.exp == FLOAT128_MAX_EXPONENT) {
@@ -349 +349 @@
                         return a;
                 }
-
+                
                 result.parts.sign = a.parts.sign;
-
+                
                 frac1_hi = a.parts.frac_hi;
                 frac1_lo = a.parts.frac_lo;
@@ -359 +359 @@
                 exp2 = b.parts.exp;
         }
-
+        
         if (exp1 == 0) {
                 /* both are denormalized */
@@ -369 +369 @@
                         return result;
                 }
-
+                
                 result.parts.exp = 0;
                 return result;
         }
-
+        
         /* add hidden bit */
         or128(frac1_hi, frac1_lo,
             FLOAT128_HIDDEN_BIT_MASK_HI, FLOAT128_HIDDEN_BIT_MASK_LO,
             &frac1_hi, &frac1_lo);
-
+        
         if (exp2 == 0) {
                 /* denormalized */
@@ -388 +388 @@
                     &frac2_hi, &frac2_lo);
         }
-
+        
         /* create some space for rounding */
         lshift128(frac1_hi, frac1_lo, 6, &frac1_hi, &frac1_lo);
         lshift128(frac2_hi, frac2_lo, 6, &frac2_hi, &frac2_lo);
-
+        
         if (expdiff > FLOAT128_FRACTION_SIZE + 1)
                 goto done;
-
+        
         rshift128(frac2_hi, frac2_lo, expdiff, &tmp_hi, &tmp_lo);
         sub128(frac1_hi, frac1_lo, tmp_hi, tmp_lo, &frac1_hi, &frac1_lo);
-
+        
 done:
         /* TODO: find first nonzero digit and shift result and detect possibly underflow */
@@ -408 +408 @@
                 lshift128(frac1_hi, frac1_lo, 1, &frac1_hi, &frac1_lo);
                 /* TODO: fix underflow - frac1 == 0 does not necessary means underflow... */
-
+                
                 lshift128(FLOAT128_HIDDEN_BIT_MASK_HI, FLOAT128_HIDDEN_BIT_MASK_LO, 6,
                     &tmp_hi, &tmp_lo);
                 and128(frac1_hi, frac1_lo, tmp_hi, tmp_lo, &tmp_hi, &tmp_lo);
         }
-
+        
         /* rounding - if first bit after fraction is set then round up */
         add128(frac1_hi, frac1_lo, 0x0ll, 0x20ll, &frac1_hi, &frac1_lo);
-
+        
         lshift128(FLOAT128_HIDDEN_BIT_MASK_HI, FLOAT128_HIDDEN_BIT_MASK_LO, 7,
            &tmp_hi, &tmp_lo);
@@ -424 +424 @@
                 rshift128(frac1_hi, frac1_lo, 1, &frac1_hi, &frac1_lo);
         }
-
+        
         /* Clear hidden bit and shift */
         rshift128(frac1_hi, frac1_lo, 6, &frac1_hi, &frac1_lo);
@@ -432 +432 @@
         result.parts.frac_hi = tmp_hi;
         result.parts.frac_lo = tmp_lo;
-
+        
         result.parts.exp = exp1;
-
+        
         return result;
 }
@@ -444 +444 @@
         float32_u ua;
         ua.val = a;
-
+        
         float32_u ub;
         ub.val = b;
-
+        
         float32_u res;
-
+        
         if (ua.data.parts.sign != ub.data.parts.sign) {
                 ub.data.parts.sign = !ub.data.parts.sign;
@@ -455 +455 @@
         } else
                 res.data = sub_float32(ua.data, ub.data);
-
+        
         return res.val;
 }
@@ -463 +463 @@
         float32_u ua;
         ua.val = a;
-
+        
         float32_u ub;
         ub.val = b;
-
+        
         float32_u res;
-
+        
         if (ua.data.parts.sign != ub.data.parts.sign) {
                 ub.data.parts.sign = !ub.data.parts.sign;
@@ -474 +474 @@
         } else
                 res.data = sub_float32(ua.data, ub.data);
-
+        
         return res.val;
 }
@@ -486 +486 @@
         float64_u ua;
         ua.val = a;
-
+        
         float64_u ub;
         ub.val = b;
-
+        
         float64_u res;
-
+        
         if (ua.data.parts.sign != ub.data.parts.sign) {
                 ub.data.parts.sign = !ub.data.parts.sign;
@@ -497 +497 @@
         } else
                 res.data = sub_float64(ua.data, ub.data);
-
+        
         return res.val;
 }
@@ -505 +505 @@
         float64_u ua;
         ua.val = a;
-
+        
         float64_u ub;
         ub.val = b;
-
+        
         float64_u res;
-
+        
         if (ua.data.parts.sign != ub.data.parts.sign) {
                 ub.data.parts.sign = !ub.data.parts.sign;
@@ -516 +516 @@
         } else
                 res.data = sub_float64(ua.data, ub.data);
-
+        
         return res.val;
 }
@@ -528 +528 @@
         float128_u ua;
         ua.val = a;
-
+        
         float128_u ub;
         ub.val = b;
-
+        
         float128_u res;
-
+        
         if (ua.data.parts.sign != ub.data.parts.sign) {
                 ub.data.parts.sign = !ub.data.parts.sign;
@@ -539 +539 @@
         } else
                 res.data = sub_float128(ua.data, ub.data);
-
+        
         return res.val;
 }