From e22c9e084fe5d3edb92abf1d6ebcf9f186e4f3e5 Mon Sep 17 00:00:00 2001 From: kronbichler Date: Thu, 16 Apr 2009 14:13:43 +0000 Subject: [PATCH] Fix the problem with constrained dofs in the GMRES solver temporarily by setting the whole vector to zero. However, we should use something like in step-31 eventually. git-svn-id: https://svn.dealii.org/trunk@18627 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-32/step-32.cc | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index 7fa328ce34..90c7d45fc7 100644 --- a/deal.II/examples/step-32/step-32.cc +++ b/deal.II/examples/step-32/step-32.cc @@ -970,7 +970,7 @@ void BoussinesqFlowProblem::assemble_stokes_system () (Point (0,0,1)) ); for (unsigned int i=0; i::solve () gmres (solver_control, SolverGMRES::AdditionalData(100)); + // currently, have a problem with + // constrained dofs that are nonzero when + // entering the solver (see step-31). set + // the vector to zero for the moment. + stokes_solution = 0; gmres.solve(stokes_matrix, stokes_solution, stokes_rhs, preconditioner); pcout << " " -- 2.39.5