From 9f61d074b5b85cb53b1e575e97edeab8968f3cdb Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 21 Mar 2024 15:07:44 +0100 Subject: [PATCH] Switch indices --- include/deal.II/non_matching/mapping_info.h | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/include/deal.II/non_matching/mapping_info.h b/include/deal.II/non_matching/mapping_info.h index 1742fd46d3..85b1864bc8 100644 --- a/include/deal.II/non_matching/mapping_info.h +++ b/include/deal.II/non_matching/mapping_info.h @@ -1960,13 +1960,13 @@ namespace NonMatching v) = mapping_data.jacobians[q * n_lanes + v][d][s]; if (update_flags_mapping & UpdateFlags::update_inverse_jacobians) - for (unsigned int d = 0; d < spacedim; ++d) - for (unsigned int s = 0; s < dim; ++s) + for (unsigned int d = 0; d < dim; ++d) + for (unsigned int s = 0; s < spacedim; ++s) dealii::internal::VectorizedArrayTrait::get( inverse_jacobians[is_interior ? 0 : 1] - [compressed_offset][s][d], + [compressed_offset][d][s], v) = - mapping_data.inverse_jacobians[q * n_lanes + v][s][d]; + mapping_data.inverse_jacobians[q * n_lanes + v][d][s]; } if (is_interior) -- 2.39.5