s_isnanf.c revision 1.5 1 1.1 jtc /* s_isnanf.c -- float version of s_isnan.c.
2 1.1 jtc * Conversion to float by Ian Lance Taylor, Cygnus Support, ian (at) cygnus.com.
3 1.1 jtc */
4 1.1 jtc
5 1.1 jtc /*
6 1.1 jtc * ====================================================
7 1.1 jtc * Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
8 1.1 jtc *
9 1.1 jtc * Developed at SunPro, a Sun Microsystems, Inc. business.
10 1.1 jtc * Permission to use, copy, modify, and distribute this
11 1.1 jtc * software is freely granted, provided that this notice
12 1.1 jtc * is preserved.
13 1.1 jtc * ====================================================
14 1.1 jtc */
15 1.1 jtc
16 1.5 lukem #include <sys/cdefs.h>
17 1.3 jtc #if defined(LIBM_SCCS) && !defined(lint)
18 1.5 lukem __RCSID("$NetBSD: s_isnanf.c,v 1.5 1997/10/09 11:32:30 lukem Exp $");
19 1.1 jtc #endif
20 1.1 jtc
21 1.1 jtc /*
22 1.1 jtc * isnanf(x) returns 1 is x is nan, else 0;
23 1.1 jtc * no branching!
24 1.1 jtc */
25 1.1 jtc
26 1.1 jtc #include "math.h"
27 1.1 jtc #include "math_private.h"
28 1.1 jtc
29 1.1 jtc #ifdef __STDC__
30 1.1 jtc int isnanf(float x)
31 1.1 jtc #else
32 1.1 jtc int isnanf(x)
33 1.1 jtc float x;
34 1.1 jtc #endif
35 1.1 jtc {
36 1.2 jtc int32_t ix;
37 1.1 jtc GET_FLOAT_WORD(ix,x);
38 1.1 jtc ix &= 0x7fffffff;
39 1.1 jtc ix = 0x7f800000 - ix;
40 1.2 jtc return (int)(((u_int32_t)(ix))>>31);
41 1.1 jtc }
42