b -> nodalIntegrals

src/one-body-sample.cc

@@ -169,7 +169,7 @@ int main() {
                                                             neumannNodes);
 
       { // constant 1D function
-        std::vector<Dune::FieldVector<double, 1>> b;
+        std::vector<Dune::FieldVector<double, 1>> nodalIntegrals;
         ConstantFunction<SmallVector, Dune::FieldVector<double, 1>>
         constantOneFunction(1);
         NeumannBoundaryAssembler<GridType, Dune::FieldVector<double, 1>>
@@ -178,7 +178,7 @@ int main() {
         BoundaryFunctionalAssembler<P1Basis> boundaryFunctionalAssembler(
             p1Basis, neumannBoundary);
         boundaryFunctionalAssembler.assemble(
-            neumannBoundaryAssembler, b,
+            neumannBoundaryAssembler, nodalIntegrals,
             true); // resize the output vector and zero all of its entries
       }