]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Improved breakdown criterion?
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Fri, 6 Oct 2000 10:41:33 +0000 (10:41 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Fri, 6 Oct 2000 10:41:33 +0000 (10:41 +0000)
git-svn-id: https://svn.dealii.org/trunk@3386 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/lac/include/lac/solver_qmrs.h

index 8cf343f4b15ebc41e93813725b3e74c5b0b74025..356a30fe13923702e08e85ab974d729dda5a2716 100644 (file)
@@ -326,10 +326,10 @@ while (state == SolverControl::iterate)
       A.vmult(t,q);
                                       // Step 2
       sigma = q*t;
-
+      
 //TODO: Find a really good breakdown criterion
 // The absolute one detects breakdown instead of convergence
-      if (fabs(sigma) < additional_data.breakdown)
+      if (fabs(sigma/rho) < additional_data.breakdown)
        return breakdown;
                                       // Step 3
       alpha = rho/sigma;

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.