Erik Schnetter's implementation of ICN with averaging, so the intermediate steps are always at t+dt.

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

1 files changed, 20 insertions, 1 deletions

diff --git a/schedule.ccl b/schedule.ccl
index 9de0d8f..4321dcf 100644
--- a/schedule.ccl
+++ b/schedule.ccl
@@ -220,6 +220,18 @@ schedule GROUP MoL_Step WHILE MoL::MoL_Intermediate_Step IN MoL_Evolution AFTER
   LANG: C
 } "The loop over the intermediate steps for the ODE integrator"
 
+####################################################
+###  The time integrator prepares the time step  ###
+####################################################
+
+if (CCTK_Equals(ODE_Method,"ICN-avg")) 
+{
+  schedule MoL_ICNAverage AS MoL_Prepare IN MoL_Step BEFORE MoL_CalcRHS
+  {
+    LANG: C
+  } "Averages the time levels for the averaging ICN method"
+}
+
 #####################################################
 ###  The group where all the physics takes place  ###
 #####################################################
@@ -247,13 +259,20 @@ else if (CCTK_Equals(ODE_Method,"RK2"))
     LANG: C
   } "Updates calculated with the efficient Runge-Kutta 2 method"
 }
-else if (CCTK_Equals(ODE_Method,"ICN")) 
+else if (CCTK_Equals(ODE_Method,"ICN"))
 {
   schedule MoL_ICNAdd AS MoL_Add IN MoL_Step AFTER MoL_CalcRHS BEFORE MoL_PostStep
   {
     LANG: C
   } "Updates calculated with the efficient ICN method"
 }
+else if (CCTK_Equals(ODE_Method,"ICN-avg"))
+{
+  schedule MoL_ICNAdd AS MoL_Add IN MoL_Step AFTER MoL_CalcRHS BEFORE MoL_PostStep
+  {
+    LANG: C
+  } "Updates calculated with the averaging ICN method"
+}
 
 ##################################################
 ###  Physics thorns can apply boundaries and   ###