X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Facoshl.c;h=4aa84acb2ea3684b0ce4c3c2199b0ab6a6d4814e;hb=a2a328bd896fa909542b492601e3689119b95de0;hp=472c71cbc81749cd22c1886fd02ebfd77b5e628f;hpb=969ddbc423238291d5c7982790bbe72720627ba4;p=musl

diff --git a/src/math/acoshl.c b/src/math/acoshl.c
index 472c71cb..4aa84acb 100644
--- a/src/math/acoshl.c
+++ b/src/math/acoshl.c
@@ -9,16 +9,14 @@ 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;
 }