From 691a71e2b209f551e466df869219c04b19487ba7 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Dennis=20Gl=C3=A4ser?= <dennis.glaeser@iws.uni-stuttgart.de>
Date: Mon, 1 Feb 2021 19:19:52 +0100
Subject: [PATCH] [brookscorey] unique way of checking regularization region

---
 dumux/material/fluidmatrixinteractions/2p/brookscorey.hh | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
index 2324200b5d..c43837329e 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 - std::numeric_limits<Scalar>::epsilon())
+        else if (swe >= 1.0)
             return pcDerivativeHighSwEnd_*(swe - 1.0) + entryPressure_;
 
         else
@@ -646,10 +646,10 @@ public:
      */
     OptionalScalar<Scalar> dpc_dswe(const Scalar swe) const
     {
-        if (swe < pcLowSwe_)
+        if (swe <= pcLowSwe_)
             return pcDerivativeLowSw_;
 
-        else if (swe > 1.0 - std::numeric_limits<Scalar>::epsilon())
+        else if (swe >= 1.0)
             return pcDerivativeHighSwEnd_;
 
         else
@@ -704,7 +704,7 @@ public:
      */
     OptionalScalar<Scalar> dkrw_dswe(const Scalar swe) const
     {
-        if (swe < 0.0)
+        if (swe <= 0.0)
             return 0.0;
         else if (swe >= 1.0)
             return 0.0;
-- 
GitLab