diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2009-09-30 16:48:02 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2009-09-30 16:48:02 +0200 |
| commit | d7a2a37a4c1728dbddbfae9d577e2e49f5edda65 (patch) | |
| tree | 00a6535b57b74664bb3ad8be76be48131e444bb3 | |
| parent | 4b04a9bcfa1e55a00053588b6a50bb000af2d044 (diff) | |
bugfix in the eigenvalue solvers (forgot to resize the eigen vectors)
| -rw-r--r-- | Eigen/src/Eigenvalues/ComplexEigenSolver.h | 1 | ||||
| -rw-r--r-- | Eigen/src/Eigenvalues/EigenSolver.h | 1 | ||||
| -rw-r--r-- | Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h | 1 |
3 files changed, 3 insertions, 0 deletions
diff --git a/Eigen/src/Eigenvalues/ComplexEigenSolver.h b/Eigen/src/Eigenvalues/ComplexEigenSolver.h index 666381949..224918795 100644 --- a/Eigen/src/Eigenvalues/ComplexEigenSolver.h +++ b/Eigen/src/Eigenvalues/ComplexEigenSolver.h @@ -92,6 +92,7 @@ void ComplexEigenSolver<MatrixType>::compute(const MatrixType& matrix) assert(matrix.cols() == matrix.rows());
int n = matrix.cols();
m_eivalues.resize(n,1);
+ m_eivec.resize(n,n);
RealScalar eps = epsilon<RealScalar>();
diff --git a/Eigen/src/Eigenvalues/EigenSolver.h b/Eigen/src/Eigenvalues/EigenSolver.h index 3fc36c080..73d240de0 100644 --- a/Eigen/src/Eigenvalues/EigenSolver.h +++ b/Eigen/src/Eigenvalues/EigenSolver.h @@ -194,6 +194,7 @@ EigenSolver<MatrixType>& EigenSolver<MatrixType>::compute(const MatrixType& matr assert(matrix.cols() == matrix.rows()); int n = matrix.cols(); m_eivalues.resize(n,1); + m_eivec.resize(n,n); MatrixType matH = matrix; RealVectorType ort(n); diff --git a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h index 84856aa66..9e155de8f 100644 --- a/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h +++ b/Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h @@ -168,6 +168,7 @@ SelfAdjointEigenSolver<MatrixType>& SelfAdjointEigenSolver<MatrixType>::compute( assert(matrix.cols() == matrix.rows()); int n = matrix.cols(); m_eivalues.resize(n,1); + m_eivec.resize(n,n); if(n==1) { |