Fix quad_no in MatrixFreeTools::compute_matrix()

diff --git a/include/deal.II/matrix_free/tools.h b/include/deal.II/matrix_free/tools.h
index 92e3f83166ea88537f3292653e6c56f80f6a83e0..48515bbf72e65bdeb1cf5c837cdd52320e75f84b 100644 (file)
--- a/include/deal.II/matrix_free/tools.h
+++ b/include/deal.II/matrix_free/tools.h
@@ -1956,7 +1956,7 @@ namespace MatrixFreeTools
       data_cell;
 
     data_cell.dof_numbers               = {dof_no};
-    data_cell.quad_numbers              = {dof_no};
+    data_cell.quad_numbers              = {quad_no};
     data_cell.first_selected_components = {first_selected_component};
     data_cell.batch_type                = {0};
 
@@ -1992,7 +1992,7 @@ namespace MatrixFreeTools
       data_face;
 
     data_face.dof_numbers               = {dof_no, dof_no};
-    data_face.quad_numbers              = {dof_no, quad_no};
+    data_face.quad_numbers              = {quad_no, quad_no};
     data_face.first_selected_components = {first_selected_component,
                                            first_selected_component};
     data_face.batch_type                = {1, 2};
@@ -2054,7 +2054,7 @@ namespace MatrixFreeTools
       data_boundary;
 
     data_boundary.dof_numbers               = {dof_no};
-    data_boundary.quad_numbers              = {dof_no};
+    data_boundary.quad_numbers              = {quad_no};
     data_boundary.first_selected_components = {first_selected_component};
     data_boundary.batch_type                = {1};