diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index afd98c30f4f167ddde4bfa0349f620d192c61190..19ea4959ef1613f028e179dc9fc51ebfe0691026 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -122,14 +122,14 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
       localA = <Av,v>
       localb = <b - Au, v>
     */
+
     MyDirectionalConvexFunction<Dune::GlobalNonlinearity<block_size>> const psi(
         localA, localb, problem.phi, u, v);
 
     Interval<double> D;
     psi.subDiff(0, D);
-    // FIXME: this should never happen to begin with
     if (D[1] >= 0)
-      return 0.0;
+      assert(false);
 
     int bisectionsteps = 0;
     Bisection bisection(0.0, 1.0, 1e-12, true, 0);                      // TODO