37 lines
731 B
C
37 lines
731 B
C
/*
|
|
* 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
|
|
{
|
|
int32_t hx,lx;
|
|
EXTRACT_WORDS(hx,lx,x);
|
|
lx |= (hx & 0x7fffffff) ^ 0x7ff00000;
|
|
lx |= -lx;
|
|
return ~(lx >> 31) & (hx >> 30);
|
|
}
|
|
weak_alias (__isinf, isinf)
|
|
#ifdef NO_LONG_DOUBLE
|
|
strong_alias (__isinf, __isinfl)
|
|
weak_alias (__isinf, isinfl)
|
|
#endif
|