diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index dc52328bebba7f4c7686ae8e6c8667e8c11a86a3..a6aceca491e5dfad06a3c5f043a03c8bf05ca0bd 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -101,10 +101,6 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
 
     v /= vnorm; // Rescale for numerical stability
 
-    VectorType tmp = problem.f;
-    Arithmetic::addProduct(tmp, -1.0, problem.A, u);
-    double const localb = tmp * v;
-
     MyDirectionalConvexFunction<
         Dune::GlobalNonlinearity<MatrixType, VectorType>> const
     psi(computeDirectionalA(problem.A, v),