]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
better breakdown?
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 3 Dec 1999 02:50:37 +0000 (02:50 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Fri, 3 Dec 1999 02:50:37 +0000 (02:50 +0000)
git-svn-id: https://svn.dealii.org/trunk@1972 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/solver_bicgstab.h

index 55ea11569a134928c3ea23d07765bbe70a793540..097b58b144e94998ab034735e3d4b2a0df0192c1 100644 (file)
@@ -213,9 +213,12 @@ SolverBicgstab<Matrix, Vector>::iterate(const Preconditioner& precondition)
       MA->vmult(v,y);
       rhobar = rbar * v;
 
-      if (fabs(rhobar) < 1.e-19) return ReturnState(breakdown);
-    
       alpha = rho/rhobar;
+
+//TODO: Find better breakdown criterion (G)
+
+      if (fabs(alpha) > 1.e10) return ReturnState(breakdown);
+    
       s.equ(1., r, -alpha, v);
       precondition(z,s);
       MA->vmult(t,z);

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.