Как проверить double==NaN?

Использую NaN для разметки области определения дырявой функции на плоскости. Думал, что операция типа

double f = numeric_limits<double>::quiet_NaN();
bool test = (0<f && f<=0);

мне поможет, но компилятор посчитал меня идиотом и оптимизировал test=false;
Как на C++ проверить double==NaN?

Воспользуйся функцией _isnan из библиотеки float.h:
int _isnan( double x );
_isnan returns a nonzero value (TRUE) if the argument x is a NaN; otherwise it returns 0 (FALSE).

[Drakosha]

или

#include <math.h>

int isinf(double value);

int isnan(double value);

int finite(double value);

Воспользовался, работает , спасибо . Кстати, заметил, что с битовой маской вместо NaN всё быстрее получается, если NaN-ов много. Он что, прерывание (или исключение вызывает)?