diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 4ccaa3ca5d19e80fcfb450a90cfee412b1e6f445..924be5995da81afe6b958f15697636ed933ecb26 100644 --- a/src/one-body-sample.cc +++ b/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;