diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index 97091a4537ab5717ee631fd9faacc3d5b5926e1a..7a28b0a567959a09d8160f2e362413bffc47b032 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -7,7 +7,7 @@
 #include <dune/common/nullptr.hh>
 #include <dune/common/parametertree.hh>
 
-#include <dune/solvers/common/staticmatrixtools.hh>
+#include <dune/fufem/arithmetic.hh>
 #include <dune/tnnmg/problem-classes/bisection.hh>
 
 #include "globalnonlinearity.hh"
@@ -129,7 +129,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
       typename MatrixType::row_type::ConstIterator it = problem.A[i].begin();
       typename MatrixType::row_type::ConstIterator end = problem.A[i].end();
       for (; it != end; ++it)
-        StaticMatrix::axpy(linearization.A[i][it.index()], 1.0, *it);
+        Arithmetic::addProduct(linearization.A[i][it.index()], 1.0, *it);
     }
 
     // compute nonlinearity part of hessian