X-Git-Url: http://nsz.repo.hu/git/?a=blobdiff_plain;f=src%2Fmath%2Fatanl.c;h=d29e6316cf9e8b4b91114d20f90f9e3e8327fb80;hb=d1a2ead878c27ac4ec600740320f8b76e1f961e9;hp=e76693e47d339baf503a08832898ff9f34d431d9;hpb=c6383b7b10303457306932584fc23f24b5153a81;p=musl

diff --git a/src/math/atanl.c b/src/math/atanl.c
index e76693e4..d29e6316 100644
--- a/src/math/atanl.c
+++ b/src/math/atanl.c
@@ -70,8 +70,8 @@ long double atanl(long double x)
 	union IEEEl2bits u;
 	long double w,s1,s2,z;
 	int id;
-	int16_t expsign, expt;
-	int32_t expman;
+	uint16_t expsign, expt;
+	uint32_t expman;
 
 	u.e = x;
 	expsign = u.xbits.expsign;
@@ -81,15 +81,16 @@ long double atanl(long double x)
 		    ((u.bits.manh&~LDBL_NBIT)|u.bits.manl)!=0)  /* NaN */
 			return x+x;
 		z = atanhi[3] + 0x1p-120f;
-		return expsign < 0 ? -z : z;
+		return expsign>>15 ? -z : z;
 	}
 	/* Extract the exponent and the first few bits of the mantissa. */
 	/* XXX There should be a more convenient way to do this. */
 	expman = (expt << 8) | ((u.bits.manh >> (LDBL_MANH_SIZE - 9)) & 0xff);
 	if (expman < ((0x3fff - 2) << 8) + 0xc0) {  /* |x| < 0.4375 */
 		if (expt < 0x3fff - 32) {   /* if |x| is small, atanl(x)~=x */
-			/* raise inexact if x!=0 */
-			FORCE_EVAL(x + 0x1p120f);
+			/* raise underflow if subnormal */
+			if (expt == 0)
+				FORCE_EVAL((float)x);
 			return x;
 		}
 		id = -1;
@@ -122,6 +123,6 @@ long double atanl(long double x)
 	if (id < 0)
 		return x - x*(s1+s2);
 	z = atanhi[id] - ((x*(s1+s2) - atanlo[id]) - x);
-	return expsign < 0 ? -z : z;
+	return expsign>>15 ? -z : z;
 }
 #endif