From 88a6aa475c757305a2fdf2adba05f6ac6110a09d Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Tue, 26 Mar 2024 17:39:48 +0100 Subject: [PATCH] FEEvaluation: Fix submit_gradient in 1D --- include/deal.II/matrix_free/fe_evaluation.h | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/include/deal.II/matrix_free/fe_evaluation.h b/include/deal.II/matrix_free/fe_evaluation.h index 83d3eb462e..857927cd79 100644 --- a/include/deal.II/matrix_free/fe_evaluation.h +++ b/include/deal.II/matrix_free/fe_evaluation.h @@ -5149,7 +5149,7 @@ FEEvaluationBase:: submit_gradient(const Tensor<2, 1, VectorizedArrayType> grad_in, const unsigned int q_point) { - submit_value(grad_in[0], q_point); + submit_gradient(grad_in[0], q_point); } @@ -5434,8 +5434,9 @@ FEEvaluationBase:: VectorizedArrayType divergence; const std::size_t nqp = this->n_quadrature_points; - if (this->data->element_type == - internal::MatrixFreeFunctions::ElementType::tensor_raviart_thomas) + if (dim > 1 && + this->data->element_type == + internal::MatrixFreeFunctions::ElementType::tensor_raviart_thomas) { VectorizedArrayType inv_det = (!is_face && -- 2.39.5