fix += return by value like operations

2 files changed, 16 insertions, 1 deletions

diff --git a/test/cholesky.cpp b/test/cholesky.cpp
index 0edf9a793..46140bb11 100644
--- a/test/cholesky.cpp
+++ b/test/cholesky.cpp
@@ -168,6 +168,21 @@ template<typename MatrixType> void cholesky(const MatrixType& m)
     }
   }
 
+  // test some special use cases of SelfCwiseBinaryOp:
+  MatrixType m1 = MatrixType::Random(rows,cols), m2(rows,cols);
+  m2 = m1;
+  m2 += symmLo.template selfadjointView<Lower>().llt().solve(matB);
+  VERIFY_IS_APPROX(m2, m1 + symmLo.template selfadjointView<Lower>().llt().solve(matB));
+  m2 = m1;
+  m2 -= symmLo.template selfadjointView<Lower>().llt().solve(matB);
+  VERIFY_IS_APPROX(m2, m1 - symmLo.template selfadjointView<Lower>().llt().solve(matB));
+  m2 = m1;
+  m2.noalias() += symmLo.template selfadjointView<Lower>().llt().solve(matB);
+  VERIFY_IS_APPROX(m2, m1 + symmLo.template selfadjointView<Lower>().llt().solve(matB));
+  m2 = m1;
+  m2.noalias() -= symmLo.template selfadjointView<Lower>().llt().solve(matB);
+  VERIFY_IS_APPROX(m2, m1 - symmLo.template selfadjointView<Lower>().llt().solve(matB));
+  
 }
 
 template<typename MatrixType> void cholesky_cplx(const MatrixType& m)
diff --git a/test/linearstructure.cpp b/test/linearstructure.cpp
index a0b8982dd..b5c58bdaa 100644
--- a/test/linearstructure.cpp
+++ b/test/linearstructure.cpp
@@ -27,7 +27,7 @@
 template<typename MatrixType> void linearStructure(const MatrixType& m)
 {
   /* this test covers the following files:
-     Sum.h Difference.h Opposite.h ScalarMultiple.h
+     CwiseUnaryOp.h, CwiseBinaryOp.h, SelfCwiseBinaryOp.h 
   */
   typedef typename MatrixType::Index Index;
   typedef typename MatrixType::Scalar Scalar;