From: kronbichler Date: Tue, 19 Aug 2008 22:26:20 +0000 (+0000) Subject: At the start of each Stokes solve, initialize solution vector to zero. This drives... X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=47359315b72d6f27339004df84c808fa959da2dc;p=dealii-svn.git At the start of each Stokes solve, initialize solution vector to zero. This drives down the number of iterations by a full one third :-) git-svn-id: https://svn.dealii.org/trunk@16598 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-31/step-31.cc b/deal.II/examples/step-31/step-31.cc index 37ccb5b971..bd8e6d2883 100644 --- a/deal.II/examples/step-31/step-31.cc +++ b/deal.II/examples/step-31/step-31.cc @@ -2104,6 +2104,7 @@ void BoussinesqFlowProblem::solve () SolverGMRES > gmres(solver_control, SolverGMRES >::AdditionalData(100)); + stokes_solution = 0; gmres.solve(stokes_matrix, stokes_solution, stokes_rhs, preconditioner); std::cout << " "