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

1 files changed, 25 insertions, 3 deletions

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