Changeset d3e241a in mainline for uspace/lib/softfloat/generic/softfloat.c

Timestamp:

2011-08-06T07:16:50Z (13 years ago)

Author:

Petr Koupy <petr.koupy@…>

Branches:

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

Children:

b705ecc, cc677f1

Parents:

b7ee0369 (diff), c67aff2 (diff)
Note: this is a merge changeset, the changes displayed below correspond to the merge itself.
Use the (diff) links above to see all the changes relative to each parent.

Message:

Merge libposix changes.

File:

• uspace/lib/softfloat/generic/softfloat.c (modified) (40 diffs)

uspace/lib/softfloat/generic/softfloat.c

@@ -1 +1 @@
 /*
  * Copyright (c) 2005 Josef Cejka
+ * Copyright (c) 2011 Petr Koupy
  * All rights reserved.
  *
@@ -32 +33 @@
  * @{
  */
-/** @file
+/** @file Softfloat API.
  */
 
@@ -83 +84 @@
 }
 
+long double __addtf3(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+        if (ta.parts.sign != tb.parts.sign) {
+                if (ta.parts.sign) {
+                        ta.parts.sign = 0;
+                        return subFloat128(tb, ta).ld;
+                };
+                tb.parts.sign = 0;
+                return subFloat128(ta, tb).ld;
+        }
+        return addFloat128(ta, tb).ld;
+}
+
 float __subsf3(float a, float b)
 {
@@ -107 +124 @@
 }
 
+long double __subtf3(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+        if (ta.parts.sign != tb.parts.sign) {
+                tb.parts.sign = !tb.parts.sign;
+                return addFloat128(ta, tb).ld;
+        }
+        return subFloat128(ta, tb).ld;
+}
+
 float __mulsf3(float a, float b) 
 {
@@ -123 +152 @@
 }
 
+long double __multf3(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+        return  mulFloat128(ta, tb).ld;
+}
+
 float __divsf3(float a, float b) 
 {
@@ -139 +176 @@
 }
 
+long double __divtf3(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+        return  divFloat128(ta, tb).ld;
+}
+
 float __negsf2(float a)
 {
@@ -149 +194 @@
 double __negdf2(double a)
 {
-        float64 fa;
-        fa.d = a;
-        fa.parts.sign = !fa.parts.sign;
-        return fa.d;
+        float64 da;
+        da.d = a;
+        da.parts.sign = !da.parts.sign;
+        return da.d;
+}
+
+long double __negtf2(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+        ta.parts.sign = !ta.parts.sign;
+        return ta.ld;
 }
 
@@ -164 +217 @@
 }
 
+long double __extendsftf2(float a)
+{
+        float32 fa;
+        fa.f = a;
+        return convertFloat32ToFloat128(fa).ld;
+}
+
+long double __extenddftf2(double a)
+{
+        float64 da;
+        da.d = a;
+        return convertFloat64ToFloat128(da).ld;
+}
+
 float __truncdfsf2(double a) 
 {
@@ -171 +238 @@
 }
 
+float __trunctfsf2(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+        return convertFloat128ToFloat32(ta).f;
+}
+
+double __trunctfdf2(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+        return convertFloat128ToFloat64(ta).d;
+}
+
 int __fixsfsi(float a)
 {
@@ -178 +259 @@
         return float32_to_int(fa);
 }
+
 int __fixdfsi(double a)
 {
@@ -184 +266 @@
         
         return float64_to_int(da);
+}
+
+int __fixtfsi(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+
+        return float128_to_int(ta);
 }
  
@@ -193 +283 @@
         return float32_to_long(fa);
 }
+
 long __fixdfdi(double a)
 {
@@ -199 +290 @@
         
         return float64_to_long(da);
+}
+
+long __fixtfdi(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+
+        return float128_to_long(ta);
 }
  
@@ -208 +307 @@
         return float32_to_longlong(fa);
 }
+
 long long __fixdfti(double a)
 {
@@ -216 +316 @@
 }
 
+long long __fixtfti(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+
+        return float128_to_longlong(ta);
+}
+
 unsigned int __fixunssfsi(float a)
 {
@@ -223 +331 @@
         return float32_to_uint(fa);
 }
+
 unsigned int __fixunsdfsi(double a)
 {
@@ -229 +338 @@
         
         return float64_to_uint(da);
+}
+
+unsigned int __fixunstfsi(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+
+        return float128_to_uint(ta);
 }
  
@@ -238 +355 @@
         return float32_to_ulong(fa);
 }
+
 unsigned long __fixunsdfdi(double a)
 {
@@ -244 +362 @@
         
         return float64_to_ulong(da);
+}
+
+unsigned long __fixunstfdi(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+
+        return float128_to_ulong(ta);
 }
  
@@ -253 +379 @@
         return float32_to_ulonglong(fa);
 }
+
 unsigned long long __fixunsdfti(double a)
 {
@@ -259 +386 @@
         
         return float64_to_ulonglong(da);
+}
+
+unsigned long long __fixunstfti(long double a)
+{
+        float128 ta;
+        ta.ld = a;
+
+        return float128_to_ulonglong(ta);
 }
  
@@ -268 +403 @@
         return fa.f;
 }
+
 double __floatsidf(int i)
 {
@@ -275 +411 @@
         return da.d;
 }
+
+long double __floatsitf(int i)
+{
+        float128 ta;
+
+        ta = int_to_float128(i);
+        return ta.ld;
+}
  
 float __floatdisf(long i)
@@ -283 +427 @@
         return fa.f;
 }
+
 double __floatdidf(long i)
 {
@@ -290 +435 @@
         return da.d;
 }
+
+long double __floatditf(long i)
+{
+        float128 ta;
+
+        ta = long_to_float128(i);
+        return ta.ld;
+}
  
 float __floattisf(long long i)
@@ -298 +451 @@
         return fa.f;
 }
+
 double __floattidf(long long i)
 {
@@ -306 +460 @@
 }
 
+long double __floattitf(long long i)
+{
+        float128 ta;
+
+        ta = longlong_to_float128(i);
+        return ta.ld;
+}
+
 float __floatunsisf(unsigned int i)
 {
@@ -313 +475 @@
         return fa.f;
 }
+
 double __floatunsidf(unsigned int i)
 {
@@ -320 +483 @@
         return da.d;
 }
+
+long double __floatunsitf(unsigned int i)
+{
+        float128 ta;
+
+        ta = uint_to_float128(i);
+        return ta.ld;
+}
  
 float __floatundisf(unsigned long i)
@@ -328 +499 @@
         return fa.f;
 }
+
 double __floatundidf(unsigned long i)
 {
@@ -335 +507 @@
         return da.d;
 }
+
+long double __floatunditf(unsigned long i)
+{
+        float128 ta;
+
+        ta = ulong_to_float128(i);
+        return ta.ld;
+}
  
 float __floatuntisf(unsigned long long i)
@@ -343 +523 @@
         return fa.f;
 }
+
 double __floatuntidf(unsigned long long i)
 {
@@ -351 +532 @@
 }
 
+long double __floatuntitf(unsigned long long i)
+{
+        float128 ta;
+
+        ta = ulonglong_to_float128(i);
+        return ta.ld;
+}
+
 /* Comparison functions */
-/* Comparison functions */
-
-/* a<b .. -1
- * a=b ..  0
- * a>b ..  1
- * */
 
 int __cmpsf2(float a, float b) 
@@ -364 +547 @@
         fa.f = a;
         fb.f = b;
-        if ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) ) {
+
+        if ((isFloat32NaN(fa)) || (isFloat32NaN(fb))) {
                 return 1; /* no special constant for unordered - maybe signaled? */
-        };
-
+        }
         
         if (isFloat32eq(fa, fb)) {
                 return 0;
-        };
+        }
         
         if (isFloat32lt(fa, fb)) {
                 return -1;
-                };
+        }
+
         return 1;
 }
 
+int __cmpdf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+
+        if ((isFloat64NaN(da)) || (isFloat64NaN(db))) {
+                return 1; /* no special constant for unordered - maybe signaled? */
+        }
+
+        if (isFloat64eq(da, db)) {
+                return 0;
+        }
+
+        if (isFloat64lt(da, db)) {
+                return -1;
+        }
+
+        return 1;
+}
+
+int __cmptf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 1; /* no special constant for unordered - maybe signaled? */
+        }
+
+        if (isFloat128eq(ta, tb)) {
+                return 0;
+        }
+
+        if (isFloat128lt(ta, tb)) {
+                return -1;
+        }
+
+        return 1;
+}
+
 int __unordsf2(float a, float b) 
 {
@@ -384 +610 @@
         fa.f = a;
         fb.f = b;
-        return ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) );
-}
-
-/** 
- * @return zero, if neither argument is a NaN and are equal
- * */
+        return ((isFloat32NaN(fa)) || (isFloat32NaN(fb)));
+}
+
+int __unorddf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+        return ((isFloat64NaN(da)) || (isFloat64NaN(db)));
+}
+
+int __unordtf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+        return ((isFloat128NaN(ta)) || (isFloat128NaN(tb)));
+}
+
 int __eqsf2(float a, float b) 
 {
@@ -395 +634 @@
         fa.f = a;
         fb.f = b;
-        if ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) ) {
-                /* TODO: sigNaNs*/
-                return 1;
-                };
+        if ((isFloat32NaN(fa)) || (isFloat32NaN(fb))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
         return isFloat32eq(fa, fb) - 1;
 }
 
-/* strange behavior, but it was in gcc documentation */
+int __eqdf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+        if ((isFloat64NaN(da)) || (isFloat64NaN(db))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+        return isFloat64eq(da, db) - 1;
+}
+
+int __eqtf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+        return isFloat128eq(ta, tb) - 1;
+}
+
 int __nesf2(float a, float b) 
 {
+        /* strange behavior, but it was in gcc documentation */
         return __eqsf2(a, b);
 }
 
-/* return value >= 0 if a>=b and neither is NaN */
+int __nedf2(double a, double b)
+{
+        /* strange behavior, but it was in gcc documentation */
+        return __eqdf2(a, b);
+}
+
+int __netf2(long double a, long double b)
+{
+        /* strange behavior, but it was in gcc documentation */
+        return __eqtf2(a, b);
+}
+
 int __gesf2(float a, float b)
 {
@@ -414 +688 @@
         fa.f = a;
         fb.f = b;
-        if ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) ) {
-                /* TODO: sigNaNs*/
-                return -1;
-                };
+
+        if ((isFloat32NaN(fa)) || (isFloat32NaN(fb))) {
+                /* TODO: sigNaNs */
+                return -1;
+        }
         
         if (isFloat32eq(fa, fb)) {
                 return 0;
-        };
+        }
         
         if (isFloat32gt(fa, fb)) {
                 return 1;
-                };
+        }
         
         return -1;
 }
 
-/** Return negative value, if a<b and neither is NaN*/
+int __gedf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+
+        if ((isFloat64NaN(da)) || (isFloat64NaN(db))) {
+                /* TODO: sigNaNs */
+                return -1;
+        }
+
+        if (isFloat64eq(da, db)) {
+                return 0;
+        }
+
+        if (isFloat64gt(da, db)) {
+                return 1;
+        }
+
+        return -1;
+}
+
+int __getf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                /* TODO: sigNaNs */
+                return -1;
+        }
+
+        if (isFloat128eq(ta, tb)) {
+                return 0;
+        }
+
+        if (isFloat128gt(ta, tb)) {
+                return 1;
+        }
+
+        return -1;
+}
+
 int __ltsf2(float a, float b)
 {
@@ -436 +754 @@
         fa.f = a;
         fb.f = b;
-        if ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) ) {
-                /* TODO: sigNaNs*/
-                return 1;
-                };
+
+        if ((isFloat32NaN(fa)) || (isFloat32NaN(fb))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+
         if (isFloat32lt(fa, fb)) {
                 return -1;
-                };
+        }
+
         return 0;
 }
 
-/* return value <= 0 if a<=b and neither is NaN */
+int __ltdf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+
+        if ((isFloat64NaN(da)) || (isFloat64NaN(db))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+
+        if (isFloat64lt(da, db)) {
+                return -1;
+        }
+
+        return 0;
+}
+
+int __lttf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+
+        if (isFloat128lt(ta, tb)) {
+                return -1;
+        }
+
+        return 0;
+}
+
 int __lesf2(float a, float b)
 {
@@ -452 +808 @@
         fa.f = a;
         fb.f = b;
-        if ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) ) {
-                /* TODO: sigNaNs*/
-                return 1;
-                };
+
+        if ((isFloat32NaN(fa)) || (isFloat32NaN(fb))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
         
         if (isFloat32eq(fa, fb)) {
                 return 0;
-        };
+        }
         
         if (isFloat32lt(fa, fb)) {
                 return -1;
-                };
+        }
         
         return 1;
 }
 
-/** Return positive value, if a>b and neither is NaN*/
+int __ledf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+
+        if ((isFloat64NaN(da)) || (isFloat64NaN(db))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+
+        if (isFloat64eq(da, db)) {
+                return 0;
+        }
+
+        if (isFloat64lt(da, db)) {
+                return -1;
+        }
+
+        return 1;
+}
+
+int __letf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                /* TODO: sigNaNs */
+                return 1;
+        }
+
+        if (isFloat128eq(ta, tb)) {
+                return 0;
+        }
+
+        if (isFloat128lt(ta, tb)) {
+                return -1;
+        }
+
+        return 1;
+}
+
 int __gtsf2(float a, float b)
 {
@@ -474 +874 @@
         fa.f = a;
         fb.f = b;
-        if ( (isFloat32NaN(fa)) || (isFloat32NaN(fb)) ) {
-                /* TODO: sigNaNs*/
-                return -1;
-                };
+
+        if ((isFloat32NaN(fa)) || (isFloat32NaN(fb))) {
+                /* TODO: sigNaNs */
+                return -1;
+        }
+
         if (isFloat32gt(fa, fb)) {
                 return 1;
-                };
+        }
+
         return 0;
 }
 
-/* Other functions */
-
-float __powisf2(float a, int b)
-{
-/* TODO: */
-        float32 fa;
-        fa.binary = FLOAT32_NAN;
-        return fa.f;
-}
+int __gtdf2(double a, double b)
+{
+        float64 da, db;
+        da.d = a;
+        db.d = b;
+
+        if ((isFloat64NaN(da)) || (isFloat64NaN(db))) {
+                /* TODO: sigNaNs */
+                return -1;
+        }
+
+        if (isFloat64gt(da, db)) {
+                return 1;
+        }
+
+        return 0;
+}
+
+int __gttf2(long double a, long double b)
+{
+        float128 ta, tb;
+        ta.ld = a;
+        tb.ld = b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                /* TODO: sigNaNs */
+                return -1;
+        }
+
+        if (isFloat128gt(ta, tb)) {
+                return 1;
+        }
+
+        return 0;
+}
+
+
+
+#ifdef SPARC_SOFTFLOAT
+
+/* SPARC quadruple-precision wrappers */
+
+void _Qp_add(long double *c, long double *a, long double *b)
+{
+        *c = __addtf3(*a, *b);
+}
+
+void _Qp_sub(long double *c, long double *a, long double *b)
+{
+        *c = __subtf3(*a, *b);
+}
+
+void _Qp_mul(long double *c, long double *a, long double *b)
+{
+        *c = __multf3(*a, *b);
+}
+
+void _Qp_div(long double *c, long double *a, long double *b)
+{
+        *c = __divtf3(*a, *b);
+}
+
+void _Qp_neg(long double *c, long double *a)
+{
+        *c = __negtf2(*a);
+}
+
+void _Qp_stoq(long double *c, float a)
+{
+        *c = __extendsftf2(a);
+}
+
+void _Qp_dtoq(long double *c, double a)
+{
+        *c = __extenddftf2(a);
+}
+
+float _Qp_qtos(long double *a)
+{
+        return __trunctfsf2(*a);
+}
+
+double _Qp_qtod(long double *a)
+{
+        return __trunctfdf2(*a);
+}
+
+int _Qp_qtoi(long double *a)
+{
+        return __fixtfsi(*a);
+}
+
+unsigned int _Qp_qtoui(long double *a)
+{
+        return __fixunstfsi(*a);
+}
+
+long _Qp_qtox(long double *a)
+{
+        return __fixtfdi(*a);
+}
+
+unsigned long _Qp_qtoux(long double *a)
+{
+        return __fixunstfdi(*a);
+}
+
+void _Qp_itoq(long double *c, int a)
+{
+        *c = __floatsitf(a);
+}
+
+void _Qp_uitoq(long double *c, unsigned int a)
+{
+        *c = __floatunsitf(a);
+}
+
+void _Qp_xtoq(long double *c, long a)
+{
+        *c = __floatditf(a);
+}
+
+void _Qp_uxtoq(long double *c, unsigned long a)
+{
+        *c = __floatunditf(a);
+}
+
+int _Qp_cmp(long double *a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 3;
+        }
+
+        if (isFloat128eq(ta, tb)) {
+                return 0;
+        }
+
+        if (isFloat128lt(ta, tb)) {
+                return 1;
+        }
+
+        return 2;
+}
+
+int _Qp_cmpe(long double *a, long double *b)
+{
+        /* strange, but is defined this way in SPARC Compliance Definition */
+        return _Qp_cmp(a, b);
+}
+
+int _Qp_feq(long double *a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 0;
+        }
+
+        return isFloat128eq(ta, tb);
+}
+
+int _Qp_fge(long double *a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 0;
+        }
+
+        return isFloat128eq(ta, tb) || isFloat128gt(ta, tb);
+}
+
+int _Qp_fgt(long double *a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 0;
+        }
+
+        return isFloat128gt(ta, tb);
+}
+
+int _Qp_fle(long double*a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 0;
+        }
+
+        return isFloat128eq(ta, tb) || isFloat128lt(ta, tb);
+}
+
+int _Qp_flt(long double *a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 0;
+        }
+
+        return isFloat128lt(ta, tb);
+}
+
+int _Qp_fne(long double *a, long double *b)
+{
+        float128 ta, tb;
+        ta.ld = *a;
+        tb.ld = *b;
+
+        if ((isFloat128NaN(ta)) || (isFloat128NaN(tb))) {
+                return 1;
+        }
+
+        return !isFloat128eq(ta, tb);
+}
+
+#endif
 
 /** @}