From: Guido Kanschat Date: Fri, 6 Oct 2000 10:41:33 +0000 (+0000) Subject: Improved breakdown criterion? X-Git-Tag: v8.0.0~20045 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=0f915045a1d64f904ea8a4a0c82ada1b65e947c6;p=dealii.git Improved breakdown criterion? git-svn-id: https://svn.dealii.org/trunk@3386 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/solver_qmrs.h b/deal.II/lac/include/lac/solver_qmrs.h index 8cf343f4b1..356a30fe13 100644 --- a/deal.II/lac/include/lac/solver_qmrs.h +++ b/deal.II/lac/include/lac/solver_qmrs.h @@ -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;