diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 489f69ad58063920ead7912edbe655dcaa81e171..b3ca36832d38dea07cf0c1fb4a138a0d3c09b14d 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -244,6 +244,7 @@ int main(int argc, char *argv[]) { ; } } + BoundaryPatch<GridView> const frictionalBoundary(leafView, frictionalNodes); // Set up functions for time-dependent boundary conditions using FunctionType = Dune::VirtualFunction<double, double>; @@ -324,9 +325,6 @@ int main(int argc, char *argv[]) { EnergyNorm<SingletonMatrixType, SingletonVectorType> const stateEnergyNorm( frictionalBoundaryMassMatrix); { - BoundaryPatch<GridView> const frictionalBoundary(leafView, - frictionalNodes); - BoundaryMassAssembler< GridType, BoundaryPatch<GridView>, P1Basis::LocalFiniteElement, P1Basis::LocalFiniteElement, SmallSingletonMatrix> const @@ -580,10 +578,10 @@ int main(int argc, char *argv[]) { relaxationWriter << std::endl; if (parset.get<bool>("io.writeVTK")) { - SingletonVectorType vonMisesStress; auto const gridDisplacement = Dune::make_shared<BasisGridFunction<P1Basis, VectorType> const>( p1Basis, u); + SingletonVectorType vonMisesStress; VonMisesStressAssembler<GridType> localStressAssembler( E, nu, gridDisplacement); p0Assembler.assembleFunctional(localStressAssembler, vonMisesStress);