x86 long double: Support pseudo numbers in isnanl

This syncs up isnanl behaviour with gcc.  Also move the isnanl
implementation to sysdeps/x86 and remove the sysdeps/x86_64 version.

Reviewed-by: Adhemerval Zanella  <adhemerval.zanella@linaro.org>

diff --git a/sysdeps/i386/fpu/s_isnanl.c b/sysdeps/x86/fpu/s_isnanl.c
similarity index 76%
rename from sysdeps/i386/fpu/s_isnanl.c
rename to sysdeps/x86/fpu/s_isnanl.c
index fb97317bc997b70a22f18ffd75af39b8ea2dbbed..45541df08acfb3ff9d1d9b3205ebb37ad250e19a 100644 (file)
--- a/sysdeps/i386/fpu/s_isnanl.c
+++ b/sysdeps/x86/fpu/s_isnanl.c
@@ -28,16 +28,18 @@ static char rcsid[] = "$NetBSD: $";
 
 int __isnanl(long double x)
 {
-       int32_t se,hx,lx;
+       int32_t se,hx,lx,pn;
        GET_LDOUBLE_WORDS(se,hx,lx,x);
        se = (se & 0x7fff) << 1;
-       /* The additional & 0x7fffffff is required because Intel's
-          extended format has the normally implicit 1 explicit
-          present.  Sigh!  */
+       /* Detect pseudo-normal numbers, i.e. exponent is non-zero and the top
+          bit of the significand is not set.   */
+       pn = (uint32_t)((~hx & 0x80000000) & (se|(-se)))>>31;
+       /* Clear the significand bit when computing mantissa.  */
        lx |= hx & 0x7fffffff;
        se |= (uint32_t)(lx|(-lx))>>31;
        se = 0xfffe - se;
-       return (int)((uint32_t)(se))>>16;
+
+       return (int)(((uint32_t)(se)) >> 16) | pn;
 }
 hidden_def (__isnanl)
 weak_alias (__isnanl, isnanl)

diff --git a/sysdeps/x86_64/fpu/s_isnanl.c b/sysdeps/x86_64/fpu/s_isnanl.c
deleted file mode 100644 (file)
index 06e69c3..0000000
--- a/sysdeps/x86_64/fpu/s_isnanl.c
+++ /dev/null
@@ -1 +0,0 @@
-#include <sysdeps/i386/fpu/s_isnanl.c>