zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ODESolvers::seulex: Corrected 'expo' and added 'table_' resize

Browse Source
This commit is contained in:
Henry Weller
2016-07-12 20:02:16 +01:00
parent ff49c1b882
commit c535f6a3c4
1 changed files with 2 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@ -215,6 +215,7 @@ bool Foam::seulex::resize()
{ {
if (ODESolver::resize()) if (ODESolver::resize())
{ {
table_.shallowResize(kMaxx_, n_);
resizeField(dfdx_); resizeField(dfdx_);
resizeMatrix(dfdy_); resizeMatrix(dfdy_);
resizeMatrix(a_); resizeMatrix(a_);
@ -330,7 +331,7 @@ void Foam::seulex::solve
break; break;
} }
errOld = min(4*err, 1); errOld = min(4*err, 1);
scalar expo = 1/(k + 1); scalar expo = 1.0/(k + 1);
scalar facmin = pow(stepFactor3_, expo); scalar facmin = pow(stepFactor3_, expo);
scalar fac; scalar fac;
if (err == 0) if (err == 0)