diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index 2ff07f0e025de697a15a69bf7c67a375985c37c7..a74eda924b7c298533658ac30070bb30885f4d41 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -130,8 +130,10 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
 
     Interval<double> D;
     psi.subDiff(0, D);
-    if (D[1] >= 0)
-      assert(false);
+    if (D[1] > 0) {
+      assert(false); // NOTE: Numerical instability can actually get us here
+      return 0;
+    }
 
     int bisectionsteps = 0;
     Bisection bisection(0.0, 1.0, 1e-12, true, 0);                      // TODO