backport fix on 3x3 triadiagonalization, this fix #149

This commit is contained in:
Gael Guennebaud 2010-07-22 21:26:09 +02:00
parent 2d0963fb00
commit 70f355b51a

View File

@ -391,7 +391,7 @@ void Tridiagonalization<MatrixType>::_decomposeInPlace3x3(MatrixType& mat, Diago
{
diag[0] = ei_real(mat(0,0));
RealScalar v1norm2 = ei_abs2(mat(0,2));
if (ei_isMuchSmallerThan(v1norm2, RealScalar(1)))
if (v1norm2==RealScalar(0))
{
diag[1] = ei_real(mat(1,1));
diag[2] = ei_real(mat(2,2));