]> https://gitweb.dealii.org/ - dealii.git/commitdiff
initialize the MGLevelObject correctly
authorBaerbel Jannsen <baerbel.janssen@gmail.com>
Tue, 9 Feb 2010 10:26:50 +0000 (10:26 +0000)
committerBaerbel Jannsen <baerbel.janssen@gmail.com>
Tue, 9 Feb 2010 10:26:50 +0000 (10:26 +0000)
git-svn-id: https://svn.dealii.org/trunk@20525 0785d39b-7218-0410-832d-ea1e28bc413d

tests/multigrid/transfer_select.cc

index 34403bf805332442c14627d7fb4725fc77d6b224..aeb3c20c06300181c8e8bdd7139a184c7f77b1f1 100644 (file)
@@ -63,6 +63,7 @@ void check_select(const FiniteElement<dim>& fe,
   
   for (unsigned int l=0;l<tr.n_levels();++l)
     DoFRenumbering::component_wise(mgdof, l, mg_target_component);
+
   std::vector<std::vector<unsigned int> > mg_ndofs(mgdof.get_tria().n_levels());
   MGTools::count_dofs_per_component(mgdof, mg_ndofs, true, mg_target_component);
 
@@ -106,9 +107,10 @@ void check_select(const FiniteElement<dim>& fe,
   for (unsigned int i=0;i<u.size();++i)
     u(i) = i;
   
-  MGLevelObject<Vector<double> > v;
-  v.resize(2,2);
-  v[2].reinit(mg_ndofs[2][mg_selected]);
+  MGLevelObject<Vector<double> > v(0,tr.n_levels()-1);
+  for(unsigned int l=0; l<tr.n_levels()-1; ++l)
+    v[l].reinit(mg_ndofs[l][mg_target_component[mg_selected]]);
+
 
   transfer.copy_to_mg(mgdof, v, u);
   for (unsigned int i=0; i<v[2].size();++i)

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.