Merge pull request #12992 from Rombur/mf_reinit_bug

Fix bug MatrixFree::reinit
dealii · Nov 25, 2021 · a0d9f66 · a0d9f66
2 parents d030bc0 + fa55227
commit a0d9f66
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/include/deal.II/matrix_free/matrix_free.h b/include/deal.II/matrix_free/matrix_free.h
@@ -3098,7 +3098,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
   std::vector<const DoFHandler<dim> *> dof_handlers;
 
   for (const auto dh : dof_handler)
-    dof_handlers.push_back(dof_handler);
+    dof_handlers.push_back(dh);
 
   this->reinit(dof_handlers, constraint, quad, additional_data);
 }
@@ -3179,7 +3179,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
   std::vector<const DoFHandler<dim> *> dof_handlers;
 
   for (const auto dh : dof_handler)
-    dof_handlers.push_back(dof_handler);
+    dof_handlers.push_back(dh);
 
   this->reinit(mapping, dof_handlers, constraint, quad, additional_data);
 }
@@ -3201,7 +3201,7 @@ MatrixFree<dim, Number, VectorizedArrayType>::reinit(
   std::vector<const DoFHandler<dim> *> dof_handlers;
 
   for (const auto dh : dof_handler)
-    dof_handlers.push_back(dof_handler);
+    dof_handlers.push_back(dh);
 
   this->reinit(dof_handlers, constraint, quad, additional_data);
 }