diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 46b763bcd4dfd9458bdee1530691c93be080b7a3..4ca32ba6d000afb7c50eac4338e88677fa8cdf6e 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -226,6 +226,7 @@ int main() { coefficientOfFriction, normalStress, nodalIntegrals); // }}} + EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix); { typedef MyConvexProblem<OperatorType, VectorType, Dune::LinearFunction> MyConvexProblemType; @@ -240,8 +241,6 @@ int main() { nonlinearGSStep.ignoreNodes_ = &ignoreNodes; nonlinearGSStep.setProblem(u1, myBlockProblem); - // FIXME: Does this make any sense? - EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix); LoopSolver<VectorType> solver(&nonlinearGSStep, solver_maxIterations, solver_tolerance, &energyNorm, // Solver::QUIET); @@ -256,8 +255,6 @@ int main() { BlockGSStep<OperatorType, VectorType> blockGSStep(stiffnessMatrix, u2, f); blockGSStep.ignoreNodes_ = &ignoreNodes; - // FIXME: Does this make any sense? - EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix); LoopSolver<VectorType> solver(&blockGSStep, solver_maxIterations, solver_tolerance, &energyNorm, // Solver::QUIET);