From: Martin Kronbichler <martin.kronbichler@rub.de>
Date: Tue, 26 Mar 2024 16:39:48 +0000 (+0100)
Subject: FEEvaluation: Fix submit_gradient in 1D
X-Git-Tag: v9.6.0-rc1~446^2~2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=88a6aa475c757305a2fdf2adba05f6ac6110a09d;p=dealii.git

FEEvaluation: Fix submit_gradient in 1D
---

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<dim, n_components_, Number, is_face, VectorizedArrayType>::
   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<dim, n_components_, Number, is_face, VectorizedArrayType>::
   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 &&