Peter Diener's RK65 and RK87 adaptive timestep integrators.

As yet not altered to do grid arrays.

As with RK45, adaptive timestepping does not work with mesh refinement.


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

1 files changed, 18 insertions, 0 deletions

diff --git a/schedule.ccl b/schedule.ccl
index add994d..3c5b82a 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -567,6 +567,24 @@ else if (CCTK_Equals(ODE_Method,"RK45"))
     LANG: C
   } "Updates calculated with the Runge-Kutta 45 method"
 }
+else if (CCTK_Equals(ODE_Method,"RK65")) 
+{
+  STORAGE: ErrorEstimate ErrorScalars
+
+  schedule MoL_RK65Add AS MoL_Add IN MoL_Step AFTER MoL_CalcRHS BEFORE MoL_PostStep
+  {
+    LANG: C
+  } "Updates calculated with the Runge-Kutta 65 method"
+}
+else if (CCTK_Equals(ODE_Method,"RK87")) 
+{
+  STORAGE: ErrorEstimate ErrorScalars
+
+  schedule MoL_RK87Add AS MoL_Add IN MoL_Step AFTER MoL_CalcRHS BEFORE MoL_PostStep
+  {
+    LANG: C
+  } "Updates calculated with the Runge-Kutta 87 method"
+}
 else if (CCTK_Equals(ODE_Method,"ICN"))
 {
   schedule MoL_ICNAdd AS MoL_Add IN MoL_Step AFTER MoL_CalcRHS BEFORE MoL_PostStep