Fix the problem with constrained dofs in the GMRES solver temporarily by setting...

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 7fa328ce34691b7afb8a24113d8cb43931b6d977..90c7d45fc7cf190dd6612976e595ac103fd59876 100644 (file)
--- 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 << "   "