From: Martin Kronbichler <martin.kronbichler@rub.de>
Date: Sun, 5 May 2024 10:02:07 +0000 (+0200)
Subject: Fix FEEvaluation gradient layout in step-76 tutorial program
X-Git-Tag: v9.6.0-rc1~312^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F16966%2Fhead;p=dealii.git

Fix FEEvaluation gradient layout in step-76 tutorial program
---

diff --git a/examples/step-76/step-76.cc b/examples/step-76/step-76.cc
index 1d9fdb1a99..c872469b27 100644
--- a/examples/step-76/step-76.cc
+++ b/examples/step-76/step-76.cc
@@ -741,17 +741,19 @@ namespace Euler_DG
               for (unsigned int c = 0; c < dim + 2; ++c)
                 {
                   if (dim >= 1 && body_force.get() == nullptr)
-                    eval.template gradients<0, false, false>(
-                      gradient_ptr + phi.static_n_q_points * 0, values_ptr);
+                    eval.template gradients<0, false, false, dim>(gradient_ptr,
+                                                                  values_ptr);
                   else if (dim >= 1)
-                    eval.template gradients<0, false, true>(
-                      gradient_ptr + phi.static_n_q_points * 0, values_ptr);
+                    eval.template gradients<0, false, true, dim>(gradient_ptr,
+                                                                 values_ptr);
                   if (dim >= 2)
-                    eval.template gradients<1, false, true>(
-                      gradient_ptr + phi.static_n_q_points * 1, values_ptr);
+                    eval.template gradients<1, false, true, dim>(gradient_ptr +
+                                                                   1,
+                                                                 values_ptr);
                   if (dim >= 3)
-                    eval.template gradients<2, false, true>(
-                      gradient_ptr + phi.static_n_q_points * 2, values_ptr);
+                    eval.template gradients<2, false, true, dim>(gradient_ptr +
+                                                                   2,
+                                                                 values_ptr);
 
                   values_ptr += phi.static_n_q_points;
                   gradient_ptr += phi.static_n_q_points * dim;