s_isnanf.c revision 1.3 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.3 jtc #if defined(LIBM_SCCS) && !defined(lint)
17 1.3 jtc static char rcsid[] = "$Id: s_isnanf.c,v 1.3 1994/09/22 16:44:47 jtc Exp $";
18 1.1 jtc #endif
19 1.1 jtc
20 1.1 jtc /*
21 1.1 jtc * isnanf(x) returns 1 is x is nan, else 0;
22 1.1 jtc * no branching!
23 1.1 jtc */
24 1.1 jtc
25 1.1 jtc #include "math.h"
26 1.1 jtc #include "math_private.h"
27 1.1 jtc
28 1.1 jtc #ifdef __STDC__
29 1.1 jtc int isnanf(float x)
30 1.1 jtc #else
31 1.1 jtc int isnanf(x)
32 1.1 jtc float x;
33 1.1 jtc #endif
34 1.1 jtc {
35 1.2 jtc int32_t ix;
36 1.1 jtc GET_FLOAT_WORD(ix,x);
37 1.1 jtc ix &= 0x7fffffff;
38 1.1 jtc ix = 0x7f800000 - ix;
39 1.2 jtc return (int)(((u_int32_t)(ix))>>31);
40 1.1 jtc }
41