diff options
| author | hawke <hawke@578cdeb0-5ea1-4b81-8215-5a3b8777ee0b> | 2005-01-27 15:10:17 +0000 |
|---|---|---|
| committer | hawke <hawke@578cdeb0-5ea1-4b81-8215-5a3b8777ee0b> | 2005-01-27 15:10:17 +0000 |
| commit | 495a7da9357fa3a487d626de33ec75dccb843100 (patch) | |
| tree | df3499e0f57dc60d5ad4fcd9f816ba69ff04bc43 /interface.ccl | |
| parent | c1088f5c4132c6a713f694a6640f1d2736794f71 (diff) | |
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
Diffstat (limited to 'interface.ccl')
| -rw-r--r-- | interface.ccl | 10 |
1 files changed, 9 insertions, 1 deletions
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" |