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

Use umv for performance

parent 607e9fb8
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 4 additions and 7 deletions
src/one-body-sample.cc
+ 4
− 7
View file @ f0de5bb3
...@@ -187,15 +187,12 @@ int main() { ...@@ -187,15 +187,12 @@ int main() {
leafView, p1Basis, frictionalNodes, nodalIntegrals); leafView, p1Basis, frictionalNodes, nodalIntegrals);
for (size_t run = 1; run <= runs; ++run) { for (size_t run = 1; run <= runs; ++run) {
f = 0; // f = neumann
VectorType neumannTerm(grid.size(grid.maxLevel(), dim));
assemble_neumann<GridType, GridType::LeafGridView, SmallVector, P1Basis>( assemble_neumann<GridType, GridType::LeafGridView, SmallVector, P1Basis>(
leafView, p1Basis, neumannNodes, neumannTerm); leafView, p1Basis, neumannNodes, f);
f += neumannTerm;
VectorType updateTerm(grid.size(grid.maxLevel(), dim)); // f += linear update
stiffnessMatrix.mv(u1, updateTerm); stiffnessMatrix.umv(u1, f);
f += updateTerm;
// {{{ Assemble terms for the nonlinearity // {{{ Assemble terms for the nonlinearity
// TODO: Random value // TODO: Random value
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment