[ContactStructuralMechanicsApplication] Minor optimization in `mesh_t…

…ying_mortar_condition.cpp`
KratosMultiphysics · Aug 25, 2023 · 1052fd0 · 1052fd0
1 parent 7e509be
commit 1052fd0
Showing 1 changed file with 8 additions and 8 deletions.
diff --git a/...s/ContactStructuralMechanicsApplication/custom_conditions/mesh_tying_mortar_condition.cpp b/...s/ContactStructuralMechanicsApplication/custom_conditions/mesh_tying_mortar_condition.cpp
@@ -550,10 +550,10 @@ void MeshTyingMortarCondition<TDim,TNumNodes, TNumNodesMaster>::CalculateLocalLH
     const IndexType initial_column_index = dof_size * (TNumNodes + TNumNodesMaster);
 
     // Iterate over the number of dofs on master side
-    for (IndexType i = 0; i < dof_size; ++i) {
-        for (IndexType j = 0; j < TNumNodesMaster; ++j) {
-            for (IndexType k = 0; k < TNumNodes; ++k) {
-                const double value = - scale_factor * r_MOperator(k, j);
+    for (IndexType j = 0; j < TNumNodesMaster; ++j) {
+        for (IndexType k = 0; k < TNumNodes; ++k) {
+            const double value = - scale_factor * r_MOperator(k, j);
+            for (IndexType i = 0; i < dof_size; ++i) {
                 rLocalLHS(initial_row_index + j * dof_size + i, initial_column_index + k * dof_size + i) = value;
                 rLocalLHS(initial_column_index + k * dof_size + i, initial_row_index + j * dof_size + i) = value;
             }
@@ -564,10 +564,10 @@ void MeshTyingMortarCondition<TDim,TNumNodes, TNumNodesMaster>::CalculateLocalLH
     initial_row_index = dof_size * TNumNodesMaster;
 
     // Iterate over the number of dofs on slave side
-    for (IndexType i = 0; i < dof_size; ++i) {
-        for (IndexType j = 0; j < TNumNodes; ++j) {
-            for (IndexType k = 0; k < TNumNodes; ++k) {
-                const double value = scale_factor * r_DOperator(k, j);
+    for (IndexType j = 0; j < TNumNodes; ++j) {
+        for (IndexType k = 0; k < TNumNodes; ++k) {
+            const double value = scale_factor * r_DOperator(k, j);
+            for (IndexType i = 0; i < dof_size; ++i) {
                 rLocalLHS(initial_row_index + j * dof_size + i, initial_column_index + k * dof_size + i) = value;
                 rLocalLHS(initial_column_index + k * dof_size + i, initial_row_index + j * dof_size + i) = value;
             }