From 873d2c9c23a617432f96ca5cb9ffe6c02661e9c1 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Wed, 6 Jun 2018 15:13:08 +0100
Subject: [PATCH] seulex: changed monotonicity check to correspond to original
 implementation

Resolves bug-report https://bugs.openfoam.org/view.php?id=2972
---
 src/ODE/ODESolvers/seulex/seulex.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/ODE/ODESolvers/seulex/seulex.C b/src/ODE/ODESolvers/seulex/seulex.C
index 90cd3077a..4df64c6c1 100644
--- a/src/ODE/ODESolvers/seulex/seulex.C
+++ b/src/ODE/ODESolvers/seulex/seulex.C
@@ -153,7 +153,7 @@ bool Foam::seulex::seul
 
             LUBacksubstitute(a_, pivotIndices_, dy_);
 
-            const scalar denom = min(1, dy1 + small);
+            const scalar denom = max(1, dy1);
             scalar dy2 = 0;
             for (label i=0; i<n_; i++)
             {