diff --git a/src/ODE/ODESolvers/seulex/seulex.C b/src/ODE/ODESolvers/seulex/seulex.C
index cc788f93082fa0dda41b16ca13ce417ecd2937a2..9dfa3498f91eee99c2a05e4a4f8f5068e9a91af0 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;