From b7c18230223e9aab0f02ae348bcff6def5c4c275 Mon Sep 17 00:00:00 2001 From: young Date: Wed, 7 Oct 2009 14:05:41 +0000 Subject: [PATCH] Added an assertion to protect when SLEPc fails to converge the number of eigenvalues requested (ie less than requested) git-svn-id: https://svn.dealii.org/trunk@19750 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/slepc_solver.h | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/deal.II/lac/include/lac/slepc_solver.h b/deal.II/lac/include/lac/slepc_solver.h index 5296b113e1..da355171f0 100644 --- a/deal.II/lac/include/lac/slepc_solver.h +++ b/deal.II/lac/include/lac/slepc_solver.h @@ -555,10 +555,10 @@ namespace SLEPcWrappers set_matrices (A,B); solve (n_eigenvectors, &n_converged); - if (n_converged > n_eigenvectors) + if (n_converged >= n_eigenvectors) n_converged = n_eigenvectors; - AssertThrow (n_converged == n_eigenvectors, ExcSLEPcWrappersUsageError()); + AssertThrow (n_converged == n_eigenvectors, ExcSLEPcWrappersUsageError()); AssertThrow (vr.size() >= 1, ExcSLEPcWrappersUsageError()); vr.resize (n_converged, vr.front()); kr.resize (n_converged); -- 2.39.5