From aea9c126f62b838becae32193884ab702fdd3869 Mon Sep 17 00:00:00 2001 From: Peter Munch Date: Fri, 20 Nov 2020 00:39:37 +0100 Subject: [PATCH] Fix bug in constructor of FEEvaluationBaseData --- include/deal.II/matrix_free/fe_evaluation.h | 35 ++++++++++++--------- 1 file changed, 21 insertions(+), 14 deletions(-) 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:: .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( -- 2.39.5