diff options
author | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2015-08-16 14:00:02 +0200 |
---|---|---|
committer | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2015-08-16 14:00:02 +0200 |
commit | d6a4805fdf6ef4cb13c05c218006cda32a8c857a (patch) | |
tree | 9490d33c1500e40c18d95faadff74e3b77b2b917 /test | |
parent | a40f6ab276a476113defe883eec00c617ff22920 (diff) |
Protect further isnan/isfinite/isinf calls
Diffstat (limited to 'test')
-rw-r--r-- | test/packetmath.cpp | 16 | ||||
-rw-r--r-- | test/stable_norm.cpp | 50 |
2 files changed, 33 insertions, 33 deletions
diff --git a/test/packetmath.cpp b/test/packetmath.cpp index 3cf82eae0..1cc7a9cd9 100644 --- a/test/packetmath.cpp +++ b/test/packetmath.cpp @@ -318,7 +318,7 @@ template<typename Scalar> void packetmath_real() data1[1] = std::numeric_limits<Scalar>::epsilon(); packet_helper<internal::packet_traits<Scalar>::HasExp,Packet> h; h.store(data2, internal::pexp(h.load(data1))); - VERIFY(numext::isnan(data2[0])); + VERIFY((numext::isnan)(data2[0])); VERIFY_IS_EQUAL(std::exp(std::numeric_limits<Scalar>::epsilon()), data2[1]); data1[0] = -std::numeric_limits<Scalar>::epsilon(); @@ -354,34 +354,34 @@ template<typename Scalar> void packetmath_real() data1[1] = std::numeric_limits<Scalar>::epsilon(); packet_helper<internal::packet_traits<Scalar>::HasLog,Packet> h; h.store(data2, internal::plog(h.load(data1))); - VERIFY(std::isnan(data2[0])); + VERIFY((numext::isnan)(data2[0])); // VERIFY_IS_EQUAL(std::log(std::numeric_limits<Scalar>::epsilon()), data2[1]); data1[0] = -std::numeric_limits<Scalar>::epsilon(); data1[1] = 0; h.store(data2, internal::plog(h.load(data1))); - VERIFY(std::isnan(data2[0])); + VERIFY((numext::isnan)(data2[0])); // VERIFY_IS_EQUAL(std::log(0), data2[1]); data1[0] = (std::numeric_limits<Scalar>::min)(); data1[1] = -(std::numeric_limits<Scalar>::min)(); h.store(data2, internal::plog(h.load(data1))); VERIFY_IS_EQUAL(std::log((std::numeric_limits<Scalar>::min)()), data2[0]); - // VERIFY(std::isnan(data2[1])); + // VERIFY((numext::isnan)(data2[1])); data1[0] = std::numeric_limits<Scalar>::denorm_min(); data1[1] = -std::numeric_limits<Scalar>::denorm_min(); h.store(data2, internal::plog(h.load(data1))); // VERIFY_IS_EQUAL(std::log(std::numeric_limits<Scalar>::denorm_min()), data2[0]); - // VERIFY(std::isnan(data2[1])); + // VERIFY((numext::isnan)(data2[1])); data1[0] = -1.0f; h.store(data2, internal::plog(h.load(data1))); - VERIFY(std::isnan(data2[0])); + VERIFY((numext::isnan)(data2[0])); #if !EIGEN_FAST_MATH h.store(data2, internal::psqrt(h.load(data1))); - VERIFY(numext::isnan(data2[0])); - VERIFY(numext::isnan(data2[1])); + VERIFY((numext::isnan)(data2[0])); + VERIFY((numext::isnan)(data2[1])); #endif } } diff --git a/test/stable_norm.cpp b/test/stable_norm.cpp index 99cb0d5a6..7561ae8be 100644 --- a/test/stable_norm.cpp +++ b/test/stable_norm.cpp @@ -40,7 +40,7 @@ template<typename MatrixType> void stable_norm(const MatrixType& m) && "the stable norm algorithm cannot be guaranteed on this computer"); Scalar inf = std::numeric_limits<RealScalar>::infinity(); - if(NumTraits<Scalar>::IsComplex && numext::isnan(inf*RealScalar(1)) ) + if(NumTraits<Scalar>::IsComplex && (numext::isnan)(inf*RealScalar(1)) ) { complex_real_product_ok = false; static bool first = true; @@ -81,18 +81,18 @@ template<typename MatrixType> void stable_norm(const MatrixType& m) RealScalar size = static_cast<RealScalar>(m.size()); // test numext::isfinite - VERIFY(!numext::isfinite( std::numeric_limits<RealScalar>::infinity())); - VERIFY(!numext::isfinite(sqrt(-abs(big)))); + VERIFY(!(numext::isfinite)( std::numeric_limits<RealScalar>::infinity())); + VERIFY(!(numext::isfinite)(sqrt(-abs(big)))); // test overflow - VERIFY(numext::isfinite(sqrt(size)*abs(big))); + VERIFY((numext::isfinite)(sqrt(size)*abs(big))); VERIFY_IS_NOT_APPROX(sqrt(copy(vbig.squaredNorm())), abs(sqrt(size)*big)); // here the default norm must fail VERIFY_IS_APPROX(vbig.stableNorm(), sqrt(size)*abs(big)); VERIFY_IS_APPROX(vbig.blueNorm(), sqrt(size)*abs(big)); VERIFY_IS_APPROX(vbig.hypotNorm(), sqrt(size)*abs(big)); // test underflow - VERIFY(numext::isfinite(sqrt(size)*abs(small))); + VERIFY((numext::isfinite)(sqrt(size)*abs(small))); VERIFY_IS_NOT_APPROX(sqrt(copy(vsmall.squaredNorm())), abs(sqrt(size)*small)); // here the default norm must fail VERIFY_IS_APPROX(vsmall.stableNorm(), sqrt(size)*abs(small)); VERIFY_IS_APPROX(vsmall.blueNorm(), sqrt(size)*abs(small)); @@ -115,39 +115,39 @@ template<typename MatrixType> void stable_norm(const MatrixType& m) { v = vrand; v(i,j) = std::numeric_limits<RealScalar>::quiet_NaN(); - VERIFY(!numext::isfinite(v.squaredNorm())); VERIFY(numext::isnan(v.squaredNorm())); - VERIFY(!numext::isfinite(v.norm())); VERIFY(numext::isnan(v.norm())); - VERIFY(!numext::isfinite(v.stableNorm())); VERIFY(numext::isnan(v.stableNorm())); - VERIFY(!numext::isfinite(v.blueNorm())); VERIFY(numext::isnan(v.blueNorm())); - VERIFY(!numext::isfinite(v.hypotNorm())); VERIFY(numext::isnan(v.hypotNorm())); + VERIFY(!(numext::isfinite)(v.squaredNorm())); VERIFY((numext::isnan)(v.squaredNorm())); + VERIFY(!(numext::isfinite)(v.norm())); VERIFY((numext::isnan)(v.norm())); + VERIFY(!(numext::isfinite)(v.stableNorm())); VERIFY((numext::isnan)(v.stableNorm())); + VERIFY(!(numext::isfinite)(v.blueNorm())); VERIFY((numext::isnan)(v.blueNorm())); + VERIFY(!(numext::isfinite)(v.hypotNorm())); VERIFY((numext::isnan)(v.hypotNorm())); } // +inf { v = vrand; v(i,j) = std::numeric_limits<RealScalar>::infinity(); - VERIFY(!numext::isfinite(v.squaredNorm())); VERIFY(isPlusInf(v.squaredNorm())); - VERIFY(!numext::isfinite(v.norm())); VERIFY(isPlusInf(v.norm())); - VERIFY(!numext::isfinite(v.stableNorm())); + VERIFY(!(numext::isfinite)(v.squaredNorm())); VERIFY(isPlusInf(v.squaredNorm())); + VERIFY(!(numext::isfinite)(v.norm())); VERIFY(isPlusInf(v.norm())); + VERIFY(!(numext::isfinite)(v.stableNorm())); if(complex_real_product_ok){ VERIFY(isPlusInf(v.stableNorm())); } - VERIFY(!numext::isfinite(v.blueNorm())); VERIFY(isPlusInf(v.blueNorm())); - VERIFY(!numext::isfinite(v.hypotNorm())); VERIFY(isPlusInf(v.hypotNorm())); + VERIFY(!(numext::isfinite)(v.blueNorm())); VERIFY(isPlusInf(v.blueNorm())); + VERIFY(!(numext::isfinite)(v.hypotNorm())); VERIFY(isPlusInf(v.hypotNorm())); } // -inf { v = vrand; v(i,j) = -std::numeric_limits<RealScalar>::infinity(); - VERIFY(!numext::isfinite(v.squaredNorm())); VERIFY(isPlusInf(v.squaredNorm())); - VERIFY(!numext::isfinite(v.norm())); VERIFY(isPlusInf(v.norm())); - VERIFY(!numext::isfinite(v.stableNorm())); + VERIFY(!(numext::isfinite)(v.squaredNorm())); VERIFY(isPlusInf(v.squaredNorm())); + VERIFY(!(numext::isfinite)(v.norm())); VERIFY(isPlusInf(v.norm())); + VERIFY(!(numext::isfinite)(v.stableNorm())); if(complex_real_product_ok) { VERIFY(isPlusInf(v.stableNorm())); } - VERIFY(!numext::isfinite(v.blueNorm())); VERIFY(isPlusInf(v.blueNorm())); - VERIFY(!numext::isfinite(v.hypotNorm())); VERIFY(isPlusInf(v.hypotNorm())); + VERIFY(!(numext::isfinite)(v.blueNorm())); VERIFY(isPlusInf(v.blueNorm())); + VERIFY(!(numext::isfinite)(v.hypotNorm())); VERIFY(isPlusInf(v.hypotNorm())); } // mix @@ -157,11 +157,11 @@ template<typename MatrixType> void stable_norm(const MatrixType& m) v = vrand; v(i,j) = -std::numeric_limits<RealScalar>::infinity(); v(i2,j2) = std::numeric_limits<RealScalar>::quiet_NaN(); - VERIFY(!numext::isfinite(v.squaredNorm())); VERIFY(numext::isnan(v.squaredNorm())); - VERIFY(!numext::isfinite(v.norm())); VERIFY(numext::isnan(v.norm())); - VERIFY(!numext::isfinite(v.stableNorm())); VERIFY(numext::isnan(v.stableNorm())); - VERIFY(!numext::isfinite(v.blueNorm())); VERIFY(numext::isnan(v.blueNorm())); - VERIFY(!numext::isfinite(v.hypotNorm())); VERIFY(numext::isnan(v.hypotNorm())); + VERIFY(!(numext::isfinite)(v.squaredNorm())); VERIFY((numext::isnan)(v.squaredNorm())); + VERIFY(!(numext::isfinite)(v.norm())); VERIFY((numext::isnan)(v.norm())); + VERIFY(!(numext::isfinite)(v.stableNorm())); VERIFY((numext::isnan)(v.stableNorm())); + VERIFY(!(numext::isfinite)(v.blueNorm())); VERIFY((numext::isnan)(v.blueNorm())); + VERIFY(!(numext::isfinite)(v.hypotNorm())); VERIFY((numext::isnan)(v.hypotNorm())); } } |