Move energyNorm

src/one-body-sample.cc

@@ -160,6 +160,7 @@ int main() {
     OperatorType stiffnessMatrix;
     OperatorAssembler<P1Basis, P1Basis>(p1Basis, p1Basis)
         .assemble(localStiffness, stiffnessMatrix);
+    EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix);
 
     Dune::BitSetVector<dim> ignoreNodes(grid.size(grid.maxLevel(), dim), false);
     Dune::BitSetVector<1> neumannNodes(grid.size(grid.maxLevel(), dim), false);
@@ -198,7 +199,6 @@ int main() {
         coefficientOfFriction, normalStress, nodalIntegrals);
     // }}}
 
-    EnergyNorm<OperatorType, VectorType> energyNorm(stiffnessMatrix);
     {
       typedef MyConvexProblem<OperatorType, VectorType, Dune::LinearFunction>
       MyConvexProblemType;