]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix the problem with constrained dofs in the GMRES solver temporarily by setting...
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 16 Apr 2009 14:13:43 +0000 (14:13 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 16 Apr 2009 14:13:43 +0000 (14:13 +0000)
git-svn-id: https://svn.dealii.org/trunk@18627 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-32/step-32.cc

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

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.