From eae928738528296c7a0bb6d6b139cda94a9695f2 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Fri, 1 Jun 2012 11:30:17 +0200
Subject: [PATCH] b = 0; A.umv(u,b) -> A.mv(u,b)
---
dune/tectonic/myblockproblem.hh | 3 +--
1 file changed, 1 insertion(+), 2 deletions(-)
diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index 2c91cfbd..a69f2a0b 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -182,8 +182,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
// compute quadratic part of gradient
linearization.b.resize(u.size());
- linearization.b = 0.0;
- problem.A.template umv<VectorType, VectorType>(u, linearization.b);
+ problem.A.mv(u, linearization.b);
linearization.b -= problem.f;
// compute nonlinearity part of gradient
--
GitLab