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, 4 insertions, 4 deletions

diff --git a/src/ICN.c b/src/ICN.c
index fb71783..af074e3 100644
--- a/src/ICN.c
+++ b/src/ICN.c
@@ -105,7 +105,7 @@ void MoL_ICNAdd(CCTK_ARGUMENTS)
   totalsize = 1;
   for (arraydim = 0; arraydim < cctk_dim; arraydim++)
   {
-    totalsize *= cctk_lsh[arraydim];
+    totalsize *= cctk_ash[arraydim];
   }
 
 #ifdef MOLDEBUG
@@ -149,7 +149,7 @@ void MoL_ICNAdd(CCTK_ARGUMENTS)
     arraytotalsize = 1;
     for (arraydim = 0; arraydim < arraydata.dim; arraydim++)
     {
-      arraytotalsize *= arraydata.lsh[arraydim];
+      arraytotalsize *= arraydata.ash[arraydim];
     }
 
 /*     CCTK_VWarn(1, __LINE__, __FILE__, CCTK_THORNSTRING,  */
@@ -256,7 +256,7 @@ void MoL_ICNAverage(CCTK_ARGUMENTS)
   totalsize = 1;
   for (arraydim = 0; arraydim < cctk_dim; arraydim++)
   {
-    totalsize *= cctk_lsh[arraydim];
+    totalsize *= cctk_ash[arraydim];
   }
 
 #ifdef MOLDEBUG
@@ -298,7 +298,7 @@ void MoL_ICNAverage(CCTK_ARGUMENTS)
     arraytotalsize = 1;
     for (arraydim = 0; arraydim < arraydata.dim; arraydim++)
     {
-      arraytotalsize *= arraydata.lsh[arraydim];
+      arraytotalsize *= arraydata.ash[arraydim];
     }
 
 #pragma omp parallel for