From f833325be258eb6321dabc62b81cdf2aded404e7 Mon Sep 17 00:00:00 2001
From: Jean-Paul Pelteret <jppelteret@gmail.com>
Date: Sat, 11 May 2019 19:22:15 +0200
Subject: [PATCH] Update indentation of assembly operations

---
 examples/step-44/step-44.cc | 20 ++++++++++----------
 1 file changed, 10 insertions(+), 10 deletions(-)

diff --git a/examples/step-44/step-44.cc b/examples/step-44/step-44.cc
index 202ab92c33..f831cde58e 100644
--- a/examples/step-44/step-44.cc
+++ b/examples/step-44/step-44.cc
@@ -2157,6 +2157,8 @@ namespace Step44
         const SymmetricTensor<4, dim> Jc  = lqph[q_point]->get_Jc();
         const double d2Psi_vol_dJ2        = lqph[q_point]->get_d2Psi_vol_dJ2();
         const double det_F                = lqph[q_point]->get_det_F();
+        const SymmetricTensor<2, dim> &I =
+          Physics::Elasticity::StandardTensors<dim>::I;
 
         // Next we define some aliases to make the assembly process easier to
         // follow
@@ -2185,11 +2187,12 @@ namespace Step44
                 // the local matrix diagonals:
                 if ((i_group == j_group) && (i_group == u_dof))
                   {
-                    data.cell_matrix(i, j) += symm_grad_Nx[i] *
-                                              Jc // The material contribution:
-                                              * symm_grad_Nx[j] * JxW;
-                    if (component_i ==
-                        component_j) // geometrical stress contribution
+                    // The material contribution:
+                    data.cell_matrix(i, j) += symm_grad_Nx[i] * Jc * //
+                                              symm_grad_Nx[j] * JxW; //
+
+                    // The geometrical stress contribution:
+                    if (component_i == component_j)
                       data.cell_matrix(i, j) += grad_Nx[i][component_i] * tau *
                                                 grad_Nx[j][component_j] * JxW;
                   }
@@ -2197,11 +2200,8 @@ namespace Step44
                 // contribution
                 else if ((i_group == p_dof) && (j_group == u_dof))
                   {
-                    data.cell_matrix(i, j) +=
-                      N[i] * det_F *
-                      (symm_grad_Nx[j] *
-                       Physics::Elasticity::StandardTensors<dim>::I) *
-                      JxW;
+                    data.cell_matrix(i, j) += N[i] * det_F *               //
+                                              (symm_grad_Nx[j] * I) * JxW; //
                   }
                 // and lastly the $\mathsf{\mathbf{k}}_{ \widetilde{J}
                 // \widetilde{p}}$ and $\mathsf{\mathbf{k}}_{ \widetilde{J}
-- 
2.39.5