diff --git a/dune/tectonic/mydirectionalconvexfunction.hh b/dune/tectonic/mydirectionalconvexfunction.hh
index e815fa48e817d0b6bfecbaccb43b6a05765f1a2f..38c73d4646437cce0b2a6a76861c4aaa1dbd50ec 100644
--- a/dune/tectonic/mydirectionalconvexfunction.hh
+++ b/dune/tectonic/mydirectionalconvexfunction.hh
@@ -45,8 +45,9 @@ template <class Nonlinearity> class MyDirectionalConvexFunction {
     Vector uxv = u;
     Arithmetic::addProduct(uxv, x, v);
     phi.directionalSubDiff(uxv, v, D);
-    D[0] += A * x - b;
-    D[1] += A * x - b;
+    auto const Axmb = A * x - b;
+    D[0] += Axmb;
+    D[1] += Axmb;
   }
 
   void domain(Dune::Solvers::Interval<double> &domain) const {