From 495a7da9357fa3a487d626de33ec75dccb843100 Mon Sep 17 00:00:00 2001 From: hawke Date: Thu, 27 Jan 2005 15:10:17 +0000 Subject: Adaptive step size control using RK45. Due to Erik Schnetter. Note that if you want to use this with Carpet you currently have to use the development (darcs) version of Carpet together with the parameter carpet::adaptive_stepsize = "yes". git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@82 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b --- interface.ccl | 10 +++++++++- 1 file changed, 9 insertions(+), 1 deletion(-) (limited to 'interface.ccl') diff --git a/interface.ccl b/interface.ccl index f9f45f5..40e1c79 100644 --- a/interface.ccl +++ b/interface.ccl @@ -194,6 +194,7 @@ CCTK_REAL RKBetaCoefficients TYPE=ARRAY DIM=1 SIZE=MoL_Intermediate_Steps DISTRI CCTK_INT MoL_Counters TYPE = SCALAR { MoL_Intermediate_Step + MoL_Stepsize_Bad } "The counter for the time integration method" CCTK_REAL MoL_Original_Time TYPE = SCALAR @@ -221,6 +222,13 @@ CCTK_REAL ArraySandRScratchSpace TYPE = ARRAY DIM = 1 SIZE = MoL_Max_Evolved_Arr #CCTK_COMPLEX ComplexArraySandRScratchSpace TYPE = ARRAY DIM = 1 SIZE = MoL_Max_Evolved_ComplexArray_Size Timelevels = 1 -#Error vector for RK45 +#Error vector and scalars for RK45 CCTK_REAL ErrorEstimate[MoL_Num_Evolved_Vars] TYPE=GF Timelevels = 1 tags='Prolongation="None"' + +CCTK_REAL ErrorScalars TYPE=SCALAR +{ + Error + Count + EstimatedDt +} "Global error estimate" -- cgit v1.2.3