Cleanup

src/one-body-sample.cc

@@ -166,11 +166,13 @@ int main(int argc, char *argv[]) {
     P0Basis const p0Basis(leafView);
     P1Basis const p1Basis(leafView);
 
-    MassAssembler<GridType, P1Basis::LocalFiniteElement,
-                  P1Basis::LocalFiniteElement> const localMass;
     MatrixType massMatrix;
     {
       timer.reset();
+
+      // Assemble mass matrix and compute density
+      MassAssembler<GridType, P1Basis::LocalFiniteElement,
+                    P1Basis::LocalFiniteElement> const localMass;
       OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis)
           .assemble(localMass, massMatrix);
 
@@ -186,7 +188,6 @@ int main(int argc, char *argv[]) {
           area *= (upperRight[i] - lowerLeft[i]);
 
       double const gravity = 9.81;
-
       double const normalStress =
           parset.get<double>("boundary.friction.normalstress");
 
@@ -201,12 +202,12 @@ int main(int argc, char *argv[]) {
     }
 
     // Assemble elastic force on the body
-    StVenantKirchhoffAssembler<GridType, P1Basis::LocalFiniteElement,
-                               P1Basis::LocalFiniteElement> const
-    localStiffness(E, nu);
     MatrixType stiffnessMatrix;
     {
       timer.reset();
+      StVenantKirchhoffAssembler<GridType, P1Basis::LocalFiniteElement,
+                                 P1Basis::LocalFiniteElement> const
+      localStiffness(E, nu);
       OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis)
           .assemble(localStiffness, stiffnessMatrix);
       if (parset.get<bool>("enable_timer"))