Fix generalized symm eigensolver (I don't know why the eigenvectors were normalized)

This commit is contained in:
Gael Guennebaud 2010-06-10 00:04:33 +02:00
parent bcf738811e
commit 8692ccc5fb

View File

@ -523,8 +523,6 @@ compute(const MatrixType& matA, const MatrixType& matB, bool computeEigenvectors
{ {
// transform back the eigen vectors: evecs = inv(U) * evecs // transform back the eigen vectors: evecs = inv(U) * evecs
cholB.matrixU().solveInPlace(m_eivec); cholB.matrixU().solveInPlace(m_eivec);
for (Index i=0; i<m_eivec.cols(); ++i)
m_eivec.col(i) = m_eivec.col(i).normalized();
} }
return *this; return *this;
} }