From 47dea7acbf15931ccf2d8b59c73eac163f0729d0 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 11 Sep 2012 17:15:34 +0200
Subject: [PATCH] Cleanup
---
src/one-body-sample.cc | 13 +++++++------
1 file changed, 7 insertions(+), 6 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 55b2c72e..d78ffe8a 100644
--- a/src/one-body-sample.cc
+++ b/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"))
--
GitLab