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

update of the eigeinsolver unit test to check complex

Browse Source
This commit is contained in:
Gael Guennebaud 2008-06-03 18:04:36 +00:00
parent a0cff1a295
commit 42ad9c4352
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
test/eigensolver.cpp
View File

@ -28,7 +28,7 @@
template<typename MatrixType> void eigensolver(const MatrixType& m) template<typename MatrixType> void eigensolver(const MatrixType& m)
{ {
/* this test covers the following files: /* this test covers the following files:
EigenSolver.h EigenSolver.h, SelfAdjointEigenSolver.h (and indirectly: Tridiagonalization.h)
*/ */
int rows = m.rows(); int rows = m.rows();
int cols = m.cols(); int cols = m.cols();
@ -39,13 +39,13 @@ template<typename MatrixType> void eigensolver(const MatrixType& m)
MatrixType covMat = a.adjoint() * a; MatrixType covMat = a.adjoint() * a;
SelfAdjointEigenSolver<MatrixType> eiSymm(covMat); SelfAdjointEigenSolver<MatrixType> eiSymm(covMat);
VERIFY_IS_APPROX(covMat * eiSymm.eigenvectors(), eiSymm.eigenvectors() * eiSymm.eigenvalues().asDiagonal()); VERIFY_IS_APPROX(covMat * eiSymm.eigenvectors(), (eiSymm.eigenvectors() * eiSymm.eigenvalues().asDiagonal().eval()));
EigenSolver<MatrixType> eiNotSymmButSymm(covMat); // EigenSolver<MatrixType> eiNotSymmButSymm(covMat);
VERIFY_IS_APPROX((covMat.template cast<Complex>()) * (eiNotSymmButSymm.eigenvectors().template cast<Complex>()), // VERIFY_IS_APPROX((covMat.template cast<Complex>()) * (eiNotSymmButSymm.eigenvectors().template cast<Complex>()),
(eiNotSymmButSymm.eigenvectors().template cast<Complex>()) * (eiNotSymmButSymm.eigenvalues().asDiagonal())); // (eiNotSymmButSymm.eigenvectors().template cast<Complex>()) * (eiNotSymmButSymm.eigenvalues().asDiagonal()));
EigenSolver<MatrixType> eiNotSymm(a); // EigenSolver<MatrixType> eiNotSymm(a);
// VERIFY_IS_APPROX(a.template cast<Complex>() * eiNotSymm.eigenvectors().template cast<Complex>(), // VERIFY_IS_APPROX(a.template cast<Complex>() * eiNotSymm.eigenvectors().template cast<Complex>(),
// eiNotSymm.eigenvectors().template cast<Complex>() * eiNotSymm.eigenvalues().asDiagonal()); // eiNotSymm.eigenvectors().template cast<Complex>() * eiNotSymm.eigenvalues().asDiagonal());
@ -57,5 +57,6 @@ void test_eigensolver()
CALL_SUBTEST( eigensolver(Matrix3f()) ); CALL_SUBTEST( eigensolver(Matrix3f()) );
CALL_SUBTEST( eigensolver(Matrix4d()) ); CALL_SUBTEST( eigensolver(Matrix4d()) );
CALL_SUBTEST( eigensolver(MatrixXd(7,7)) ); CALL_SUBTEST( eigensolver(MatrixXd(7,7)) );
CALL_SUBTEST( eigensolver(MatrixXcd(6,6)) );
} }
} }