From 70f355b51a8e98d999711e6458034603dea1cf5a Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Thu, 22 Jul 2010 21:26:09 +0200 Subject: [PATCH] backport fix on 3x3 triadiagonalization, this fix #149 --- Eigen/src/QR/Tridiagonalization.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Eigen/src/QR/Tridiagonalization.h b/Eigen/src/QR/Tridiagonalization.h index 69b395fe9..398535ba1 100644 --- a/Eigen/src/QR/Tridiagonalization.h +++ b/Eigen/src/QR/Tridiagonalization.h @@ -391,7 +391,7 @@ void Tridiagonalization::_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));