diff options
author | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2012-07-24 15:17:59 +0100 |
---|---|---|
committer | Jitse Niesen <jitse@maths.leeds.ac.uk> | 2012-07-24 15:17:59 +0100 |
commit | ba5eecae53aa038374d1708573cf03a2df3f76f3 (patch) | |
tree | b1bb109190fc9b31cbbda8c227343170e2058307 /test/eigensolver_generic.cpp | |
parent | b7ac053b9c1196bf775ca0ec45765f6262648c43 (diff) |
Allow user to specify max number of iterations (bug #479).
Diffstat (limited to 'test/eigensolver_generic.cpp')
-rw-r--r-- | test/eigensolver_generic.cpp | 12 |
1 files changed, 11 insertions, 1 deletions
diff --git a/test/eigensolver_generic.cpp b/test/eigensolver_generic.cpp index 0b55ccd93..1b83afbc4 100644 --- a/test/eigensolver_generic.cpp +++ b/test/eigensolver_generic.cpp @@ -2,7 +2,7 @@ // for linear algebra. // // Copyright (C) 2008 Gael Guennebaud <gael.guennebaud@inria.fr> -// Copyright (C) 2010 Jitse Niesen <jitse@maths.leeds.ac.uk> +// Copyright (C) 2010,2012 Jitse Niesen <jitse@maths.leeds.ac.uk> // // This Source Code Form is subject to the terms of the Mozilla // Public License v. 2.0. If a copy of the MPL was not distributed @@ -45,6 +45,16 @@ template<typename MatrixType> void eigensolver(const MatrixType& m) VERIFY_IS_APPROX(ei1.eigenvectors().colwise().norm(), RealVectorType::Ones(rows).transpose()); VERIFY_IS_APPROX(a.eigenvalues(), ei1.eigenvalues()); + EigenSolver<MatrixType> ei2; + ei2.compute(a, true, RealSchur<MatrixType>::m_maxIterations * rows); + VERIFY_IS_EQUAL(ei2.info(), Success); + VERIFY_IS_EQUAL(ei2.eigenvectors(), ei1.eigenvectors()); + VERIFY_IS_EQUAL(ei2.eigenvalues(), ei1.eigenvalues()); + if (rows > 2) { + ei2.compute(a, true, 1); + VERIFY_IS_EQUAL(ei2.info(), NoConvergence); + } + EigenSolver<MatrixType> eiNoEivecs(a, false); VERIFY_IS_EQUAL(eiNoEivecs.info(), Success); VERIFY_IS_APPROX(ei1.eigenvalues(), eiNoEivecs.eigenvalues()); |