From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Mon, 29 Jun 2015 21:24:14 +0000 (-0500)
Subject: Handle the case of lucky breakdown better, without testing whether 1/a
X-Git-Tag: v8.3.0-rc1~77^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F1062%2Fhead;p=dealii.git

Handle the case of lucky breakdown better, without testing whether 1/a
is actually finite.
---

diff --git a/include/deal.II/lac/solver_gmres.h b/include/deal.II/lac/solver_gmres.h
index 0e5c6af011..ce2bdf292a 100644
--- a/include/deal.II/lac/solver_gmres.h
+++ b/include/deal.II/lac/solver_gmres.h
@@ -977,7 +977,7 @@ SolverFGMRES<VECTOR>::solve (
 
       for (unsigned int j=0; j<basis_size; ++j)
         {
-          if (numbers::is_finite(1./a)) // treat lucky breakdown
+          if (a != 0) // treat lucky breakdown
             v(j,x).equ(1./a, *aux);
           else
             v(j,x) = 0.;