From a8fbd505eab8afc48539e4a40b869ab46d5c31a8 Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Thu, 11 Jun 2020 23:23:19 +0200
Subject: [PATCH] [fix][geometry] Fix base epsilon for 1d case bbox-bbox
---
dumux/common/geometry/boundingboxtree.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dumux/common/geometry/boundingboxtree.hh b/dumux/common/geometry/boundingboxtree.hh
index 2e74524da2..3541867287 100644
--- a/dumux/common/geometry/boundingboxtree.hh
+++ b/dumux/common/geometry/boundingboxtree.hh
@@ -406,7 +406,7 @@ template<int dimworld, class ctypea, class ctypeb, typename std::enable_if_t<dim
inline bool intersectsBoundingBoxBoundingBox(const ctypea* a, const ctypeb* b)
{
using ctype = typename Dune::PromotionTraits<ctypea, ctypeb>::PromotedType;
- static constexpr ctype eps_ = 1.0e-1;
+ static constexpr ctype eps_ = 1.0e-7;
const ctype eps0 = eps_*std::max(b[1]-b[0], a[1]-a[0]);
return b[0] - eps0 <= a[1] && a[0] <= b[1] + eps0;
}
--
GitLab