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