From 5937c5b7fd5cc99717ef096458accaa4a6527d28 Mon Sep 17 00:00:00 2001 From: Elias Pipping <elias.pipping@fu-berlin.de> Date: Wed, 2 Nov 2011 17:44:44 +0100 Subject: [PATCH] Move 1D boundary assembly into a block --- src/one-body-sample.cc | 24 +++++++++++++----------- 1 file changed, 13 insertions(+), 11 deletions(-) diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc index 830939b5..4f7359d5 100644 --- a/src/one-body-sample.cc +++ b/src/one-body-sample.cc @@ -158,17 +158,19 @@ int main() { BoundaryPatch<GridType::LeafGridView> neumannBoundary(grid.leafView(), neumannNodes); - std::vector<Dune::FieldVector<double, 1>> b; - ConstantFunction<SmallVector, Dune::FieldVector<double, 1>> - constantOneFunction(1); - NeumannBoundaryAssembler<GridType, Dune::FieldVector<double, 1>> - neumannBoundaryAssembler(constantOneFunction); - - BoundaryFunctionalAssembler<P1Basis> boundaryFunctionalAssembler( - p1Basis, neumannBoundary); - boundaryFunctionalAssembler.assemble( - neumannBoundaryAssembler, b, - true); // resize the output vector and zero all of its entries + { // constant 1D function + std::vector<Dune::FieldVector<double, 1>> b; + ConstantFunction<SmallVector, Dune::FieldVector<double, 1>> + constantOneFunction(1); + NeumannBoundaryAssembler<GridType, Dune::FieldVector<double, 1>> + neumannBoundaryAssembler(constantOneFunction); + + BoundaryFunctionalAssembler<P1Basis> boundaryFunctionalAssembler( + p1Basis, neumannBoundary); + boundaryFunctionalAssembler.assemble( + neumannBoundaryAssembler, b, + true); // resize the output vector and zero all of its entries + } } { // experiment with convex problems and the like -- GitLab