diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2010-11-26 15:31:47 +0100 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2010-11-26 15:31:47 +0100 |
| commit | 0d6321225786bc3d95f1dbe5236f07c5e5e96179 (patch) | |
| tree | 001ac35c045bc44f7f59d13f1e0afbaf76f7e3ed /test/eigensolver_selfadjoint.cpp | |
| parent | 421b2b5ff7796c28868dc3b370ca8ceeeda20772 (diff) | |
add a TridiagonalizationMatrixTReturnType class to make Tridiagonalization::matrixT() more efficient and future proof.
Diffstat (limited to 'test/eigensolver_selfadjoint.cpp')
| -rw-r--r-- | test/eigensolver_selfadjoint.cpp | 5 |
1 files changed, 5 insertions, 0 deletions
diff --git a/test/eigensolver_selfadjoint.cpp b/test/eigensolver_selfadjoint.cpp index 806fe5831..b85bcc289 100644 --- a/test/eigensolver_selfadjoint.cpp +++ b/test/eigensolver_selfadjoint.cpp @@ -155,6 +155,11 @@ template<typename MatrixType> void selfadjointeigensolver(const MatrixType& m) VERIFY_RAISES_ASSERT(eiSymmUninitialized.operatorSqrt()); VERIFY_RAISES_ASSERT(eiSymmUninitialized.operatorInverseSqrt()); + // test Tridiagonalization's methods + Tridiagonalization<MatrixType> tridiag(symmA); + // FIXME tridiag.matrixQ().adjoint() does not work + VERIFY_IS_APPROX(MatrixType(symmA.template selfadjointView<Lower>()), tridiag.matrixQ() * tridiag.matrixT().eval() * MatrixType(tridiag.matrixQ()).adjoint()); + if (rows > 1) { // Test matrix with NaN |