Fix bug MatrixFree::reinit

diff --git a/include/deal.II/matrix_free/matrix_free.h b/include/deal.II/matrix_free/matrix_free.h
index c63110990ec579948ae4b8876b68d07db85bb2be..54cc39e52613cb0edc2f998159eda140ac83ea81 100644 (file)
--- a/include/deal.II/matrix_free/matrix_free.h
+++ b/include/deal.II/matrix_free/matrix_free.h
@@ -3087,7 +3087,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);
 }
@@ -3168,7 +3168,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);
 }
@@ -3190,7 +3190,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);
 }