math: add dummy implementations of 128 bit long double functions

This is in preparation for the aarch64 port only to have the long
double math symbols available on ld128 platforms. The implementations
should be fixed up later once we have proper tests for these functions.

Added bigendian handling for ld128 bit manipulations too.

diff --git a/src/internal/libm.h b/src/internal/libm.h
index 88a7eb479e53e9253adbe08f79481d64f66946b6..df8641116af2ee442525828bfc27bc2de75ee508 100644 (file)
--- a/src/internal/libm.h
+++ b/src/internal/libm.h
@@ -42,6 +42,20 @@ union ldshape {
                uint64_t hi;
        } i2;
 };
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384 && __BYTE_ORDER == __BIG_ENDIAN
+union ldshape {
+       long double f;
+       struct {
+               uint16_t se;
+               uint16_t top;
+               uint32_t mid;
+               uint64_t lo;
+       } i;
+       struct {
+               uint64_t hi;
+               uint64_t lo;
+       } i2;
+};
 #else
 #error Unsupported long double representation
 #endif

diff --git a/src/math/acoshl.c b/src/math/acoshl.c
index 4aa84acb2ea3684b0ce4c3c2199b0ab6a6d4814e..8d4b43f64cdd899877bf39f43c7657f1f3b3f17c 100644 (file)
--- a/src/math/acoshl.c
+++ b/src/math/acoshl.c
@@ -20,4 +20,10 @@ long double acoshl(long double x)
                return logl(2*x - 1/(x+sqrtl(x*x-1)));
        return logl(x) + 0.693147180559945309417232121458176568L;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double acoshl(long double x)
+{
+       return acosh(x);
+}
 #endif

diff --git a/src/math/asinhl.c b/src/math/asinhl.c
index e5f3175121cc6145d8d6dc7590d1ddb189065894..8635f52e82381bdaff12badbd7cd170b2e05bf62 100644 (file)
--- a/src/math/asinhl.c
+++ b/src/math/asinhl.c
@@ -32,4 +32,10 @@ long double asinhl(long double x)
        }
        return s ? -x : x;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double asinhl(long double x)
+{
+       return asinh(x);
+}
 #endif

diff --git a/src/math/atanhl.c b/src/math/atanhl.c
index f63d60b153faec6881f5ecc7007e5d2a418adbe1..87cd1cdb5f325f01cf12d19594efbcead93472b1 100644 (file)
--- a/src/math/atanhl.c
+++ b/src/math/atanhl.c
@@ -5,7 +5,7 @@ long double atanhl(long double x)
 {
        return atanh(x);
 }
-#elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
+#elif (LDBL_MANT_DIG == 64 || LDBL_MANT_DIG == 113) && LDBL_MAX_EXP == 16384
 /* atanh(x) = log((1+x)/(1-x))/2 = log1p(2x/(1-x))/2 ~= x + x^3/3 + o(x^5) */
 long double atanhl(long double x)
 {

diff --git a/src/math/coshl.c b/src/math/coshl.c
index 080e5eb0425c0f8bd1ddcf1b0539c2ac9ba081d1..06a56fe3be85d346cd98de913c091e900046a938 100644 (file)
--- a/src/math/coshl.c
+++ b/src/math/coshl.c
@@ -38,4 +38,10 @@ long double coshl(long double x)
        t = expl(0.5*x);
        return 0.5*t*t;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double coshl(long double x)
+{
+       return cosh(x);
+}
 #endif

diff --git a/src/math/erfl.c b/src/math/erfl.c
index 96b74dee6652e827b87d1b4f8357d29ca314b7c3..e267c231059cb85c126b8e86a10d43254ff5aaf0 100644 (file)
--- a/src/math/erfl.c
+++ b/src/math/erfl.c
@@ -340,4 +340,14 @@ long double erfcl(long double x)
        y = 0x1p-16382L;
        return sign ? 2 - y : y*y;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double erfl(long double x)
+{
+       return erf(x);
+}
+long double erfcl(long double x)
+{
+       return erfc(x);
+}
 #endif

diff --git a/src/math/expl.c b/src/math/expl.c
index b62980fa4065f3043f57a76eeb9b4ad9f93c0c28..0a7f44f685ee3b708962dc9683698ffcd8074775 100644 (file)
--- a/src/math/expl.c
+++ b/src/math/expl.c
@@ -119,4 +119,10 @@ long double expl(long double x)
        x = 1.0 + 2.0 * x;
        return scalbnl(x, k);
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double expl(long double x)
+{
+       return exp(x);
+}
 #endif

diff --git a/src/math/expm1l.c b/src/math/expm1l.c
index 21a86c00518783a9f994eaf13aa119214bb2387a..d17150785282f9ce2e566de629b0061cd8848207 100644 (file)
--- a/src/math/expm1l.c
+++ b/src/math/expm1l.c
@@ -114,4 +114,10 @@ long double expm1l(long double x)
        x = px * qx + (px - 1.0);
        return x;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double expm1l(long double x)
+{
+       return expm1(x);
+}
 #endif

diff --git a/src/math/lgammal.c b/src/math/lgammal.c
index 55ec5325ed1382166baec88a264b6fe9940faf03..2b354a7c13245a4dba257a298581a0fa24bac4cf 100644 (file)
--- a/src/math/lgammal.c
+++ b/src/math/lgammal.c
@@ -340,9 +340,16 @@ long double __lgammal_r(long double x, int *sg) {
                r = nadj - r;
        return r;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+double __lgamma_r(double x, int *sg);
+
+long double __lgammal_r(long double x, int *sg)
+{
+       return __lgamma_r(x, sg);
+}
 #endif
 
-#if (LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024) || (LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384)
 extern int __signgam;
 
 long double lgammal(long double x)
@@ -351,4 +358,3 @@ long double lgammal(long double x)
 }
 
 weak_alias(__lgammal_r, lgammal_r);
-#endif

diff --git a/src/math/log10l.c b/src/math/log10l.c
index c7aacf909a45156d1f60513f60ee237ad36e8272..63dcc286d6211d455f291dc9fff2704792a2d236 100644 (file)
--- a/src/math/log10l.c
+++ b/src/math/log10l.c
@@ -182,4 +182,10 @@ done:
        z += e * (L102A);
        return z;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double log10l(long double x)
+{
+       return log10(x);
+}
 #endif

diff --git a/src/math/log1pl.c b/src/math/log1pl.c
index 37da46d23914b311f3ae75ea64ec29fed817fb10..141b5f0b0c295a38d7ab1e1793bdbe2ec35fa2ae 100644 (file)
--- a/src/math/log1pl.c
+++ b/src/math/log1pl.c
@@ -168,4 +168,10 @@ long double log1pl(long double xm1)
        z = z + e * C1;
        return z;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double log1pl(long double x)
+{
+       return log1p(x);
+}
 #endif

diff --git a/src/math/log2l.c b/src/math/log2l.c
index d00531d5976e26796ab126d616fb3057ce280029..722b451a026ae2ccd1b70381c4a87b7936ddb3c2 100644 (file)
--- a/src/math/log2l.c
+++ b/src/math/log2l.c
@@ -173,4 +173,10 @@ done:
        z += e;
        return z;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double log2l(long double x)
+{
+       return log2(x);
+}
 #endif

diff --git a/src/math/logl.c b/src/math/logl.c
index 03c5188fd9d79389d938cff1b6900ea90fa4aaf3..5d5365929e62a351c4544bd9e6aace6e02834292 100644 (file)
--- a/src/math/logl.c
+++ b/src/math/logl.c
@@ -166,4 +166,10 @@ long double logl(long double x)
        z = z + e * C1; /* This sum has an error of 1/2 lsb. */
        return z;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double logl(long double x)
+{
+       return log(x);
+}
 #endif

diff --git a/src/math/powl.c b/src/math/powl.c
index ce6274cf73e00c5233acd4954e3d7d2dccf6ff70..a765706d9b34c4121a5d62750f158681ce39be58 100644 (file)
--- a/src/math/powl.c
+++ b/src/math/powl.c
@@ -513,5 +513,10 @@ static long double powil(long double x, int nn)
                y = 1.0/y;
        return y;
 }
-
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double powl(long double x, long double y)
+{
+       return pow(x, y);
+}
 #endif

diff --git a/src/math/sinhl.c b/src/math/sinhl.c
index 4864ddfa6a77c215a2dc1158ca1bb75137730063..b305d4d2f3e9c5965bdbf4ddb445adef0168fd85 100644 (file)
--- a/src/math/sinhl.c
+++ b/src/math/sinhl.c
@@ -34,4 +34,10 @@ long double sinhl(long double x)
        t = expl(0.5*absx);
        return h*t*t;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double sinhl(long double x)
+{
+       return sinh(x);
+}
 #endif

diff --git a/src/math/tanhl.c b/src/math/tanhl.c
index f594b85f3d6ce3a6722ebcf49eb6584410a41ad0..4e1aa9f87dde30bc14f43bed200fe2a968e45108 100644 (file)
--- a/src/math/tanhl.c
+++ b/src/math/tanhl.c
@@ -39,4 +39,10 @@ long double tanhl(long double x)
        }
        return sign ? -t : t;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double tanhl(long double x)
+{
+       return tanh(x);
+}
 #endif

diff --git a/src/math/tgammal.c b/src/math/tgammal.c
index 5c1a02a6399692ca60037f346e4e1be1f8163f19..5336c5b194e86c57797562191ad39cd4af7fafe4 100644 (file)
--- a/src/math/tgammal.c
+++ b/src/math/tgammal.c
@@ -272,4 +272,10 @@ small:
                q = z / (x * __polevll(x, S, 8));
        return q;
 }
+#elif LDBL_MANT_DIG == 113 && LDBL_MAX_EXP == 16384
+// TODO: broken implementation to make things compile
+long double tgammal(long double x)
+{
+       return tgamma(x);
+}
 #endif