>From 53c656d9cb03e53182ff01004d1f2bf18854635a Mon Sep 17 00:00:00 2001
From: Szabolcs Nagy <nsz@port70.net>
Date: Tue, 10 Mar 2015 20:01:20 +0000
Subject: [PATCH 3/4] 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.
---
 src/internal/libm.h |   14 ++++++++++++++
 src/math/acoshl.c   |    6 ++++++
 src/math/asinhl.c   |    6 ++++++
 src/math/atanhl.c   |    2 +-
 src/math/coshl.c    |    6 ++++++
 src/math/erfl.c     |   10 ++++++++++
 src/math/expl.c     |    6 ++++++
 src/math/expm1l.c   |    6 ++++++
 src/math/lgammal.c  |   10 ++++++++--
 src/math/log10l.c   |    6 ++++++
 src/math/log1pl.c   |    6 ++++++
 src/math/log2l.c    |    6 ++++++
 src/math/logl.c     |    6 ++++++
 src/math/powl.c     |    7 ++++++-
 src/math/sinhl.c    |    6 ++++++
 src/math/tanhl.c    |    6 ++++++
 src/math/tgammal.c  |    6 ++++++
 17 files changed, 111 insertions(+), 4 deletions(-)

diff --git a/src/internal/libm.h b/src/internal/libm.h
index 88a7eb4..df86411 100644
--- 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 4aa84ac..8d4b43f 100644
--- 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 e5f3175..8635f52 100644
--- 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 f63d60b..87cd1cd 100644
--- 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 080e5eb..06a56fe 100644
--- 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 96b74de..e267c23 100644
--- 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 b62980f..0a7f44f 100644
--- 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 21a86c0..d171507 100644
--- 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 55ec532..2b354a7 100644
--- 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 c7aacf9..63dcc28 100644
--- 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 37da46d..141b5f0 100644
--- 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 d00531d..722b451 100644
--- 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 03c5188..5d53659 100644
--- 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 ce6274c..a765706 100644
--- 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 4864ddf..b305d4d 100644
--- 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 f594b85..4e1aa9f 100644
--- 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 5c1a02a..5336c5b 100644
--- 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
-- 
1.7.9.5