diff options
Diffstat (limited to 'test')
-rw-r--r-- | test/basicstuff.cpp | 33 |
1 files changed, 25 insertions, 8 deletions
diff --git a/test/basicstuff.cpp b/test/basicstuff.cpp index 1d61b3aaa..c21a58979 100644 --- a/test/basicstuff.cpp +++ b/test/basicstuff.cpp @@ -34,26 +34,38 @@ template<typename MatrixType> void basicStuff(const MatrixType& m) MatrixType m1 = MatrixType::random(rows, cols), m2 = MatrixType::random(rows, cols), - m3; - VectorType v1 = VectorType::random(rows, 1), - v2 = VectorType::random(rows, 1); + m3, + mzero = MatrixType::zero(rows, cols), + identity = Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime> + ::identity(rows, rows), + square = Matrix<Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime> + ::random(rows, rows); + VectorType v1 = VectorType::random(rows), + v2 = VectorType::random(rows), + vzero = VectorType::zero(rows); Scalar s1 = NumTraits<Scalar>::random(), s2 = NumTraits<Scalar>::random(); QVERIFY(v1.isApprox(v1)); - QVERIFY((v1-v1).isMuchSmallerThan(v1)); - QVERIFY((v1-v1).isMuchSmallerThan(v1.norm())); + QVERIFY(!v1.isApprox(2*v1)); + + QVERIFY(vzero.isMuchSmallerThan(v1)); + QVERIFY(vzero.isMuchSmallerThan(v1.norm())); + QVERIFY(!v1.isMuchSmallerThan(v1)); QVERIFY(m1.isApprox(m1)); - QVERIFY((m1-m1).isMuchSmallerThan(m1)); + QVERIFY(!m1.isApprox(2*m1)); + QVERIFY(mzero.isMuchSmallerThan(m1)); + QVERIFY(!m1.isMuchSmallerThan(m1)); + + QVERIFY(vzero.isApprox(v1-v1)); + QVERIFY(mzero.isApprox(m1-m1)); QVERIFY((m1+m1).isApprox(2 * 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)); - - m3 = m2; QVERIFY((m3 += m1).isApprox(m1 + m2)); m3 = m2; @@ -70,6 +82,11 @@ template<typename MatrixType> void basicStuff(const MatrixType& m) m3 *= (m1.transpose() * m2); QVERIFY(m3.isApprox(m1 * (m1.transpose() * m2))); QVERIFY(m3.isApprox(m1.lazyProduct(m1.transpose() * m2))); + + QVERIFY(m1.isApprox(identity * m1)); + QVERIFY(v1.isApprox(identity * v1)); + + QVERIFY((square * (m1 + m2)).isApprox(square * m1 + square * m2)); } void EigenTest::testBasicStuff() |