elemental function sisnan(x) result(res)
real(kind(1.e0)), intent(in) :: x

short black skirt, taurus,

good choice!,

logical :: res
res = ieor(ibclr(transfer(x,sNan),SPSB), sNaN) == 0
end function