]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Make things work.
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 6 Jan 2010 17:07:14 +0000 (17:07 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 6 Jan 2010 17:07:14 +0000 (17:07 +0000)
git-svn-id: https://svn.dealii.org/trunk@20322 0785d39b-7218-0410-832d-ea1e28bc413d

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

index 492afefae080ea4267c72d312c6d6becee7becc9..2d35914611be0fa6a7f03ff6e04b5e250dd2b7cc 100644 (file)
@@ -337,12 +337,12 @@ template <class Preconditioner>
 void SchurComplement<Preconditioner>::vmult (Vector<double>       &dst,
                                             const Vector<double> &src) const
 {
-//  std::cout << " SRC " << src.l2_norm() << std::endl;
   system_matrix->block(0,1).vmult (tmp1, src);
   A_inverse->name = "in schur";
   A_inverse->vmult (tmp2, tmp1);
   system_matrix->block(1,0).vmult (dst, tmp2);
-//  std::cout << " DST " << dst.l2_norm() << std::endl;
+
+  system_matrix->block(1,1).vmult_add (dst, src);
 }
 
 
@@ -669,6 +669,7 @@ void StokesProblem<dim>::assemble_multigrid ()
   mg_A_preconditioner.resize (triangulation.n_levels());
   for (unsigned int level=0; level<triangulation.n_levels(); ++level)
     {
+      Assert (
       mg_A_preconditioner[level]
        = std_cxx1x::shared_ptr<typename InnerPreconditioner<dim>::type>(new typename InnerPreconditioner<dim>::type());
       mg_A_preconditioner[level]

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.