diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 5f105be4b4947495301d46945cc6908dd0e357c5..3b648758e9b2dd4129aae170d703b89d48c6ab2d 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -171,19 +171,18 @@ int main() { BoundaryPatch<GridType::LeafGridView> neumannBoundary(grid.leafView(), neumannNodes); - { // constant 1D function - std::vector<Dune::FieldVector<double, 1>> nodalIntegrals; - ConstantFunction<SmallVector, Dune::FieldVector<double, 1>> - constantOneFunction(1); - NeumannBoundaryAssembler<GridType, Dune::FieldVector<double, 1>> - neumannBoundaryAssembler(constantOneFunction); - - BoundaryFunctionalAssembler<P1Basis> boundaryFunctionalAssembler( - p1Basis, neumannBoundary); - boundaryFunctionalAssembler.assemble( - neumannBoundaryAssembler, nodalIntegrals, - true); // resize the output vector and zero all of its entries - } + // constant 1D function + std::vector<Dune::FieldVector<double, 1>> nodalIntegrals; + ConstantFunction<SmallVector, Dune::FieldVector<double, 1>> + constantOneFunction(1); + NeumannBoundaryAssembler<GridType, Dune::FieldVector<double, 1>> + neumannBoundaryAssembler(constantOneFunction); + + BoundaryFunctionalAssembler<P1Basis> boundaryFunctionalAssembler( + p1Basis, neumannBoundary); + boundaryFunctionalAssembler.assemble( + neumannBoundaryAssembler, nodalIntegrals, + true); // resize the output vector and zero all of its entries { // constant 2D function std::vector<SmallVector> b; @@ -200,9 +199,8 @@ int main() { neumannBoundaryAssembler, b, true); // resize the output vector and zero all of its entries } - } - { // experiment with convex problems and the like + // experiment with convex problems and the like typedef ZeroNonlinearity<SmallVector, SmallMatrix> NonlinearityType; NonlinearityType phi;