Skip to content
Snippets Groups Projects
Commit 5a4a5567 authored by Elias Pipping's avatar Elias Pipping Committed by Elias Pipping
Browse files

Move blocks around

parent 0f650ca6
No related branches found
No related tags found
No related merge requests found
......@@ -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;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment