zhyang-dev/OpenFOAM-6
1
0
Fork 0
mirror of https://github.com/OpenFOAM/OpenFOAM-6.git synced 2025-12-08 06:57:46 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

seulex: Reverting change to stability/monotonicity check

Browse Source 
This change cause some cases to fail, in particular
tutorials/lagrangian/sprayFoam/aachenBomb

Reopened https://bugs.openfoam.org/view.php?id=2972
This commit is contained in:
Henry Weller
2018-06-12 10:11:50 +01:00
parent 786318cabd
commit ee6a53adee
1 changed files with 3 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
src/ODE/ODESolvers/seulex/seulex.C
View File

@ -153,7 +153,8 @@ bool Foam::seulex::seul
LUBacksubstitute(a_, pivotIndices_, dy_); LUBacksubstitute(a_, pivotIndices_, dy_);
const scalar denom = max(1, dy1); // const scalar denom = max(1, dy1);
const scalar denom = min(1, dy1 + small);
scalar dy2 = 0; scalar dy2 = 0;
for (label i=0; i<n_; i++) for (label i=0; i<n_; i++)
{ {
@ -286,7 +287,7 @@ void Foam::seulex::solve
if (mag(dx) <= mag(x)*sqr(small)) if (mag(dx) <= mag(x)*sqr(small))
{ {
WarningInFunction WarningInFunction
<< "step size underflow :" << dx << endl; << "step size underflow :" << dx << endl;
} }
scalar errOld = 0; scalar errOld = 0;