Resolve aliasing issue (relates #49)

lisitsyn · Feb 1, 2019 · ec35140 · ec35140
1 parent 1509ae8
commit ec35140
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/include/tapkee/routines/eigendecomposition.hpp b/include/tapkee/routines/eigendecomposition.hpp
@@ -57,7 +57,7 @@ EigendecompositionResult eigendecomposition_impl_dense(const MatrixType& wm, Ind
 	timed_context context("Eigen library dense eigendecomposition");
 
 	DenseSymmetricMatrix dense_wm = wm;
-	dense_wm += dense_wm.transpose();
+	dense_wm += dense_wm.transpose().eval();
 	dense_wm /= 2.0;
 	DenseSelfAdjointEigenSolver solver(dense_wm);