From: Denis Davydov <davydden@gmail.com>
Date: Fri, 7 Nov 2014 17:24:57 +0000 (+0100)
Subject: switch to EPSGetErrorEstimate as others might give inconsistent results
X-Git-Tag: v8.2.0-rc1~35^2~3
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f1898414e9df8fd4148e0d39910e7aa62c05619d;p=dealii.git

switch to EPSGetErrorEstimate as others might give inconsistent results
---

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));
         }