From: guido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Wed, 1 Sep 2004 17:04:56 +0000 (+0000)
Subject: Wolfgang's changes reverted; were not helpful and need further iscussion
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ca3a757d34a5e757cf5069554e63e112fa6a740c;p=dealii-svn.git

Wolfgang's changes reverted; were not helpful and need further iscussion


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

diff --git a/deal.II/deal.II/include/multigrid/mg_dof_tools.h b/deal.II/deal.II/include/multigrid/mg_dof_tools.h
index 7c2a46105a..5b1271df7c 100644
--- a/deal.II/deal.II/include/multigrid/mg_dof_tools.h
+++ b/deal.II/deal.II/include/multigrid/mg_dof_tools.h
@@ -185,7 +185,6 @@ class MGTools
 		     const std::vector<bool>& selected = std::vector<bool>(),
 		     const std::vector<unsigned int>& target_component
 		     = std::vector<unsigned int>());
-
 				     /**
 				      * Adjust vectors on all levels
 				      * to correct size.  Count the
@@ -208,7 +207,7 @@ class MGTools
 				      *
 				      * The argument
 				      * @p target_component allows to
-				      * re-sort and group components
+				      * re-sort and groupt components
 				      * as in
 				      * DoFRenumbering::component_wise.
 				      */
@@ -218,7 +217,7 @@ class MGTools
 		     MGLevelObject<Vector<number> > &v,
 		     const std::vector<bool> &selected,
 		     const std::vector<unsigned int> &target_component,
-		     const std::vector<std::vector<unsigned int> >& cached_sizes
+		     std::vector<std::vector<unsigned int> >& cached_sizes
 		     = std::vector<std::vector<unsigned int> >());
 };
 
diff --git a/deal.II/deal.II/source/multigrid/mg_dof_tools.cc b/deal.II/deal.II/source/multigrid/mg_dof_tools.cc
index 2ce07b5122..268613979e 100644
--- a/deal.II/deal.II/source/multigrid/mg_dof_tools.cc
+++ b/deal.II/deal.II/source/multigrid/mg_dof_tools.cc
@@ -20,7 +20,6 @@
 #include <lac/compressed_sparsity_pattern.h>
 #include <lac/sparsity_pattern.h>
 #include <lac/block_vector.h>
-#include <lac/vector.h>
 #include <grid/tria.h>
 #include <grid/tria_iterator.h>
 #include <multigrid/mg_dof_handler.h>
@@ -465,7 +464,7 @@ MGTools::reinit_vector (const MGDoFHandler<dim> &mg_dof,
                         MGLevelObject<Vector<number> > &v,
                         const std::vector<bool> &selected,
 			const std::vector<unsigned int> &target_component,
-			const std::vector<std::vector<unsigned int> >& cached_dofs)
+			std::vector<std::vector<unsigned int> >& ndofs)
 {
   Assert (selected.size() == target_component.size(),
 	  ExcDimensionMismatch(selected.size(), target_component.size()));
@@ -483,19 +482,19 @@ MGTools::reinit_vector (const MGDoFHandler<dim> &mg_dof,
   while (!selected[selected_block])
     ++selected_block;
 
-  std::vector<std::vector<unsigned int> >
-    n_dofs(mg_dof.get_tria().n_levels(),
-	   std::vector<unsigned int>(target_component.size()));
-
-  if (cached_dofs.size() == 0)
-    count_dofs_per_component (mg_dof, n_dofs, target_component);
-  else
-    n_dofs = cached_dofs;
+  if (ndofs.size() == 0)
+    {
+      std::vector<std::vector<unsigned int> >
+	new_dofs(mg_dof.get_tria().n_levels(),
+		 std::vector<unsigned int>(target_component.size()));
+      std::swap(ndofs, new_dofs);
+      count_dofs_per_component (mg_dof, ndofs, target_component);
+    }
   
   for (unsigned int level=v.get_minlevel();
        level<=v.get_maxlevel();++level)
     {
-      v[level].reinit(n_dofs[level][selected_block]);
+      v[level].reinit(ndofs[level][selected_block]);
     }
 }
 
@@ -617,13 +616,13 @@ template void MGTools::reinit_vector<deal_II_dimension> (
   MGLevelObject<Vector<double> >&,
   const std::vector<bool>&,
   const std::vector<unsigned int>&,
-  const std::vector<std::vector<unsigned int> >&);
+  std::vector<std::vector<unsigned int> >&);
 template void MGTools::reinit_vector<deal_II_dimension> (
   const MGDoFHandler<deal_II_dimension>&,
   MGLevelObject<Vector<float> >&,
   const std::vector<bool>&,
   const std::vector<unsigned int>&,
-  const std::vector<std::vector<unsigned int> >&);
+  std::vector<std::vector<unsigned int> >&);
 
 
 template void MGTools::count_dofs_per_component<deal_II_dimension> (