home *** CD-ROM | disk | FTP | other *** search
- /*
- * Written by J.T. Conklin <jtc@netbsd.org>.
- * Changed to return -1 for -Inf by Ulrich Drepper <drepper@cygnus.com>.
- * Public domain.
- */
-
- #if defined(LIBM_SCCS) && !defined(lint)
- static char rcsid[] = "$NetBSD: s_isinf.c,v 1.3 1995/05/11 23:20:14 jtc Exp $";
- #endif
-
- /*
- * isinf(x) returns 1 is x is inf, -1 if x is -inf, else 0;
- * no branching!
- */
-
- #include "math.h"
- #include "math_private.h"
-
- #ifdef __STDC__
- int __isinf(double x)
- #else
- int __isinf(x)
- double x;
- #endif
- {
- u_int32_t hx;
- int32_t lx;
- EXTRACT_WORDS(hx,lx,x);
- lx |= (hx & 0x7fffffff) ^ 0x7ff00000;
- lx |= -lx;
- return ~(lx >> 31) & (1 - ((hx >> 30) & 2));
- }
- weak_alias (__isinf, isinf)
- #ifdef NO_LONG_DOUBLE
- strong_alias (__isinf, __isinfl)
- weak_alias (__isinf, isinfl)
- #endif
-