Check reduced residual for equality as well.

git-svn-id: https://svn.dealii.org/trunk@4861 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/lac/source/solver_control.cc b/deal.II/lac/source/solver_control.cc
index 3660cfd05d9e421e693af5105ce817cb159ae1bf..9894feaa8b3dffd58004d53ba1e58def59835829 100644 (file)
--- a/deal.II/lac/source/solver_control.cc
+++ b/deal.II/lac/source/solver_control.cc
@@ -216,13 +216,20 @@ SolverControl::State
 ReductionControl::check (const unsigned int step,
                         const double check_value)
 {
+                                  // if this is the first time we
+                                  // come here, then store the
+                                  // residual for later comparisons
   if (step==0)
     {
       initial_val = check_value;
       reduced_tol = check_value * reduce;
     };
 
-  if (check_value < reduced_tol)
+                                  // check whether desired reduction
+                                  // has been achieved. also check
+                                  // for equality in case initial
+                                  // residual already was zero
+  if (check_value <= reduced_tol)
     {
       if (_log_result)
        deallog << "Convergence step " << step