From d76f1c26585fb23abdedbda233d98be51c6854a8 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Sun, 14 Jul 2013 15:33:27 +0200
Subject: [PATCH] [Cleanup] Use Axy()
---
dune/tectonic/mydirectionalconvexfunction.hh | 5 ++---
1 file changed, 2 insertions(+), 3 deletions(-)
diff --git a/dune/tectonic/mydirectionalconvexfunction.hh b/dune/tectonic/mydirectionalconvexfunction.hh
index 389eb410..ec870c3a 100644
--- a/dune/tectonic/mydirectionalconvexfunction.hh
+++ b/dune/tectonic/mydirectionalconvexfunction.hh
@@ -4,6 +4,7 @@
#ifndef MY_DIRECTIONAL_CONVEX_FUNCTION_HH
#define MY_DIRECTIONAL_CONVEX_FUNCTION_HH
+#include <dune/fufem/arithmetic.hh>
#include <dune/fufem/interval.hh>
/*
@@ -16,9 +17,7 @@
template <class MatrixType, class VectorType>
double computeDirectionalA(MatrixType const &A, VectorType const &v) {
- VectorType tmp(v.size());
- A.mv(v, tmp);
- return tmp * v;
+ return Arithmetic::Axy(A, v, v);
}
template <class MatrixType, class VectorType>
--
GitLab