Cleanup

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;