diff options
| author |  Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2014-09-24 14:37:13 +0200 |
|---|---|---|
| committer | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2014-09-24 14:37:13 +0200 |
| commit | 27d6b4daf90488e105d70f45e5d5684b8f1a906d (patch) | |
| tree | 53d709495df6ad1b25f7159da29a1f237d07a126 /test/eigensolver_selfadjoint.cpp | |
| parent | 446001ef51be920649cb4f3c07848967b6788532 (diff) | |
Tridiagonalization::diagonal() and ::subDiagonal() did not work. Added unit-test
Diffstat (limited to 'test/eigensolver_selfadjoint.cpp')
| -rw-r--r-- | test/eigensolver_selfadjoint.cpp | 11 |
1 files changed, 10 insertions, 1 deletions
diff --git a/test/eigensolver_selfadjoint.cpp b/test/eigensolver_selfadjoint.cpp index 3851f9df2..41b598795 100644 --- a/test/eigensolver_selfadjoint.cpp +++ b/test/eigensolver_selfadjoint.cpp @@ -111,8 +111,17 @@ template<typename MatrixType> void selfadjointeigensolver(const MatrixType& m) // test Tridiagonalization's methods Tridiagonalization<MatrixType> tridiag(symmC); - // FIXME tridiag.matrixQ().adjoint() does not work + VERIFY_IS_APPROX(tridiag.diagonal(), tridiag.matrixT().template diagonal()); + VERIFY_IS_APPROX(tridiag.subDiagonal(), tridiag.matrixT().template diagonal<-1>()); + MatrixType T = tridiag.matrixT(); + if(rows>1 && cols>1) { + // FIXME check that upper and lower part are 0: + //VERIFY(T.topRightCorner(rows-2, cols-2).template triangularView<Upper>().isZero()); + } + VERIFY_IS_APPROX(tridiag.diagonal(), T.diagonal().real()); + VERIFY_IS_APPROX(tridiag.subDiagonal(), T.template diagonal<1>().real()); VERIFY_IS_APPROX(MatrixType(symmC.template selfadjointView<Lower>()), tridiag.matrixQ() * tridiag.matrixT().eval() * MatrixType(tridiag.matrixQ()).adjoint()); + VERIFY_IS_APPROX(MatrixType(symmC.template selfadjointView<Lower>()), tridiag.matrixQ() * tridiag.matrixT() * tridiag.matrixQ().adjoint()); // Test computation of eigenvalues from tridiagonal matrix if(rows > 1) |