From bacd5318629f7a003763d85cff95a7cc1971fa90 Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Thu, 9 Dec 2010 19:56:20 +0100
Subject: [PATCH] fix bug #128 : tridiagonalization failed for 1x1 matrices

---
 Eigen/src/Eigenvalues/Tridiagonalization.h | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/Eigen/src/Eigenvalues/Tridiagonalization.h b/Eigen/src/Eigenvalues/Tridiagonalization.h
index 9a46fc861..66bec20d9 100644
--- a/Eigen/src/Eigenvalues/Tridiagonalization.h
+++ b/Eigen/src/Eigenvalues/Tridiagonalization.h
@@ -356,12 +356,13 @@ namespace internal {
 template<typename MatrixType, typename CoeffVectorType>
 void tridiagonalization_inplace(MatrixType& matA, CoeffVectorType& hCoeffs)
 {
-  eigen_assert(matA.rows()==matA.cols());
-  eigen_assert(matA.rows()==hCoeffs.size()+1);
   typedef typename MatrixType::Index Index;
   typedef typename MatrixType::Scalar Scalar;
   typedef typename MatrixType::RealScalar RealScalar;
   Index n = matA.rows();
+  eigen_assert(n==matA.cols());
+  eigen_assert(n==hCoeffs.size()+1 || n==1);
+  
   for (Index i = 0; i<n-1; ++i)
   {
     Index remainingSize = n-i-1;