GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-13 04:09:10 +08:00
Code Issues Packages Projects Releases Wiki Activity

bugfix in the eigenvalue solvers (forgot to resize the eigen vectors)

Browse Source
This commit is contained in:
Gael Guennebaud 2009-09-30 16:48:02 +02:00
parent 4b04a9bcfa
commit d7a2a37a4c
3 changed files with 3 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
Eigen/src/Eigenvalues/ComplexEigenSolver.h
View File

@ -92,6 +92,7 @@ void ComplexEigenSolver<MatrixType>::compute(const MatrixType& matrix)
assert(matrix.cols() == matrix.rows()); assert(matrix.cols() == matrix.rows());
int n = matrix.cols(); int n = matrix.cols();
m_eivalues.resize(n,1); m_eivalues.resize(n,1);
m_eivec.resize(n,n);
RealScalar eps = epsilon<RealScalar>(); RealScalar eps = epsilon<RealScalar>();

1
Eigen/src/Eigenvalues/EigenSolver.h
View File

@ -194,6 +194,7 @@ EigenSolver<MatrixType>& EigenSolver<MatrixType>::compute(const MatrixType& matr
assert(matrix.cols() == matrix.rows()); assert(matrix.cols() == matrix.rows());
int n = matrix.cols(); int n = matrix.cols();
m_eivalues.resize(n,1); m_eivalues.resize(n,1);
m_eivec.resize(n,n);
MatrixType matH = matrix; MatrixType matH = matrix;
RealVectorType ort(n); RealVectorType ort(n);

1
Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h
View File

@ -168,6 +168,7 @@ SelfAdjointEigenSolver<MatrixType>& SelfAdjointEigenSolver<MatrixType>::compute(
assert(matrix.cols() == matrix.rows()); assert(matrix.cols() == matrix.rows());
int n = matrix.cols(); int n = matrix.cols();
m_eivalues.resize(n,1); m_eivalues.resize(n,1);
m_eivec.resize(n,n);
if(n==1) if(n==1)
{ {