Chris@366: 
Chris@366: #ifndef NAN_INF_H
Chris@366: #define NAN_INF_H
Chris@366: 
Chris@366: #define ISNAN(x) (sizeof(x) == sizeof(double) ? ISNANd(x) : ISNANf(x))
Chris@366: static inline int ISNANf(float x) { return x != x; }
Chris@366: static inline int ISNANd(double x) { return x != x; }
Chris@366:           
Chris@366: #define ISINF(x) (sizeof(x) == sizeof(double) ? ISINFd(x) : ISINFf(x))
Chris@366: static inline int ISINFf(float x) { return !ISNANf(x) && ISNANf(x - x); }
Chris@366: static inline int ISINFd(double x) { return !ISNANd(x) && ISNANd(x - x); }
Chris@366: 
Chris@366: #endif