Adaptive step size control using RK45. Due to Erik Schnetter.

Note that if you want to use this with Carpet you currently have to
use the development (darcs) version of Carpet together with the
parameter

carpet::adaptive_stepsize = "yes".


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@82 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b

1 files changed, 12 insertions, 1 deletions

diff --git a/src/ParamCheck.c b/src/ParamCheck.c
index bd91e5e..3bd1c21 100644
--- a/src/ParamCheck.c
+++ b/src/ParamCheck.c
@@ -128,7 +128,6 @@ int MoL_ParamCheck(CCTK_ARGUMENTS)
                    "number of intermediate steps must be 3");
   }
 
-
   if ( (CCTK_Equals(ODE_Method, "RK45")) &&
        ( !((MoL_Intermediate_Steps == 6)||(MoL_Num_Scratch_Levels > 5)) ) )
   {
@@ -137,6 +136,18 @@ int MoL_ParamCheck(CCTK_ARGUMENTS)
                    "and the number of scratch levels at least 6.");
   }
 
+  if (adaptive_stepsize)
+  {
+    if (CCTK_Equals(ODE_Method, "RK45"))
+    {
+      /* everything is fine, do nothing */
+    }
+    else
+    {
+      CCTK_PARAMWARN("Adaptive time step sizes are only possible with the RK45 solver");
+    }
+  }
+
   return 0;
   
 }