From 2075f04a865be8ddead6587cc21170894a7ea964 Mon Sep 17 00:00:00 2001
From: Timo Koch <timo.koch@iws.uni-stuttgart.de>
Date: Mon, 1 Feb 2021 13:26:44 +0100
Subject: [PATCH] [brookscorey] Add epsilon to regularization query
---
dumux/material/fluidmatrixinteractions/2p/brookscorey.hh | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
index 9b7a19caea..2324200b5d 100644
--- a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
+++ b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
@@ -634,7 +634,7 @@ public:
if (swe <= pcLowSwe_)
return pcLowSwePcValue_ + pcDerivativeLowSw_*(swe - pcLowSwe_);
- else if (swe > 1.0)
+ else if (swe > 1.0 - std::numeric_limits<Scalar>::epsilon())
return pcDerivativeHighSwEnd_*(swe - 1.0) + entryPressure_;
else
@@ -649,7 +649,7 @@ public:
if (swe < pcLowSwe_)
return pcDerivativeLowSw_;
- else if (swe > 1.0)
+ else if (swe > 1.0 - std::numeric_limits<Scalar>::epsilon())
return pcDerivativeHighSwEnd_;
else
--
GitLab