sysdeps/libm-i387/s_isnanl.c

   1 /* s_isnanl.c -- long double version for i387 of s_isnan.c.
   2  * Conversion to long double by Ulrich Drepper,
   3  * Cygnus Support, drepper@cygnus.com.
   4  */
   5
   6 /*
   7  * ====================================================
   8  * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
   9  *
  10  * Developed at SunPro, a Sun Microsystems, Inc. business.
  11  * Permission to use, copy, modify, and distribute this
  12  * software is freely granted, provided that this notice
  13  * is preserved.
  14  * ====================================================
  15  */
  16
  17 #if defined(LIBM_SCCS) && !defined(lint)
  18 static char rcsid[] = "$NetBSD: $";
  19 #endif
  20
  21 /*
  22  * isnanl(x) returns 1 is x is nan, else 0;
  23  * no branching!
  24  */
  25
  26 #include "math.h"
  27 #include "math_private.h"
  28
  29 #ifdef __STDC__
  30         int __isnanl(long double x)
  31 #else
  32         int __isnanl(x)
  33         long double x;
  34 #endif
  35 {
  36         int32_t se,hx,lx;
  37         GET_LDOUBLE_WORDS(se,hx,lx,x);
  38         se = (se & 0x7fff) << 1;
  39         /* The additional & 0x7fffffff is required because Intel's
  40            extended format has the normally implicit 1 explicit
  41            present.  Sigh!  */
  42         lx |= hx & 0x7fffffff;
  43         se |= (u_int32_t)(lx|(-lx))>>31;
  44         se = 0xfffe - se;
  45         return (int)((u_int32_t)(se))>>16;
  46 }
  47 weak_alias (__isnanl, isnanl)