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

Do not write to f directly; add to it

parent 1b9f9b0d
No related branches found
No related tags found
No related merge requests found
......@@ -173,8 +173,11 @@ int main() {
VectorType u2 = u1;
VectorType f(grid.size(grid.maxLevel(), dim));
f = 0;
VectorType neumannTerm(grid.size(grid.maxLevel(), dim));
assemble_neumann<GridType, GridType::LeafGridView, SmallVector, P1Basis>(
leafView, p1Basis, neumannNodes, f);
leafView, p1Basis, neumannNodes, neumannTerm);
f += neumannTerm;
// {{{ Assemble terms for the nonlinearity
std::vector<Dune::FieldVector<double, 1>> nodalIntegrals;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment