From 946105499f02d308f9aba1accfa8880e42de9e5d Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 14 Dec 2011 16:48:52 +0100
Subject: [PATCH] Use usmv
---
dune/tectonic/samplefunctional.hh | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/dune/tectonic/samplefunctional.hh b/dune/tectonic/samplefunctional.hh
index c6e0bc4e..1192ad6a 100644
--- a/dune/tectonic/samplefunctional.hh
+++ b/dune/tectonic/samplefunctional.hh
@@ -26,9 +26,8 @@ template <int dim> class SampleFunctional {
: A(A), b(b), phi(phi), ignore(ignore) {}
double operator()(SmallVector const &v) const {
- SmallVector y;
- A.mv(v, y); // Av
- y /= 2; // 1/2 Av
+ SmallVector y(0);
+ A.usmv(0.5, v, y); // 1/2 Av
y -= b; // 1/2 Av - b
return y * v + (*phi)(v); // <1/2 Av - b,v> + H(|v|)
}
--
GitLab