int willusmath_is_nan(double x) {
  int isz, isgt, islt;
  isz = (x == 0.);
  isgt = (x > 0.);
  islt = (x < 0.);
  return ((isz && isgt) || (isz && islt) || (isgt && islt) ||
          (!isz && !isgt && !islt));
}


x86 assembly (clang, -O0):
<willusmath_is_nan>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
sete   %al
setnp  %cl
and    %cl,%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0xc(%rbp)
movsd  -0x8(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x10(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x8(%rbp),%xmm0
seta   %al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0xc(%rbp)
je     1168 <willusmath_is_nan+0x68>
mov    $0x1,%al
cmpl   $0x0,-0x10(%rbp)
mov    %al,-0x15(%rbp)
jne    11ca <willusmath_is_nan+0xca>
cmpl   $0x0,-0xc(%rbp)
je     1181 <willusmath_is_nan+0x81>
mov    $0x1,%al
cmpl   $0x0,-0x14(%rbp)
mov    %al,-0x15(%rbp)
jne    11ca <willusmath_is_nan+0xca>
cmpl   $0x0,-0x10(%rbp)
je     119a <willusmath_is_nan+0x9a>
mov    $0x1,%al
cmpl   $0x0,-0x14(%rbp)
mov    %al,-0x15(%rbp)
jne    11ca <willusmath_is_nan+0xca>
xor    %eax,%eax
cmpl   $0x0,-0xc(%rbp)
mov    %al,-0x16(%rbp)
jne    11c4 <willusmath_is_nan+0xc4>
xor    %eax,%eax
cmpl   $0x0,-0x10(%rbp)
mov    %al,-0x16(%rbp)
jne    11c4 <willusmath_is_nan+0xc4>
cmpl   $0x0,-0x14(%rbp)
setne  %al
xor    $0xff,%al
mov    %al,-0x16(%rbp)
mov    -0x16(%rbp),%al
mov    %al,-0x15(%rbp)
mov    -0x15(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<willusmath_is_nan>:
xor    %eax,%eax
ucomisd %xmm0,%xmm0
setp   %al
ret

LLVM IR (-O0):
define i32 @willusmath_is_nan(double %0) {
%2 = alloca double
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
store double %0, double* %2
%6 = load double, double* %2
%7 = fcmp oeq double %6, 0.000000e+00
%8 = zext i1 %7 to i32
store i32 %8, i32* %3
%9 = load double, double* %2
%10 = fcmp ogt double %9, 0.000000e+00
%11 = zext i1 %10 to i32
store i32 %11, i32* %4
%12 = load double, double* %2
%13 = fcmp olt double %12, 0.000000e+00
%14 = zext i1 %13 to i32
store i32 %14, i32* %5
%15 = load i32, i32* %3
%16 = icmp ne i32 %15, 0
br i1 %16, label %17, label %20
17: ; preds = %1
%18 = load i32, i32* %4
%19 = icmp ne i32 %18, 0
br i1 %19, label %44, label %20
20: ; preds = %17, %1
%21 = load i32, i32* %3
%22 = icmp ne i32 %21, 0
br i1 %22, label %23, label %26
23: ; preds = %20
%24 = load i32, i32* %5
%25 = icmp ne i32 %24, 0
br i1 %25, label %44, label %26
26: ; preds = %23, %20
%27 = load i32, i32* %4
%28 = icmp ne i32 %27, 0
br i1 %28, label %29, label %32
29: ; preds = %26
%30 = load i32, i32* %5
%31 = icmp ne i32 %30, 0
br i1 %31, label %44, label %32
32: ; preds = %29, %26
%33 = load i32, i32* %3
%34 = icmp ne i32 %33, 0
br i1 %34, label %42, label %35
35: ; preds = %32
%36 = load i32, i32* %4
%37 = icmp ne i32 %36, 0
br i1 %37, label %42, label %38
38: ; preds = %35
%39 = load i32, i32* %5
%40 = icmp ne i32 %39, 0
%41 = xor i1 %40, true
br label %42
42: ; preds = %38, %35, %32
%43 = phi i1 [ false, %35 ], [ false, %32 ], [ %41, %38 ]
br label %44
44: ; preds = %42, %29, %23, %17
%45 = phi i1 [ true, %29 ], [ true, %23 ], [ true, %17 ], [ %43, %42 ]
%46 = zext i1 %45 to i32
ret i32 %46
}

LLVM IR (-O2):
define i32 @willusmath_is_nan(double %0) {
%2 = fcmp uno double %0, 0.000000e+00
%3 = zext i1 %2 to i32
ret i32 %3
}