aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/Core/products/SelfadjointMatrixMatrix.h
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2010-07-17 11:57:54 +0200
committerGravatar Gael Guennebaud <g.gael@free.fr>2010-07-17 11:57:54 +0200
commit2d78023815e623b026852d8369433711bf18b07e (patch)
tree85906d697088838c30502aca94d2e3b21ed8fc47 /Eigen/src/Core/products/SelfadjointMatrixMatrix.h
parentcbd6fe323c77ae93c3fbd7ee300a2deeebf39554 (diff)
fix hemm to not use the imaginary part of the diagonal entries
Diffstat (limited to 'Eigen/src/Core/products/SelfadjointMatrixMatrix.h')
-rw-r--r--Eigen/src/Core/products/SelfadjointMatrixMatrix.h2
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/Core/products/SelfadjointMatrixMatrix.h b/Eigen/src/Core/products/SelfadjointMatrixMatrix.h
index d8fa1bd9c..4a95ad5e1 100644
--- a/Eigen/src/Core/products/SelfadjointMatrixMatrix.h
+++ b/Eigen/src/Core/products/SelfadjointMatrixMatrix.h
@@ -136,7 +136,7 @@ struct ei_symm_pack_rhs
for (Index w=0 ; w<h; ++w)
blockB[count+w] = alpha*rhs(k,j2+w);
- blockB[count+h] = alpha*rhs(k,k);
+ blockB[count+h] = alpha*ei_real(rhs(k,k));
// transpose
for (Index w=h+1 ; w<nr; ++w)