From 89d4713a181b1b5bde43f846c0a85c3d71417d6d Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 29 Oct 2020 12:17:39 +0100 Subject: [PATCH] Switch SymmetricTensor -> Tensor for better performance. --- include/deal.II/fe/mapping_q_internal.h | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/include/deal.II/fe/mapping_q_internal.h b/include/deal.II/fe/mapping_q_internal.h index 9ea8c9c231..bb3c28ca59 100644 --- a/include/deal.II/fe/mapping_q_internal.h +++ b/include/deal.II/fe/mapping_q_internal.h @@ -912,7 +912,7 @@ namespace internal return; } - SymmetricTensor<2, n_functions> matrix; + Tensor<2, n_functions> matrix; std::array shape_values; for (unsigned int q = 0; q < unit_support_points.size(); ++q) { @@ -935,7 +935,7 @@ namespace internal // reference point positions which sets up an inverse // interpolation. for (unsigned int i = 0; i < n_functions; ++i) - for (unsigned int j = 0; j <= i; ++j) + for (unsigned int j = 0; j < n_functions; ++j) matrix[i][j] += shape_values[i] * shape_values[j]; for (unsigned int i = 0; i < n_functions; ++i) coefficients[i] += shape_values[i] * unit_support_points[q]; -- 2.39.5