From 47359315b72d6f27339004df84c808fa959da2dc Mon Sep 17 00:00:00 2001 From: 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