From: David Wells <drwells@email.unc.edu>
Date: Mon, 21 Aug 2023 01:56:24 +0000 (-0400)
Subject: Explicitly spell out a loop.
X-Git-Tag: relicensing~577^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F15895%2Fhead;p=dealii.git

Explicitly spell out a loop.
---

diff --git a/source/fe/fe_values_views_internal.cc b/source/fe/fe_values_views_internal.cc
index ab6abb8a6e..030db4d57a 100644
--- a/source/fe/fe_values_views_internal.cc
+++ b/source/fe/fe_values_views_internal.cc
@@ -355,9 +355,13 @@ namespace FEValuesViews
                 &shape_gradients[snc][0];
               for (unsigned int q_point = 0; q_point < n_quadrature_points;
                    ++q_point)
-                symmetric_gradients[q_point] +=
-                  value * dealii::SymmetricTensor<2, spacedim>(
-                            symmetrize_single_row(comp, *shape_gradient_ptr++));
+                {
+                  for (unsigned int d = 0; d < dim; ++d)
+                    symmetric_gradients[q_point][comp][d] +=
+                      0.5 * value * (*shape_gradient_ptr)[d];
+                  symmetric_gradients[q_point][comp][comp] +=
+                    0.5 * value * (*shape_gradient_ptr++)[comp];
+                }
             }
           else
             for (unsigned int q_point = 0; q_point < n_quadrature_points;