From f9cbf38c30a4ebe578280a772df7b739beda3759 Mon Sep 17 00:00:00 2001 From: kanschat Date: Fri, 3 Dec 1999 02:50:37 +0000 Subject: [PATCH] better breakdown? git-svn-id: https://svn.dealii.org/trunk@1972 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/lac/include/lac/solver_bicgstab.h | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/deal.II/lac/include/lac/solver_bicgstab.h b/deal.II/lac/include/lac/solver_bicgstab.h index 55ea11569a..097b58b144 100644 --- a/deal.II/lac/include/lac/solver_bicgstab.h +++ b/deal.II/lac/include/lac/solver_bicgstab.h @@ -213,9 +213,12 @@ SolverBicgstab::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); -- 2.39.5