/*

 * 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>

#include <ldbl-classify-compat.h>

#include <shlib-compat.h>

int

__isinf (double x)

{

  int32_t hx, lx;

  EXTRACT_WORDS (hx, lx, x);

  lx |= (hx & 0x7fffffff) ^ 0x7ff00000;

  lx |= -lx;

  return ~(lx >> 31) & (hx >> 30);

}

hidden_def (__isinf)

weak_alias (__isinf, isinf)

#ifdef NO_LONG_DOUBLE

# if LDBL_CLASSIFY_COMPAT && SHLIB_COMPAT (libc, GLIBC_2_0, GLIBC_2_23)

compat_symbol (libc, __isinf, __isinfl, GLIBC_2_0);

# endif

weak_alias (__isinf, isinfl)

#endif