diff options
author | 2009-09-22 01:41:09 -0400 | |
---|---|---|
committer | 2009-09-22 01:41:09 -0400 | |
commit | 176c26feb55c6f3faf7f4402f7dc65431d92c8aa (patch) | |
tree | f9a701c7c5fb5d3bd9ef7044c58eb2bc43ada823 /test | |
parent | 4f9e27034392d4f6744509e52f7b7d829e9070ce (diff) |
allow to do xpr = solve(b) etc... just by adding a dummy MatrixBase::resize()
Diffstat (limited to 'test')
-rw-r--r-- | test/lu.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/test/lu.cpp b/test/lu.cpp index f34c941f7..a7217f5b9 100644 --- a/test/lu.cpp +++ b/test/lu.cpp @@ -40,7 +40,6 @@ template<typename MatrixType> void lu_non_invertible() typename ei_lu_kernel_impl<MatrixType>::ReturnMatrixType m1kernel = lu.kernel(); typename ei_lu_image_impl <MatrixType>::ReturnMatrixType m1image = lu.image(); - std::cout << "rows:" << rows << " cols:" << cols << " | " << rank << " ----- " << lu.rank() << std::endl; VERIFY(rank == lu.rank()); VERIFY(cols - lu.rank() == lu.dimensionOfKernel()); VERIFY(!lu.isInjective()); @@ -54,7 +53,8 @@ template<typename MatrixType> void lu_non_invertible() m2 = MatrixType::Random(cols,cols2); m3 = m1*m2; m2 = MatrixType::Random(cols,cols2); - m2 = lu.solve(m3); + // test that the code, which does resize(), may be applied to an xpr + m2.block(0,0,cols,cols2) = lu.solve(m3); VERIFY_IS_APPROX(m3, m1*m2); typedef Matrix<typename MatrixType::Scalar, MatrixType::RowsAtCompileTime, MatrixType::RowsAtCompileTime> SquareMatrixType; |