fix hemm to not use the imaginary part of the diagonal entries

This commit is contained in:
Gael Guennebaud 2010-07-17 11:57:54 +02:00
parent cbd6fe323c
commit 2d78023815

View File

@ -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)