Add verbose output for adaptive step size control

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

1 files changed, 24 insertions, 0 deletions

diff --git a/src/SetTime.c b/src/SetTime.c
index 32e893c..969c95a 100644
--- a/src/SetTime.c
+++ b/src/SetTime.c
@@ -123,6 +123,14 @@ void MoL_SetTime(CCTK_ARGUMENTS)
   
   CCTK_REAL beta;
 
+  if (adaptive_stepsize && ! CCTK_EQUALS(verbose, "none"))
+  {
+    CCTK_VInfo (CCTK_THORNSTRING, "Integrating from %g to %g with step size %g",
+                (double)(cctkGH->cctk_time - cctkGH->cctk_delta_time),
+                (double)cctkGH->cctk_time,
+                (double)cctkGH->cctk_delta_time);
+  }
+
   *Original_Time = cctkGH->cctk_time;
   *Original_Delta_Time = cctkGH->cctk_delta_time;
   cctkGH->cctk_time -= cctkGH->cctk_delta_time / cctkGH->cctk_timefac;
@@ -298,6 +306,10 @@ void MoL_ResetTime(CCTK_ARGUMENTS)
          + ((alpha_array87[substep] - 1)
             * (* Original_Delta_Time) / cctkGH->cctk_timefac));
   }
+  else if (CCTK_EQUALS(ODE_Method,"AB"))
+  {
+    CCTK_WARN (0, "internal error");
+  }
   else if (CCTK_EQUALS(ODE_Method,"RK2-MR-2:1"))
   {
     const int substep = MoL_Intermediate_Steps - (* MoL_Intermediate_Step);
@@ -363,6 +375,12 @@ void MoL_ResetTime(CCTK_ARGUMENTS)
 
   free(previous_times);
   previous_times = NULL;
+
+  if (adaptive_stepsize && ! CCTK_EQUALS(verbose, "none"))
+  {
+    CCTK_VInfo (CCTK_THORNSTRING, "Evaluating RHS at %g",
+                (double)cctkGH->cctk_time);
+  }
 }
 
  /*@@
@@ -471,6 +489,12 @@ void MoL_ResetDeltaTime(CCTK_ARGUMENTS)
          cctkGH->cctk_delta_time/cctkGH->cctk_timefac);
   fflush(stdout);
 #endif
+
+  if (adaptive_stepsize && ! CCTK_EQUALS(verbose, "none"))
+  {
+    CCTK_VInfo (CCTK_THORNSTRING, "Evaluating RHS with a time step of %g",
+                (double)cctkGH->cctk_delta_time);
+  }
 }
 
 /********************************************************************