X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;ds=sidebyside;f=src%2Fmath%2Facoshl.c;h=8d4b43f64cdd899877bf39f43c7657f1f3b3f17c;hb=9b132e556774c744f9052581d2d8d0fab417e97c;hp=472c71cbc81749cd22c1886fd02ebfd77b5e628f;hpb=482ccd2f7497a79ca83e998f54e823e7cedaaa6e;p=musl

diff --git a/src/math/acoshl.c b/src/math/acoshl.c
index 472c71cb..8d4b43f6 100644
--- a/src/math/acoshl.c
+++ b/src/math/acoshl.c
@@ -9,17 +9,21 @@ long double acoshl(long double x)
 /* acosh(x) = log(x + sqrt(x*x-1)) */
 long double acoshl(long double x)
 {
-	union {
-		long double f;
-		struct{uint64_t m; int16_t se; uint16_t pad;} i;
-	} u = {.f = x};
+	union ldshape u = {x};
+	int e = u.i.se & 0x7fff;
 
-	if (u.i.se < 0x3fff + 1)
-		/* x < 2, invalid if x < 1 or nan */
+	if (e < 0x3fff + 1)
+		/* |x| < 2, invalid if x < 1 or nan */
 		return log1pl(x-1 + sqrtl((x-1)*(x-1)+2*(x-1)));
-	if (u.i.se < 0x3fff + 32)
-		/* x < 0x1p32 */
+	if (e < 0x3fff + 32)
+		/* |x| < 0x1p32 */
 		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