From f1898414e9df8fd4148e0d39910e7aa62c05619d Mon Sep 17 00:00:00 2001 From: Denis Davydov Date: Fri, 7 Nov 2014 18:24:57 +0100 Subject: [PATCH] switch to EPSGetErrorEstimate as others might give inconsistent results --- source/lac/slepc_solver.cc | 8 +++----- 1 file changed, 3 insertions(+), 5 deletions(-) diff --git a/source/lac/slepc_solver.cc b/source/lac/slepc_solver.cc index 5bfcd6cac1..0aea6ad118 100644 --- a/source/lac/slepc_solver.cc +++ b/source/lac/slepc_solver.cc @@ -211,11 +211,9 @@ namespace SLEPcWrappers // only if at least one eigenvector has converged. if ((*n_converged)>0) { - // EPSComputeRelativeError is consistent with the one - // used during the solution process. Given EPS_CONV_ABS set above, - // this can be either the l2 norm or the mass-matrix induced norm - // when EPS_GHEP is set. - ierr = EPSComputeRelativeError (solver_data->eps, 0, &residual_norm); + // EPSGetErrorEstimate is consistent with the residual norm + // used during the solution process. + ierr = EPSGetErrorEstimate (solver_data->eps, 0, &residual_norm); AssertThrow (ierr == 0, ExcSLEPcError(ierr)); } -- 2.39.5