From: Baerbel Jannsen <baerbel.janssen@gmail.com>
Date: Fri, 12 Mar 2010 00:09:56 +0000 (+0000)
Subject: add output for debugging
X-Git-Tag: v8.0.0~6339
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=07577ce664f4489ab396ded0b1c2f96990331279;p=dealii.git

add output for debugging


git-svn-id: https://svn.dealii.org/trunk@20803 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/deal.II/include/multigrid/multigrid.h b/deal.II/deal.II/include/multigrid/multigrid.h
index 097526fcd9..8442d467ef 100644
--- a/deal.II/deal.II/include/multigrid/multigrid.h
+++ b/deal.II/deal.II/include/multigrid/multigrid.h
@@ -664,14 +664,30 @@ PreconditionMG<dim, VECTOR, TRANSFER>::vmult (
   VECTOR2& dst,
   const VECTOR2& src) const
 {
+#ifdef baerbel_mg_test
+  const unsigned int min_l = multigrid->solution.get_minlevel();
+  const unsigned int max_l = multigrid->solution.get_maxlevel();
+  deallog << "PreconditionMG in " << src.l2_norm() << std::endl;
+#endif
   transfer->copy_to_mg(*mg_dof_handler,
 		       multigrid->defect,
 		       src);
+#ifdef baerbel_mg_test
+  for(unsigned int l=min_l; l<max_l; ++l)
+  deallog << "defect [" << l << "] " << multigrid->defect[l].l2_norm() << std::endl;
+#endif
   multigrid->cycle();
+
+#ifdef baerbel_mg_test
+  for(unsigned int l=min_l; l<max_l; ++l)
+  deallog << "solution [" << l << "] " << multigrid->solution[l].l2_norm() << std::endl;
+#endif
   transfer->copy_from_mg(*mg_dof_handler,
 			 dst,
 			 multigrid->solution);
-
+#ifdef baerbel_mg_test
+  deallog << "PreconditionMG out " << dst.l2_norm() << std::endl;
+#endif
   /*
   DataOut<dim> data_out;
   data_out.attach_dof_handler (*mg_dof_handler);