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

b = 0; A.umv(u,b) -> A.mv(u,b)

parent 33837ef6
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
dune/tectonic/myblockproblem.hh
+ 1
− 2
View file @ eae92873
...@@ -182,8 +182,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem { ...@@ -182,8 +182,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
// compute quadratic part of gradient // compute quadratic part of gradient
linearization.b.resize(u.size()); linearization.b.resize(u.size());
linearization.b = 0.0; problem.A.mv(u, linearization.b);
problem.A.template umv<VectorType, VectorType>(u, linearization.b);
linearization.b -= problem.f; linearization.b -= problem.f;
// compute nonlinearity part of gradient // compute nonlinearity part of gradient
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment