Changeset c0c38c7c in mainline for uspace/lib/softfloat/mul.c

Timestamp:

2015-03-14T21:36:44Z (9 years ago)

Author:

Martin Decky <martin@…>

Branches:

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

Children:

7f9d97f3

Parents:

2c7fdaa

Message:

software floating point overhaul
use proper type mapping
fix cosine calculation

File:

• uspace/lib/softfloat/mul.c (modified) (9 diffs)

uspace/lib/softfloat/mul.c

@@ -34 +34 @@
  */
 
-#include "sftypes.h"
 #include "mul.h"
 #include "comparison.h"
@@ -62 +61 @@
                         return result;
                 }
+                
                 if (is_float32_signan(b)) { /* TODO: fix SigNaN */
                         result.parts.fraction = b.parts.fraction;
@@ -67 +67 @@
                         return result;
                 }
+                
                 /* set NaN as result */
                 result.bin = FLOAT32_NAN;
@@ -78 +79 @@
                         return result;
                 }
+                
                 result.parts.fraction = a.parts.fraction;
                 result.parts.exp = a.parts.exp;
@@ -89 +91 @@
                         return result;
                 }
+                
                 result.parts.fraction = b.parts.fraction;
                 result.parts.exp = b.parts.exp;
@@ -106 +109 @@
         }
         
-        if (exp < 0) { 
+        if (exp < 0) {
                 /* FIXME: underflow */
                 /* return signed zero */
@@ -164 +167 @@
                 /* denormalized number */
                 frac1 >>= 1; /* denormalize */
+                
                 while ((frac1 > 0) && (exp < 0)) {
                         frac1 >>= 1;
                         ++exp;
                 }
+                
                 if (frac1 == 0) {
                         /* FIXME : underflow */
@@ -175 +180 @@
                 }
         }
+        
         result.parts.exp = exp; 
         result.parts.fraction = frac1 & ((1 << FLOAT32_FRACTION_SIZE) - 1);
@@ -380 +386 @@
 }
 
+#ifdef float32_t
+
+float32_t __mulsf3(float32_t a, float32_t b)
+{
+        float32_u ua;
+        ua.val = a;
+        
+        float32_u ub;
+        ub.val = b;
+        
+        float32_u res;
+        res.data = mul_float32(ua.data, ub.data);
+        
+        return res.val;
+}
+
+float32_t __aeabi_fmul(float32_t a, float32_t b)
+{
+        float32_u ua;
+        ua.val = a;
+        
+        float32_u ub;
+        ub.val = b;
+        
+        float32_u res;
+        res.data = mul_float32(ua.data, ub.data);
+        
+        return res.val;
+}
+
+#endif
+
+#ifdef float64_t
+
+float64_t __muldf3(float64_t a, float64_t b)
+{
+        float64_u ua;
+        ua.val = a;
+        
+        float64_u ub;
+        ub.val = b;
+        
+        float64_u res;
+        res.data = mul_float64(ua.data, ub.data);
+        
+        return res.val;
+}
+
+float64_t __aeabi_dmul(float64_t a, float64_t b)
+{
+        float64_u ua;
+        ua.val = a;
+        
+        float64_u ub;
+        ub.val = b;
+        
+        float64_u res;
+        res.data = mul_float64(ua.data, ub.data);
+        
+        return res.val;
+}
+
+#endif
+
+#ifdef float128_t
+
+float128_t __multf3(float128_t a, float128_t b)
+{
+        float128_u ua;
+        ua.val = a;
+        
+        float128_u ub;
+        ub.val = b;
+        
+        float128_u res;
+        res.data = mul_float128(ua.data, ub.data);
+        
+        return res.val;
+}
+
+void _Qp_mul(float128_t *c, float128_t *a, float128_t *b)
+{
+        *c = __multf3(*a, *b);
+}
+
+#endif
+
 /** @}
  */