diff --git a/src/compute_state.cc b/src/compute_state.cc
index 599c0de4d75336d84f9835cb4b84f14d2d450841..4580c45903eec3c26e5f6842d4cca2417f141102 100644
--- a/src/compute_state.cc
+++ b/src/compute_state.cc
@@ -4,9 +4,10 @@
 #include <dune/common/fmatrix.hh>
 #include <dune/fufem/interval.hh>
 
-#include <dune/tectonic/mydirectionalconvexfunction.hh>
 #include <dune/tnnmg/problem-classes/bisection.hh>
 
+#include <dune/tectonic/mydirectionalconvexfunction.hh>
+
 #include "compute_state.hh"
 
 // The nonlinearity exp(-x)
@@ -34,8 +35,11 @@ class DecayingExponential {
 
 double compute_state_update_bisection(double h, double unorm, double L,
                                       double old_state) {
+  DecayingExponential::VectorType const start(0);
+  DecayingExponential::VectorType const direction(1);
+  DecayingExponential const phi(h);
   MyDirectionalConvexFunction<DecayingExponential> const J(
-      1.0, old_state - unorm / L, DecayingExponential(h), 0, 1);
+      1.0, old_state - unorm / L, phi, start, direction);
   int bisectionsteps = 0;
   Bisection const bisection(0.0, 1.0, 1e-12, true, 0);    // TODO
   return bisection.minimize(J, 0.0, 0.0, bisectionsteps); // TODO