diff options
Diffstat (limited to 'test/stable_norm.cpp')
-rw-r--r-- | test/stable_norm.cpp | 11 |
1 files changed, 8 insertions, 3 deletions
diff --git a/test/stable_norm.cpp b/test/stable_norm.cpp index ee5f91674..008e35d87 100644 --- a/test/stable_norm.cpp +++ b/test/stable_norm.cpp @@ -161,8 +161,12 @@ template<typename MatrixType> void stable_norm(const MatrixType& m) // mix { - Index i2 = internal::random<Index>(0,rows-1); - Index j2 = internal::random<Index>(0,cols-1); + // Ensure unique indices otherwise inf may be overwritten by NaN. + Index i2, j2; + do { + i2 = internal::random<Index>(0,rows-1); + j2 = internal::random<Index>(0,cols-1); + } while (i2 == i && j2 == j); v = vrand; v(i,j) = -std::numeric_limits<RealScalar>::infinity(); v(i2,j2) = std::numeric_limits<RealScalar>::quiet_NaN(); @@ -170,7 +174,8 @@ template<typename MatrixType> void stable_norm(const MatrixType& m) 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())); + // hypot propagates inf over NaN. + VERIFY(!(numext::isfinite)(v.hypotNorm())); VERIFY((numext::isinf)(v.hypotNorm())); } // stableNormalize[d] |