Assemble a 2D function

src/one-body-sample.cc

@@ -171,6 +171,22 @@ int main() {
             neumannBoundaryAssembler, b,
             true); // resize the output vector and zero all of its entries
       }
+
+      { // constant 2D function
+        std::vector<SmallVector> b;
+        SmallVector SampleVector;
+        SampleVector[0] = 1;
+        SampleVector[1] = 2;
+        ConstantFunction<SmallVector, SmallVector> fNeumann(SampleVector);
+        NeumannBoundaryAssembler<GridType, SmallVector>
+        neumannBoundaryAssembler(fNeumann);
+
+        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