Subversion Repositories or1k

 
/* @(#)z_isnanf.c 1.0 98/08/13 */
/******************************************************************
 * isnanf
 *
 * Input:
 *   x - pointer to a floating point value
 *
 * Output:
 *   An integer that indicates if the number is NaN.
 *
 * Description:
 *   This routine returns an integer that indicates if the number
 *   passed in is NaN (1) or is finite (0).
 *
 *****************************************************************/
 
#include "fdlibm.h"
#include "zmath.h"
 
int
_DEFUN (isnanf, (float),
        float x)
{
  __int32_t wx;
  int exp;
 
  GET_FLOAT_WORD (wx, x);
  exp = (wx & 0x7f800000) >> 23;
 
  if ((exp == 0x7f8) && (wx & 0x7fffff))
    return (1);
  else
    return (0);
}
 
 
#ifdef _DOUBLE_IS_32BITS
 
int
_DEFUN (isnan, (double),
        double x)
{
  return isnanf((float) x);
}
 
#endif /* defined(_DOUBLE_IS_32BITS) */