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, 12 insertions, 0 deletions

diff --git a/src/SetTime.c b/src/SetTime.c
index 927f5fa..b7634c9 100644
--- a/src/SetTime.c
+++ b/src/SetTime.c
@@ -80,6 +80,10 @@ int MoL_SetTime(CCTK_ARGUMENTS)
   {
     cctkGH->cctk_delta_time = 0.5*(*Original_Delta_Time);
   }
+  else if (CCTK_EQUALS(ODE_Method,"ICN-avg"))
+  {
+    cctkGH->cctk_delta_time = *Original_Delta_Time;
+  }
   else if (CCTK_EQUALS(ODE_Method,"Generic"))
   {
     beta = RKBetaCoefficients[0];
@@ -130,6 +134,10 @@ int MoL_ResetTime(CCTK_ARGUMENTS)
   {
     cctkGH->cctk_time = (*Original_Time)-0.5*(*Original_Delta_Time)/cctkGH->cctk_timefac;
   }
+  else if (CCTK_EQUALS(ODE_Method,"ICN-avg"))
+  {
+    cctkGH->cctk_time = (*Original_Time);
+  }
   else if (CCTK_EQUALS(ODE_Method,"Generic"))
   {
     previous_times[0] = (*Original_Time) -
@@ -223,6 +231,10 @@ int MoL_ResetDeltaTime(CCTK_ARGUMENTS)
       cctkGH->cctk_delta_time = 0.5*(*Original_Delta_Time);
     }
   }
+  else if (CCTK_EQUALS(ODE_Method,"ICN-avg"))
+  {
+    cctkGH->cctk_delta_time = (*Original_Delta_Time);
+  }
   else if (CCTK_EQUALS(ODE_Method,"Generic"))
   {
     cctkGH->cctk_delta_time = RKBetaCoefficients[MoL_Intermediate_Steps -