diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index d9b01fb0641d46627ac7107ec072d1816ab0cd52..4de7a709ec01e03900f7148582d877ce89a43552 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -149,7 +149,7 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
     for (size_t i = 0; i < problem.A.N(); ++i) {
       auto const end = problem.A[i].end();
       for (auto it = problem.A[i].begin(); it != end; ++it)
-        Arithmetic::addProduct(linearization.A[i][it.index()], 1.0, *it);
+        linearization.A[i][it.index()] += *it;
     }
 
     // compute nonlinearity part of hessian