From c27d741a23d57fb32f86555f537197180cb3bd1a Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Tue, 19 Aug 2008 22:26:20 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-31/step-31.cc | 1 + 1 file changed, 1 insertion(+) 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 << " " -- 2.39.5