From 5feeabd16b7152fc64bc67b77921d1ad31fe423c Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Thu, 3 Nov 2011 14:19:03 +0100
Subject: [PATCH] Clean up/comment nodal basis and elastic assembly
---
src/one-body-sample.cc | 12 +++++-------
1 file changed, 5 insertions(+), 7 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 562c4cb7..3501cda8 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -162,20 +162,17 @@ int main() {
0); // zero overlap (whatever that is)
grid.globalRefine(refinements);
+ // Set up nodal basis
typedef P1NodalBasis<GridType::LeafGridView, double> P1Basis;
P1Basis const p1Basis(grid.leafView());
- OperatorAssembler<P1Basis, P1Basis> const globalAssembler(p1Basis, p1Basis);
-
+ // Assemble elastic force on the body
StVenantKirchhoffAssembler<GridType, P1Basis::LocalFiniteElement,
P1Basis::LocalFiniteElement> const
localStiffness(E, nu);
-
OperatorType stiffnessMatrix;
- globalAssembler.assemble(localStiffness, stiffnessMatrix);
-
- // Filled in some block
- VectorType f(grid.size(grid.maxLevel(), dim));
+ OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis)
+ .assemble(localStiffness, stiffnessMatrix);
Dune::BitSetVector<dim> ignoreNodes(grid.size(grid.maxLevel(), dim), false);
Dune::BitSetVector<1> neumannNodes(grid.size(grid.maxLevel(), dim), false);
@@ -198,6 +195,7 @@ int main() {
normalStress.resize(grid.size(grid.maxLevel(), dim));
std::fill(normalStress.begin(), normalStress.end(), 0.0);
+ VectorType f(grid.size(grid.maxLevel(), dim));
assemble_neumann<GridType, SmallVector, P1Basis>(grid, p1Basis,
neumannNodes, f);
--
GitLab