From 346c00f4c88ecfd35d94dbc19eed300021a0580d Mon Sep 17 00:00:00 2001 From: Benoit Jacob Date: Mon, 3 Dec 2007 10:23:08 +0000 Subject: Rework the unit-tests to use lower precision, so as to eliminate false positives. Also some cleanup in the fuzzy compares. --- src/Core/Fuzzy.h | 6 ++-- src/Core/MatrixBase.h | 17 +--------- test/adjoint.cpp | 41 +++++++++++------------ test/basicstuff.cpp | 90 ++++++++++++++++++++++++++------------------------- test/main.cpp | 13 +++----- test/main.h | 68 ++++++++++++++++++++++++++++++++++++++ 6 files changed, 143 insertions(+), 92 deletions(-) diff --git a/src/Core/Fuzzy.h b/src/Core/Fuzzy.h index 8aa6dd8b2..c1e4ddd74 100644 --- a/src/Core/Fuzzy.h +++ b/src/Core/Fuzzy.h @@ -50,18 +50,18 @@ bool MatrixBase::isApprox( template bool MatrixBase::isMuchSmallerThan( - const Scalar& other, + const typename NumTraits::Real& other, const typename NumTraits::Real& prec ) const { if(IsVector) { - return(norm2() <= abs2(other) * prec * prec); + return(norm2() <= abs2(other * prec)); } else { for(int i = 0; i < cols(); i++) - if(col(i).norm2() > abs2(other) * prec * prec) + if(col(i).norm2() > abs2(other * prec)) return false; return true; } diff --git a/src/Core/MatrixBase.h b/src/Core/MatrixBase.h index f9a668b06..c4588dd88 100644 --- a/src/Core/MatrixBase.h +++ b/src/Core/MatrixBase.h @@ -34,21 +34,6 @@ template class MatrixBase template void _copy_helper(const MatrixBase& other); - template - bool _isApprox_helper( - const OtherDerived& other, - const typename NumTraits::Real& prec = precision() - ) const; - bool _isMuchSmallerThan_helper( - const Scalar& other, - const typename NumTraits::Real& prec = precision() - ) const; - template - bool _isMuchSmallerThan_helper( - const MatrixBase& other, - const typename NumTraits::Real& prec = precision() - ) const; - public: static const int SizeAtCompileTime = RowsAtCompileTime == Dynamic || ColsAtCompileTime == Dynamic @@ -124,7 +109,7 @@ template class MatrixBase const typename NumTraits::Real& prec = precision() ) const; bool isMuchSmallerThan( - const Scalar& other, + const typename NumTraits::Real& other, const typename NumTraits::Real& prec = precision() ) const; template diff --git a/test/adjoint.cpp b/test/adjoint.cpp index d34137c10..a86225408 100644 --- a/test/adjoint.cpp +++ b/test/adjoint.cpp @@ -55,37 +55,38 @@ template void adjoint(const MatrixType& m) s2 = random(); // check involutivity of adjoint, transpose, conjugate - QVERIFY(m1.transpose().transpose().isApprox(m1)); - QVERIFY(m1.conjugate().conjugate().isApprox(m1)); - QVERIFY(m1.adjoint().adjoint().isApprox(m1)); + VERIFY_IS_APPROX(m1.transpose().transpose(), m1); + VERIFY_IS_APPROX(m1.conjugate().conjugate(), m1); + VERIFY_IS_APPROX(m1.adjoint().adjoint(), m1); // check basic compatibility of adjoint, transpose, conjugate - QVERIFY(m1.transpose().conjugate().adjoint().isApprox(m1)); - QVERIFY(m1.adjoint().conjugate().transpose().isApprox(m1)); - if(!NumTraits::IsComplex) QVERIFY(m1.adjoint().transpose().isApprox(m1)); + VERIFY_IS_APPROX(m1.transpose().conjugate().adjoint(), m1); + VERIFY_IS_APPROX(m1.adjoint().conjugate().transpose(), m1); + if(!NumTraits::IsComplex) + VERIFY_IS_APPROX(m1.adjoint().transpose(), m1); // check multiplicative behavior - QVERIFY((m1.transpose() * m2).transpose().isApprox(m2.transpose() * m1)); - QVERIFY((m1.adjoint() * m2).adjoint().isApprox(m2.adjoint() * m1)); - QVERIFY((m1.transpose() * m2).conjugate().isApprox(m1.adjoint() * m2.conjugate())); - QVERIFY((s1 * m1).transpose().isApprox(s1 * m1.transpose())); - QVERIFY((s1 * m1).conjugate().isApprox(conj(s1) * m1.conjugate())); - QVERIFY((s1 * m1).adjoint().isApprox(conj(s1) * m1.adjoint())); + VERIFY_IS_APPROX((m1.transpose() * m2).transpose(), m2.transpose() * m1); + VERIFY_IS_APPROX((m1.adjoint() * m2).adjoint(), m2.adjoint() * m1); + VERIFY_IS_APPROX((m1.transpose() * m2).conjugate(), m1.adjoint() * m2.conjugate()); + VERIFY_IS_APPROX((s1 * m1).transpose(), s1 * m1.transpose()); + VERIFY_IS_APPROX((s1 * m1).conjugate(), conj(s1) * m1.conjugate()); + VERIFY_IS_APPROX((s1 * m1).adjoint(), conj(s1) * m1.adjoint()); // check basic properties of dot, norm, norm2 typedef typename NumTraits::Real RealScalar; - QVERIFY(isApprox((s1 * v1 + s2 * v2).dot(v3), s1 * v1.dot(v3) + s2 * v2.dot(v3))); - QVERIFY(isApprox(v3.dot(s1 * v1 + s2 * v2), conj(s1) * v3.dot(v1) + conj(s2) * v3.dot(v2))); - QVERIFY(isApprox(conj(v1.dot(v2)), v2.dot(v1))); - QVERIFY(isApprox(abs(v1.dot(v1)), v1.norm2())); + VERIFY_IS_APPROX((s1 * v1 + s2 * v2).dot(v3), s1 * v1.dot(v3) + s2 * v2.dot(v3)); + VERIFY_IS_APPROX(v3.dot(s1 * v1 + s2 * v2), conj(s1)*v3.dot(v1)+conj(s2)*v3.dot(v2)); + VERIFY_IS_APPROX(conj(v1.dot(v2)), v2.dot(v1)); + VERIFY_IS_APPROX(abs(v1.dot(v1)), v1.norm2()); if(NumTraits::HasFloatingPoint) - QVERIFY(isApprox(v1.norm2(), v1.norm() * v1.norm())); - QVERIFY(isMuchSmallerThan(abs(vzero.dot(v1)), static_cast(1))); + VERIFY_IS_APPROX(v1.norm2(), v1.norm() * v1.norm()); + VERIFY_IS_MUCH_SMALLER_THAN(abs(vzero.dot(v1)), static_cast(1)); if(NumTraits::HasFloatingPoint) - QVERIFY(isMuchSmallerThan(vzero.norm(), static_cast(1))); + VERIFY_IS_MUCH_SMALLER_THAN(vzero.norm(), static_cast(1)); // check compatibility of dot and adjoint - QVERIFY(isApprox(v1.dot(square * v2), (square.adjoint() * v1).dot(v2))); + VERIFY_IS_APPROX(v1.dot(square * v2), (square.adjoint() * v1).dot(v2)); } void EigenTest::testAdjoint() diff --git a/test/basicstuff.cpp b/test/basicstuff.cpp index f28477a68..03f555749 100644 --- a/test/basicstuff.cpp +++ b/test/basicstuff.cpp @@ -62,64 +62,66 @@ template void basicStuff(const MatrixType& m) s2 = random(); // test Fuzzy.h and Zero.h. - QVERIFY(v1.isApprox(v1)); - QVERIFY(!v1.isApprox(2*v1)); - QVERIFY(vzero.isMuchSmallerThan(v1)); + VERIFY_IS_APPROX( v1, v1); + VERIFY_IS_NOT_APPROX( v1, 2*v1); + VERIFY_IS_MUCH_SMALLER_THAN( vzero, v1); if(NumTraits::HasFloatingPoint) - QVERIFY(vzero.isMuchSmallerThan(v1.norm())); - QVERIFY(!v1.isMuchSmallerThan(v1)); - QVERIFY(vzero.isApprox(v1-v1)); - QVERIFY(m1.isApprox(m1)); - QVERIFY(!m1.isApprox(2*m1)); - QVERIFY(mzero.isMuchSmallerThan(m1)); - QVERIFY(!m1.isMuchSmallerThan(m1)); - QVERIFY(mzero.isApprox(m1-m1)); + VERIFY_IS_MUCH_SMALLER_THAN( vzero, v1.norm()); + VERIFY_IS_NOT_MUCH_SMALLER_THAN(v1, v1); + VERIFY_IS_APPROX( vzero, v1-v1); + VERIFY_IS_APPROX( m1, m1); + VERIFY_IS_NOT_APPROX( m1, 2*m1); + VERIFY_IS_MUCH_SMALLER_THAN( mzero, m1); + VERIFY_IS_NOT_MUCH_SMALLER_THAN(m1, m1); + VERIFY_IS_APPROX( mzero, m1-m1); // test the linear structure, i.e. the following files: // Sum.h Difference.h Opposite.h ScalarMultiple.h - QVERIFY((m1+m1).isApprox(2 * m1)); - QVERIFY((m1+m2-m1).isApprox(m2)); - QVERIFY((-m2+m1+m2).isApprox(m1)); - QVERIFY((m1 * s1).isApprox(s1 * m1)); - QVERIFY(((m1 + m2) * s1).isApprox(s1 * m1 + s1 * m2)); - QVERIFY(((s1 + s2) * m1).isApprox(m1 * s1 + m1 * s2)); - QVERIFY(((m1 - m2) * s1).isApprox(s1 * m1 - s1 * m2)); - QVERIFY(((s1 - s2) * m1).isApprox(m1 * s1 - m1 * s2)); - QVERIFY(((-m1 + m2) * s1).isApprox(-s1 * m1 + s1 * m2)); - QVERIFY(((-s1 + s2) * m1).isApprox(-m1 * s1 + m1 * s2)); - m3 = m2; - QVERIFY((m3 += m1).isApprox(m1 + m2)); - m3 = m2; - QVERIFY((m3 -= m1).isApprox(-m1 + m2)); - m3 = m2; - QVERIFY((m3 *= s1).isApprox(s1 * m2)); - m3 = m2; - if(NumTraits::HasFloatingPoint - && s1 != static_cast(0)) - QVERIFY((m3 /= s1).isApprox(m2 / s1)); + VERIFY_IS_APPROX(-(-m1), m1); + VERIFY_IS_APPROX(m1+m1, 2*m1); + VERIFY_IS_APPROX(m1+m2-m1, m2); + VERIFY_IS_APPROX(-m2+m1+m2, m1); + VERIFY_IS_APPROX(m1*s1, s1*m1); + VERIFY_IS_APPROX((m1+m2)*s1, s1*m1+s1*m2); + VERIFY_IS_APPROX((s1+s2)*m1, m1*s1+m1*s2); + VERIFY_IS_APPROX((m1-m2)*s1, s1*m1-s1*m2); + VERIFY_IS_APPROX((s1-s2)*m1, m1*s1-m1*s2); + VERIFY_IS_APPROX((-m1+m2)*s1, -s1*m1+s1*m2); + VERIFY_IS_APPROX((-s1+s2)*m1, -m1*s1+m1*s2); + m3 = m2; m3 += m1; + VERIFY_IS_APPROX(m3, m1+m2); + m3 = m2; m3 -= m1; + VERIFY_IS_APPROX(m3, m2-m1); + m3 = m2; m3 *= s1; + VERIFY_IS_APPROX(m3, s1*m2); + if(NumTraits::HasFloatingPoint) + { + m3 = m2; m3 /= s1; + VERIFY_IS_APPROX(m3, m2/s1); + } // begin testing Product.h: only associativity for now // (we use Transpose.h but this doesn't count as a test for it) - QVERIFY(((m1 * m1.transpose()) * m2).isApprox(m1 * (m1.transpose() * m2))); + VERIFY_IS_APPROX((m1*m1.transpose())*m2, m1*(m1.transpose()*m2)); m3 = m1; m3 *= (m1.transpose() * m2); - QVERIFY(m3.isApprox(m1 * (m1.transpose() * m2))); - QVERIFY(m3.isApprox(m1.lazyProduct(m1.transpose() * m2))); + VERIFY_IS_APPROX(m3, m1*(m1.transpose()*m2)); + VERIFY_IS_APPROX(m3, m1.lazyProduct(m1.transpose()*m2)); // continue testing Product.h: distributivity - QVERIFY((square * (m1 + m2)).isApprox(square * m1 + square * m2)); - QVERIFY((square * (m1 - m2)).isApprox(square * m1 - square * m2)); + VERIFY_IS_APPROX(square*(m1 + m2), square*m1+square*m2); + VERIFY_IS_APPROX(square*(m1 - m2), square*m1-square*m2); // continue testing Product.h: compatibility with ScalarMultiple.h - QVERIFY((s1 * (square * m1)).isApprox((s1 * square) * m1)); - QVERIFY((s1 * (square * m1)).isApprox(square * (m1 * s1))); + VERIFY_IS_APPROX(s1*(square*m1), (s1*square)*m1); + VERIFY_IS_APPROX(s1*(square*m1), square*(m1*s1)); // continue testing Product.h: lazyProduct - QVERIFY(square.lazyProduct(m1).isApprox(square * m1)); + VERIFY_IS_APPROX(square.lazyProduct(m1), square*m1); // test Product.h together with Identity.h. This does test Identity.h. - QVERIFY(m1.isApprox(identity * m1)); - QVERIFY(v1.isApprox(identity * v1)); + VERIFY_IS_APPROX(m1, identity*m1); + VERIFY_IS_APPROX(v1, identity*v1); // test FromArray.h Scalar* array1 = new Scalar[rows]; @@ -127,9 +129,9 @@ template void basicStuff(const MatrixType& m) Matrix::fromArray(array1, rows) = Matrix::random(rows); Matrix::fromArray(array2, rows) = Matrix::fromArray(array1, rows); - bool b = Matrix::fromArray(array1, rows) - .isApprox(Matrix::fromArray(array2, rows)); - QVERIFY(b); + Matrix ma1 = Matrix::fromArray(array1, rows); + Matrix ma2 = Matrix::fromArray(array2, rows); + VERIFY_IS_APPROX(ma1, ma2); delete[] array1; delete[] array2; } diff --git a/test/main.cpp b/test/main.cpp index 3a63653a4..91f83e396 100644 --- a/test/main.cpp +++ b/test/main.cpp @@ -31,7 +31,7 @@ int main(int argc, char *argv[]) bool has_set_repeat = false; bool has_set_seed = false; - bool want_help = false; + bool need_help = false; unsigned int seed; int repeat; @@ -70,23 +70,18 @@ int main(int argc, char *argv[]) return 1; } } - else if(arg == "h" || arg == "-h" || arg.contains("help", Qt::CaseInsensitive)) - { - want_help = true; - } else { - qDebug() << "Invalid command-line argument" << arg; - return 1; + need_help = true; } } - if(want_help) + if(need_help) { qDebug() << "This test application takes the following optional arguments:"; qDebug() << " rN Repeat each test N times (default:" << DEFAULT_REPEAT << ")"; qDebug() << " sN Use N as seed for random numbers (default: based on current time)"; - return 0; + return 1; } if(!has_set_seed) seed = (unsigned int) time(NULL); diff --git a/test/main.h b/test/main.h index 4bb7b1b15..00048d237 100644 --- a/test/main.h +++ b/test/main.h @@ -37,6 +37,74 @@ namespace Eigen { +template inline typename NumTraits::Real test_precision(); +template<> inline int test_precision() { return 0; } +template<> inline float test_precision() { return 1e-2; } +template<> inline double test_precision() { return 1e-5; } +template<> inline float test_precision >() { return test_precision(); } +template<> inline double test_precision >() { return test_precision(); } + +inline bool test_isApprox(const int& a, const int& b) +{ return isApprox(a, b, test_precision()); } +inline bool test_isMuchSmallerThan(const int& a, const int& b) +{ return isMuchSmallerThan(a, b, test_precision()); } +inline bool test_isApproxOrLessThan(const int& a, const int& b) +{ return isApproxOrLessThan(a, b, test_precision()); } + +inline bool test_isApprox(const float& a, const float& b) +{ return isApprox(a, b, test_precision()); } +inline bool test_isMuchSmallerThan(const float& a, const float& b) +{ return isMuchSmallerThan(a, b, test_precision()); } +inline bool test_isApproxOrLessThan(const float& a, const float& b) +{ return isApproxOrLessThan(a, b, test_precision()); } + +inline bool test_isApprox(const double& a, const double& b) +{ return isApprox(a, b, test_precision()); } +inline bool test_isMuchSmallerThan(const double& a, const double& b) +{ return isMuchSmallerThan(a, b, test_precision()); } +inline bool test_isApproxOrLessThan(const double& a, const double& b) +{ return isApproxOrLessThan(a, b, test_precision()); } + +inline bool test_isApprox(const std::complex& a, const std::complex& b) +{ return isApprox(a, b, test_precision >()); } +inline bool test_isMuchSmallerThan(const std::complex& a, const std::complex& b) +{ return isMuchSmallerThan(a, b, test_precision >()); } + +inline bool test_isApprox(const std::complex& a, const std::complex& b) +{ return isApprox(a, b, test_precision >()); } +inline bool test_isMuchSmallerThan(const std::complex& a, const std::complex& b) +{ return isMuchSmallerThan(a, b, test_precision >()); } + +template +inline bool test_isApprox(const MatrixBase& m1, + const MatrixBase& m2) +{ + return m1.isApprox(m2, test_precision()); +} + +template +inline bool test_isMuchSmallerThan(const MatrixBase& m1, + const MatrixBase& m2) +{ + return m1.isMuchSmallerThan(m2, test_precision()); +} + +template +inline bool test_isMuchSmallerThan(const MatrixBase& m, + const typename NumTraits::Real& s) +{ + return m.isMuchSmallerThan(s, test_precision()); +} + + +#define VERIFY(a) QVERIFY(a) +#define VERIFY_IS_APPROX(a, b) QVERIFY(test_isApprox(a, b)) +#define VERIFY_IS_NOT_APPROX(a, b) QVERIFY(!test_isApprox(a, b)) +#define VERIFY_IS_MUCH_SMALLER_THAN(a, b) QVERIFY(test_isMuchSmallerThan(a, b)) +#define VERIFY_IS_NOT_MUCH_SMALLER_THAN(a, b) QVERIFY(!test_isMuchSmallerThan(a, b)) +#define VERIFY_IS_APPROX_OR_LESS_THAN(a, b) QVERIFY(test_isApproxOrLessThan(a, b)) +#define VERIFY_IS_NOT_APPROX_OR_LESS_THAN(a, b) QVERIFY(!test_isApproxOrLessThan(a, b)) + class EigenTest : public QObject { Q_OBJECT -- cgit v1.2.3