From 4f2cd011eed684c14eafb8dd47df6e98026df28f Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Wed, 12 Feb 2020 15:46:25 +0100 Subject: [PATCH] Ensure to have valid quadrature point data for all SIMD lanes --- include/deal.II/matrix_free/mapping_info.templates.h | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/include/deal.II/matrix_free/mapping_info.templates.h b/include/deal.II/matrix_free/mapping_info.templates.h index 307abf4348..7c02c202c4 100644 --- a/include/deal.II/matrix_free/mapping_info.templates.h +++ b/include/deal.II/matrix_free/mapping_info.templates.h @@ -1457,6 +1457,11 @@ namespace internal face_data.normal_vectors[q][d][v] = face_data.normal_vectors[q][d][0]; } + if (fe_face_values.get_update_flags() & + update_quadrature_points) + for (unsigned int d = 0; d < dim; ++d) + face_data.quadrature_points[q][d][v] = + face_data.quadrature_points[q][d][0]; } } if (is_boundary_face == false && -- 2.39.5