diff --git a/dune/tectonic/localfriction.hh b/dune/tectonic/localfriction.hh
index bfac47f69358aa513b02643ce3144778641f355a..b9cded4b4cd785dc32b522531f8707cf3fa36061 100644
--- a/dune/tectonic/localfriction.hh
+++ b/dune/tectonic/localfriction.hh
@@ -8,6 +8,7 @@
 #include <dune/common/fmatrix.hh>
 #include <dune/common/shared_ptr.hh>
 
+#include <dune/fufem/arithmetic.hh>
 #include <dune/fufem/interval.hh>
 
 #include "frictionpotential.hh"
@@ -120,7 +121,7 @@ template <int dimension> class LocalFriction {
       return;
 
     if (xnorm >= smp)
-      y.axpy(func->differential(xnorm) / xnorm, x);
+      Arithmetic::addProduct(y, func->differential(xnorm) / xnorm, x);
   }
 
   void directionalDomain(VectorType const &, VectorType const &,