Move energyNorm

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);