diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2010-07-17 11:57:54 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2010-07-17 11:57:54 +0200 |
| commit | 2d78023815e623b026852d8369433711bf18b07e (patch) | |
| tree | 85906d697088838c30502aca94d2e3b21ed8fc47 /Eigen/src/Core/products/SelfadjointMatrixMatrix.h | |
| parent | cbd6fe323c77ae93c3fbd7ee300a2deeebf39554 (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.h | 2 |
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) |