diff options
author | Gael Guennebaud <g.gael@free.fr> | 2015-10-09 14:57:19 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2015-10-09 14:57:19 +0200 |
commit | ae38910693d9ff96eb5be3ac6cd9a075fc8f7da9 (patch) | |
tree | b713a5280abf54df3e057342d6bf9b4457de5c06 /test | |
parent | 515ecddb97da7f239747c8d547bfb9ed5d63ae0f (diff) |
The evalautor of Solve was missing the EvalBeforeNestingBit flag.
Diffstat (limited to 'test')
-rw-r--r-- | test/nesting_ops.cpp | 26 |
1 files changed, 21 insertions, 5 deletions
diff --git a/test/nesting_ops.cpp b/test/nesting_ops.cpp index 15b231d2e..1e96075f8 100644 --- a/test/nesting_ops.cpp +++ b/test/nesting_ops.cpp @@ -47,6 +47,7 @@ template <typename MatrixType> void run_nesting_ops_1(const MatrixType& _m) template <typename MatrixType> void run_nesting_ops_2(const MatrixType& _m) { + typedef typename MatrixType::Scalar Scalar; Index rows = _m.rows(); Index cols = _m.cols(); MatrixType m1 = MatrixType::Random(rows,cols); @@ -55,7 +56,7 @@ template <typename MatrixType> void run_nesting_ops_2(const MatrixType& _m) { VERIFY_EVALUATION_COUNT( use_n_times<10>(m1), 0 ); - if(!NumTraits<typename MatrixType::Scalar>::IsComplex) + if(!NumTraits<Scalar>::IsComplex) { VERIFY_EVALUATION_COUNT( use_n_times<3>(2*m1), 0 ); VERIFY_EVALUATION_COUNT( use_n_times<4>(2*m1), 1 ); @@ -69,11 +70,21 @@ template <typename MatrixType> void run_nesting_ops_2(const MatrixType& _m) VERIFY_EVALUATION_COUNT( use_n_times<3>(m1+m1), 1 ); VERIFY_EVALUATION_COUNT( use_n_times<1>(m1*m1.transpose()), 1 ); VERIFY_EVALUATION_COUNT( use_n_times<2>(m1*m1.transpose()), 1 ); + + VERIFY_EVALUATION_COUNT( use_n_times<1>(m1 + m1*m1), 2 ); // FIXME should already be 1 thanks the already existing rule + VERIFY_EVALUATION_COUNT( use_n_times<10>(m1 + m1*m1), 2 ); + + VERIFY_EVALUATION_COUNT( use_n_times<1>(m1.template triangularView<Lower>().solve(m1.col(0))), 1 ); + VERIFY_EVALUATION_COUNT( use_n_times<10>(m1.template triangularView<Lower>().solve(m1.col(0))), 1 ); + + VERIFY_EVALUATION_COUNT( use_n_times<1>(Scalar(2)*m1.template triangularView<Lower>().solve(m1.col(0))), 2 ); // FIXME could be one by applying the scaling in-place on the solve result + VERIFY_EVALUATION_COUNT( use_n_times<1>(m1.col(0)+m1.template triangularView<Lower>().solve(m1.col(0))), 2 ); // FIXME could be one by adding m1.col() inplace + VERIFY_EVALUATION_COUNT( use_n_times<10>(m1.col(0)+m1.template triangularView<Lower>().solve(m1.col(0))), 2 ); } { VERIFY( verify_eval_type<10>(m1, m1) ); - if(!NumTraits<typename MatrixType::Scalar>::IsComplex) + if(!NumTraits<Scalar>::IsComplex) { VERIFY( verify_eval_type<3>(2*m1, 2*m1) ); VERIFY( verify_eval_type<4>(2*m1, m1) ); @@ -88,6 +99,10 @@ template <typename MatrixType> void run_nesting_ops_2(const MatrixType& _m) VERIFY( verify_eval_type<1>(m1*m1.transpose(), m1) ); VERIFY( verify_eval_type<1>(m1*(m1+m1).transpose(), m1) ); VERIFY( verify_eval_type<2>(m1*m1.transpose(), m1) ); + VERIFY( verify_eval_type<1>(m1+m1*m1, m1) ); + + VERIFY( verify_eval_type<1>(m1.template triangularView<Lower>().solve(m1), m1) ); + VERIFY( verify_eval_type<1>(m1+m1.template triangularView<Lower>().solve(m1), m1) ); } } @@ -99,9 +114,10 @@ void test_nesting_ops() CALL_SUBTEST_3(run_nesting_ops_1(Matrix4f::Random())); CALL_SUBTEST_4(run_nesting_ops_1(Matrix2d::Random())); - CALL_SUBTEST_1( run_nesting_ops_2(MatrixXf(internal::random<int>(1,EIGEN_TEST_MAX_SIZE), internal::random<int>(1,EIGEN_TEST_MAX_SIZE))) ); - CALL_SUBTEST_2( run_nesting_ops_2(MatrixXcd(internal::random<int>(1,EIGEN_TEST_MAX_SIZE), internal::random<int>(1,EIGEN_TEST_MAX_SIZE))) ); + Index s = internal::random<int>(1,EIGEN_TEST_MAX_SIZE); + CALL_SUBTEST_1( run_nesting_ops_2(MatrixXf(s,s)) ); + CALL_SUBTEST_2( run_nesting_ops_2(MatrixXcd(s,s)) ); CALL_SUBTEST_3( run_nesting_ops_2(Matrix4f()) ); CALL_SUBTEST_4( run_nesting_ops_2(Matrix2d()) ); - + TEST_SET_BUT_UNUSED_VARIABLE(s) } |