isnan() returns unexpected result for std::complex

Hi,

consider this minimal example( tested with gcc 9 and clang 9):

#include <blaze/math/DynamicVector.h>
#include <iostream>
#include <limits>

int main(void) {
    blaze::DynamicVector<double> v(1);
    v[0] = std::numeric_limits<double>::quiet_NaN();
    // 1
    std::cout << blaze::isnan(v) << std::endl;

    blaze::DynamicVector<std::complex<double>> w(1);
    w[0] = std::complex<double>(std::numeric_limits<double>::quiet_NaN(),
                                std::numeric_limits<double>::quiet_NaN());
    // 0
    std::cout << blaze::isnan(w) << std::endl;

    return 0;
}

Is isnan() expected to work with std::complex<double>? I believe according to the documentation it should work.

Tim

Comments (4)