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

localBlended: Added support for "corrected" schemes

Browse Source 
Resolves bug report: http://www.openfoam.com/mantisbt/view.php?id=30
This commit is contained in:
Henry
2010-09-17 15:10:28 +01:00
parent 4e7fa4989b
commit 1cac1aaecc
1 changed files with 14 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
src/ODE/ODESolvers/ODESolver/ODESolver.C
View File

@ -60,6 +60,8 @@ void Foam::ODESolver::solve
scalar x = xStart; scalar x = xStart;
scalar h = hEst; scalar h = hEst;
scalar hNext = 0;
scalar hPrev = 0;
for (label nStep=0; nStep<MAXSTP; nStep++) for (label nStep=0; nStep<MAXSTP; nStep++)
{ {
@ -73,14 +75,24 @@ void Foam::ODESolver::solve
if ((x + h - xEnd)*(x + h - xStart) > 0.0) if ((x + h - xEnd)*(x + h - xStart) > 0.0)
{ {
h = xEnd - x; h = xEnd - x;
hPrev = hNext;
} }
scalar hNext, hDid; hNext = 0;
scalar hDid;
solve(ode, x, y, dydx_, eps, yScale_, h, hDid, hNext); solve(ode, x, y, dydx_, eps, yScale_, h, hDid, hNext);
if ((x - xEnd)*(xEnd - xStart) >= 0.0) if ((x - xEnd)*(xEnd - xStart) >= 0.0)
{ {
hEst = hNext; if (hPrev != 0)
{
hEst = hPrev;
}
else
{
hEst = hNext;
}
return; return;
} }