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 0e5c6af01189e6f82f5f5b935e4d328adbd67574..ce2bdf292a608d54505c61c1de1a3b7a9f8fbdeb 100644 (file)
--- 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.;