From: Peter Munch <peterrmuench@gmail.com>
Date: Thu, 19 Nov 2020 23:39:37 +0000 (+0100)
Subject: Fix bug in constructor of FEEvaluationBaseData
X-Git-Tag: v9.3.0-rc1~880^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F11213%2Fhead;p=dealii.git

Fix bug in constructor of FEEvaluationBaseData
---

diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h
index c19f24dcfe..2a2cebe968 100644
--- a/include/deal.II/matrix_free/fe_evaluation.h
+++ b/include/deal.II/matrix_free/fe_evaluation.h
@@ -3266,20 +3266,27 @@ inline FEEvaluationBaseData<dim, Number, is_face, VectorizedArrayType>::
                                      .cell_data[quad_no_in]
                                      .quad_index_from_n_q_points(n_q_points)) :
                         0)
-  , n_quadrature_points(fe_degree != numbers::invalid_unsigned_int ?
-                          n_q_points :
-                          (is_face ? data_in
-                                       .get_shape_info(dof_no,
-                                                       quad_no_in,
-                                                       active_fe_index,
-                                                       active_quad_index)
-                                       .n_q_points_face :
-                                     data_in
-                                       .get_shape_info(dof_no,
-                                                       quad_no_in,
-                                                       active_fe_index,
-                                                       active_quad_index)
-                                       .n_q_points))
+  , n_quadrature_points(
+      fe_degree != numbers::invalid_unsigned_int ?
+        n_q_points :
+        (is_face ? data_in
+                     .get_shape_info(
+                       dof_no,
+                       quad_no_in,
+                       data_in.get_dof_info(dof_no)
+                         .component_to_base_index[first_selected_component],
+                       active_fe_index,
+                       active_quad_index)
+                     .n_q_points_face :
+                   data_in
+                     .get_shape_info(
+                       dof_no,
+                       quad_no_in,
+                       data_in.get_dof_info(dof_no)
+                         .component_to_base_index[first_selected_component],
+                       active_fe_index,
+                       active_quad_index)
+                     .n_q_points))
   , matrix_info(&data_in)
   , dof_info(&data_in.get_dof_info(dof_no))
   , mapping_data(