From 6b88ad63317a989fbf66b16131a1300f68fe434b Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 26 Oct 2011 13:26:26 +0200
Subject: [PATCH] Rearrange
---
src/one-body-sample.cc | 9 +++------
1 file changed, 3 insertions(+), 6 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 690a0ed0..94a97492 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -83,13 +83,8 @@ int main() {
for (size_t j = 0; j < dim; ++j)
f[i][j] = i * j;
- // TODO: Why does blockGSStep even provide a default constructor?
- BlockGSStep<OperatorType, VectorType> blockGSStep(stiffnessMatrix, u, f);
- EnergyNorm<OperatorType, VectorType> energyNorm(blockGSStep);
-
Dune::BitSetVector<VectorType::block_type::dimension> ignoreNodes(
grid.size(grid.maxLevel(), dim), false);
-
{ // Play around with the boundary
typedef GridType::LeafGridView GridView;
GridView leafView = grid.leafView();
@@ -136,9 +131,11 @@ int main() {
std::cout << "Number of bounding nodes: " << bounding_nodes << std::endl;
std::cout << "Number of extremal nodes: " << extremal_nodes << std::endl;
}
-
+ // TODO: Why does blockGSStep even provide a default constructor?
+ BlockGSStep<OperatorType, VectorType> blockGSStep(stiffnessMatrix, u, f);
blockGSStep.ignoreNodes_ = &ignoreNodes;
+ EnergyNorm<OperatorType, VectorType> energyNorm(blockGSStep);
LoopSolver<VectorType> solver(&blockGSStep, solver_maxIterations,
solver_tolerance, &energyNorm, Solver::FULL);
--
GitLab