From b1ce18096cbebfd539b01ff13cf16134059356fd Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Mon, 16 Feb 2015 22:16:50 +0000
Subject: [PATCH] seulex: Reduce the underflow time-step
---
src/ODE/ODESolvers/seulex/seulex.C | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/ODE/ODESolvers/seulex/seulex.C b/src/ODE/ODESolvers/seulex/seulex.C
index cc788f93082..9dfa3498f91 100644
--- a/src/ODE/ODESolvers/seulex/seulex.C
+++ b/src/ODE/ODESolvers/seulex/seulex.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2013-2015 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -250,7 +250,7 @@ void Foam::seulex::solve
firstk = false;
step.reject = false;
- if (mag(dx) <= mag(x)*SMALL)
+ if (mag(dx) <= mag(x)*sqr(SMALL))
{
WarningIn("seulex::solve(scalar& x, scalarField& y, stepState&")
<< "step size underflow :" << dx << endl;
--
GitLab