From: Guido Kanschat Date: Thu, 22 Feb 2001 10:38:02 +0000 (+0000) Subject: residual not used anymore X-Git-Tag: v8.0.0~19698 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=44105d05d0f9c88889c708db68a1822b9fbb216f;p=dealii.git residual not used anymore git-svn-id: https://svn.dealii.org/trunk@4008 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/lac/include/lac/solver_richardson.h b/deal.II/lac/include/lac/solver_richardson.h index 7850c5c403..847f2b4097 100644 --- a/deal.II/lac/include/lac/solver_richardson.h +++ b/deal.II/lac/include/lac/solver_richardson.h @@ -185,7 +185,11 @@ SolverRichardson::solve (const MATRIX &A, // Main loop for(int iter=0; conv==SolverControl::iterate; iter++) { - res=A.residual(r,x,b); + // Do not use residual, + // but do it in 2 steps + A.vmult(r,x); + r.sadd(-1.,1.,b); + res=sqrt(r*r); conv = control().check (iter, criterion()); if (conv != SolverControl::iterate)