seulex: changed monotonicity check to correspond to original implementation

Resolves bug-report https://bugs.openfoam.org/view.php?id=2972
This commit is contained in:
Henry Weller
2018-06-06 15:11:07 +01:00
parent 68e8507efb
commit 7f7d351b74

View File

@ -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<n_; i++)
{