From cbde98101223340c6bb9fb2fde9587f39d653d76 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 25 Oct 2011 12:31:44 +0200
Subject: [PATCH] Cleanup
---
src/one-body-sample.cc | 9 ++++-----
1 file changed, 4 insertions(+), 5 deletions(-)
diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index 2ab6ddc9..7c08e4aa 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -40,6 +40,8 @@ int main() {
double const E = 1e8;
double const nu = 0.3;
int const refinements = 5;
+ size_t const solver_maxIterations = 1000;
+ double const solver_tolerance = 1e-5;
typedef Dune::YaspGrid<dim> GridType;
@@ -81,14 +83,11 @@ int main() {
BlockGSStep<OperatorType, VectorType> blockGSStep(stiffnessMatrix, u, f);
EnergyNorm<OperatorType, VectorType> energyNorm(blockGSStep);
- // TODO: Why does blockGSStep not warn about this?
- // Why is this not documented?
Dune::BitSetVector<VectorType::block_type::dimension> ignoreNodes(false);
blockGSStep.ignoreNodes_ = &ignoreNodes;
- LoopSolver<VectorType> solver(&blockGSStep, 1000, // maxIterations
- 1e-5, // tolerance
- &energyNorm, Solver::FULL);
+ LoopSolver<VectorType> solver(&blockGSStep, solver_maxIterations,
+ solver_tolerance, &energyNorm, Solver::FULL);
solver.solve();
return 0;
--
GitLab