treat tke lucky breakdown, otherwise we had devision by zero

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

diff --git a/deal.II/include/deal.II/lac/solver_gmres.h b/deal.II/include/deal.II/lac/solver_gmres.h
index 3639ae9d9fdde7999152a7fbc972c1a2da255d17..2ca953ca831d701365b3a0226aedb645c4ae9a60 100644 (file)
--- a/deal.II/include/deal.II/lac/solver_gmres.h
+++ b/deal.II/include/deal.II/lac/solver_gmres.h
@@ -706,9 +706,9 @@ SolverGMRES<VECTOR>::solve (const MATRIX         &A,
 
           const double s = vv.l2_norm();
           h(inner_iteration+1) = s;
-//TODO: s=0 is a lucky breakdown. Handle this somehow decently
-
-          vv *= 1./s;
+          //s=0 is a lucky breakdown. Handle this somehow decently
+          if(s!=0)
+            vv *= 1./s;
 
                                            /*  Transformation into
                                                triagonal structure  */