From 12c3025d205683733326cbfba384f69f665fcd62 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Tue, 12 Jul 2016 20:02:16 +0100
Subject: [PATCH] ODESolvers::seulex: Corrected 'expo' and added 'table_'
 resize

---
 src/ODE/ODESolvers/seulex/seulex.C | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

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