From b8b048e80fbf40d4f89fd832b6479c2819f3c2e4 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Tue, 28 Feb 2012 09:16:21 +0100
Subject: [PATCH] Do not assume that we find a direction of descent
---
dune/tectonic/myblockproblem.hh | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/dune/tectonic/myblockproblem.hh b/dune/tectonic/myblockproblem.hh
index 96fa97ac..352fdb0c 100644
--- a/dune/tectonic/myblockproblem.hh
+++ b/dune/tectonic/myblockproblem.hh
@@ -131,7 +131,8 @@ template <class MyConvexProblemTypeTEMPLATE> class MyBlockProblem {
Interval<double> D;
psi.subDiff(0, D);
if (D[1] > 0) {
- assert(false); // NOTE: Numerical instability can actually get us here
+ // NOTE: Numerical instability can actually get us here
+ assert(abs(D[1]) < 1e-15);
return 0;
}
--
GitLab