diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 2ab6ddc9eee2ec47aa310a1893859a0470371deb..7c08e4aa1b060670ceadf01d0fac6646647ab10b 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;