From fb30868c8bf1234fa22900f2940751efde306679 Mon Sep 17 00:00:00 2001 From: hawke Date: Mon, 6 Sep 2004 17:12:28 +0000 Subject: Runge-Kutta (Fehlberg) 45 with error estimation. Fourth order accurate evolution with an additional fifth order step for error estimation. How much sense the error makes is unclear to me, but hey. For the moment the error is stored in an internal MoL array ErrorEstimate; there is one per evolved variable. At a later point this may be moved out to user thorns who can register their own etc. As the implementation uses 6 evaluations of the RHS (necessary) and 6 levels of scratch space (one more than necessary - laziness kicked in) then this is very expensive. This is a partial fix for PR/1840. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@74 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b --- interface.ccl | 4 ++++ 1 file changed, 4 insertions(+) (limited to 'interface.ccl') diff --git a/interface.ccl b/interface.ccl index ac744ed..f9f45f5 100644 --- a/interface.ccl +++ b/interface.ccl @@ -220,3 +220,7 @@ CCTK_REAL ArraySandRScratchSpace TYPE = ARRAY DIM = 1 SIZE = MoL_Max_Evolved_Arr #CCTK_COMPLEX ComplexArrayScratchSpace[MoL_Num_Scratch_Levels] TYPE = ARRAY DIM = 1 SIZE = MoL_Max_Evolved_ComplexArray_Size Timelevels = 1 #CCTK_COMPLEX ComplexArraySandRScratchSpace TYPE = ARRAY DIM = 1 SIZE = MoL_Max_Evolved_ComplexArray_Size Timelevels = 1 + +#Error vector for RK45 + +CCTK_REAL ErrorEstimate[MoL_Num_Evolved_Vars] TYPE=GF Timelevels = 1 tags='Prolongation="None"' -- cgit v1.2.3