From 1cac1aaecc68b14a54df1c6f4b5b448b9fabfac3 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Fri, 17 Sep 2010 15:10:28 +0100
Subject: [PATCH] localBlended: Added support for "corrected" schemes Resolves
 bug report: http://www.openfoam.com/mantisbt/view.php?id=30

---
 src/ODE/ODESolvers/ODESolver/ODESolver.C | 16 ++++++++++++++--
 1 file changed, 14 insertions(+), 2 deletions(-)

diff --git a/src/ODE/ODESolvers/ODESolver/ODESolver.C b/src/ODE/ODESolvers/ODESolver/ODESolver.C
index b4341894b5..f0c04518be 100644
--- a/src/ODE/ODESolvers/ODESolver/ODESolver.C
+++ b/src/ODE/ODESolvers/ODESolver/ODESolver.C
@@ -60,6 +60,8 @@ void Foam::ODESolver::solve
 
     scalar x = xStart;
     scalar h = hEst;
+    scalar hNext = 0;
+    scalar hPrev = 0;
 
     for (label nStep=0; nStep<MAXSTP; nStep++)
     {
@@ -73,14 +75,24 @@ void Foam::ODESolver::solve
         if ((x + h - xEnd)*(x + h - xStart) > 0.0)
         {
             h = xEnd - x;
+            hPrev = hNext;
         }
 
-        scalar hNext, hDid;
+        hNext = 0;
+        scalar hDid;
         solve(ode, x, y, dydx_, eps, yScale_, h, hDid, hNext);
 
         if ((x - xEnd)*(xEnd - xStart) >= 0.0)
         {
-            hEst = hNext;
+            if (hPrev != 0)
+            {
+                hEst = hPrev;
+            }
+            else
+            {
+                hEst = hNext;
+            }
+
             return;
         }