From a4ed7cfd73db6ef8418ffb47f4237cc529f8c710 Mon Sep 17 00:00:00 2001 From: eschnett Date: Tue, 22 Jan 2013 21:00:28 +0000 Subject: MoL Update New integrator Euler. This is an explicit, first-order method, mostly useful for debugging. New integrators AB (Adams-Bashforth) with various orders. These are explicit integrators using several past timelevels to provide higher-order integration with a single RHS evaluation each. Introduce LinearCombination, a generic routine to calculate linear combinations of grid functions. This simplifies existing code, and can be overloaded if MoL should run on a device (e.g. with OpenCL). Replace cctk_lsh with cctk_ash where necessary. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@190 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b --- param.ccl | 28 +++++++++++++++++++++++++--- 1 file changed, 25 insertions(+), 3 deletions(-) (limited to 'param.ccl') diff --git a/param.ccl b/param.ccl index 43a1fd4..cce525f 100644 --- a/param.ccl +++ b/param.ccl @@ -91,6 +91,7 @@ KEYWORD ODE_Method "The ODE method use by MoL to do time integration" "Generic" :: "Generic Shu-Osher Runga-Kutta type" "ICN" :: "Iterative Crank Nicholson" "ICN-avg" :: "Iterative Crank Nicholson with averaging" + "Euler" :: "Euler" "RK2" :: "Efficient RK2" "RK3" :: "Efficient RK3" "RK4" :: "Efficient RK4" @@ -98,6 +99,7 @@ KEYWORD ODE_Method "The ODE method use by MoL to do time integration" "RK45CK" :: "RK45CK (Cash-Karp) with error estimation" "RK65" :: "RK65 with error estimation" "RK87" :: "RK87 with error estimation" + "AB" :: "Adams-Bashforth" "RK2-MR-2:1" :: "2nd order 2:1 multirate RK scheme based on RK2 due to Schlegel et al 2009" "RK4-MR-2:1" :: "3rd order 2:1 multirate RK scheme based on RK43 due to Schlegel et al 2009" "RK4-RK2" :: "RK4 as fast method and RK2 as slow method" @@ -120,6 +122,19 @@ BOOLEAN ICN_avg_swapped "Use swapped averages in ICN method?" { } "no" +KEYWORD AB_Type "If using the the AB method, which sort" +{ + "1" :: "same as forward Euler" + "2" :: "second order" + "3" :: "third order" + "4" :: "fourth order" + "5" :: "fifth order" +} "1" + +BOOLEAN AB_initially_reduce_order "Reduce order of accuracy initially so that no past timelevels of initial data are required" +{ +} "yes" + CCTK_INT MoL_Intermediate_Steps "Number of intermediate steps taken by the ODE method" { 1:* :: "Anything greater than 1" @@ -152,6 +167,11 @@ BOOLEAN copy_ID_after_MoL_PostStep \ { } "no" +BOOLEAN set_ID_boundaries "Should boundaries be overwritten (via synchronization, prolongation, boundary conditions) by MoL?" +{ +} "yes" + + # The default for this parameter corresponds to generic RK2 STRING Generic_Method_Descriptor "A string used to create a table containing the description of the generic method" @@ -221,6 +241,8 @@ KEYWORD verbose "How verbose should MoL be?" "extreme" :: "Everything you never wanted to know" } "normal" -BOOLEAN set_ID_boundaries "Should boundaries be overwritten (via synchronization, prolongation, boundary conditions) by MoL?" -{ -} "yes" + + +shares: Cactus + +USES CCTK_REAL cctk_initial_time -- cgit v1.2.3