From 7f7d351b741bf6406366a043cac98de56d2d44dd Mon Sep 17 00:00:00 2001 From: Henry Weller Date: Wed, 6 Jun 2018 15:11:07 +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 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/ODE/ODESolvers/seulex/seulex.C b/src/ODE/ODESolvers/seulex/seulex.C index e28f4c0f5..c6402a414 100644 --- a/src/ODE/ODESolvers/seulex/seulex.C +++ b/src/ODE/ODESolvers/seulex/seulex.C @@ -2,7 +2,7 @@ ========= | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox \\ / O peration | - \\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation + \\ / A nd | Copyright (C) 2013-2018 OpenFOAM Foundation \\/ M anipulation | ------------------------------------------------------------------------------- License @@ -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