X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2F__fpclassifyl.c;h=481c0b949974ad43a59e95846bf9fb83ff5cbf23;hb=0c8bc102f287d3993751d80ba2dffb01e0c8bc7f;hp=6365c58854c8011327f69f55bda30a6adcceb376;hpb=af5f6d9556441487e5c66a7a4cfeddf4ed354aa7;p=musl

diff --git a/src/math/__fpclassifyl.c b/src/math/__fpclassifyl.c
index 6365c588..481c0b94 100644
--- a/src/math/__fpclassifyl.c
+++ b/src/math/__fpclassifyl.c
@@ -1,15 +1,20 @@
 #include "libm.h"
 
 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
-
+int __fpclassifyl(long double x)
+{
+	return __fpclassify(x);
+}
 #elif LDBL_MANT_DIG == 64 && LDBL_MAX_EXP == 16384
-/* invalid representations (top bit of u.i.m is wrong) are not handled */
 int __fpclassifyl(long double x)
 {
 	union ldshape u = {x};
 	int e = u.i.se & 0x7fff;
-	if (!e)
+	int msb = u.i.m>>63;
+	if (!e && !msb)
 		return u.i.m ? FP_SUBNORMAL : FP_ZERO;
+	if (!msb)
+		return FP_NAN;
 	if (e == 0x7fff)
 		return u.i.m << 1 ? FP_NAN : FP_INFINITE;
 	return FP_NORMAL;
@@ -19,12 +24,11 @@ int __fpclassifyl(long double x)
 {
 	union ldshape u = {x};
 	int e = u.i.se & 0x7fff;
+	u.i.se = 0;
 	if (!e)
 		return u.i2.lo | u.i2.hi ? FP_SUBNORMAL : FP_ZERO;
-	if (e == 0x7fff) {
-		u.i.se = 0;
+	if (e == 0x7fff)
 		return u.i2.lo | u.i2.hi ? FP_NAN : FP_INFINITE;
-	}
 	return FP_NORMAL;
 }
 #endif