From: kronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 16 Apr 2009 14:13:43 +0000 (+0000)
Subject: Fix the problem with constrained dofs in the GMRES solver temporarily by setting... 
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=e22c9e084fe5d3edb92abf1d6ebcf9f186e4f3e5;p=dealii-svn.git

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
---

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<dim>::assemble_stokes_system ()
 					 (Point<dim> (0,0,1)) );
 	    for (unsigned int i=0; i<dofs_per_cell; ++i)
 	      local_rhs(i) += (EquationData::Rayleigh_number *
-			      gravity * phi_u[i] * old_temperature)*
+			       gravity * phi_u[i] * old_temperature)*
 			      stokes_fe_values.JxW(q);
 	  }
   
@@ -1376,6 +1376,11 @@ void BoussinesqFlowProblem<dim>::solve ()
       gmres (solver_control,
 	     SolverGMRES<TrilinosWrappers::MPI::BlockVector >::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 << "   "