]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Handle the case of lucky breakdown better, without testing whether 1/a 1062/head
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 29 Jun 2015 21:24:14 +0000 (16:24 -0500)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Mon, 29 Jun 2015 23:11:05 +0000 (18:11 -0500)
is actually finite.

include/deal.II/lac/solver_gmres.h

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.