diff options
-rw-r--r-- | src/SetTime.c | 24 |
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); + } } /******************************************************************** |