From: David Wells 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;