diff --git a/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh b/dumux/material/fluidmatrixinteractions/2p/brookscorey.hh
index 2324200b5d62d791e425a0bcea57bfee7da45015..c43837329e9d5c40d1e3c217ab24fe227d69d89a 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;