diff --git a/src/ODE/ODESolvers/seulex/seulex.C b/src/ODE/ODESolvers/seulex/seulex.C index b07a9f498..8e6f6ac46 100644 --- a/src/ODE/ODESolvers/seulex/seulex.C +++ b/src/ODE/ODESolvers/seulex/seulex.C @@ -306,7 +306,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)