Add two methods:

RK3. The optimized version of the TVD RK3 solver. Requires no scratch space so is about as efficient as ICN, but third order.

Generic method from a parameter table. By specifying the number of intermediate steps and the alpha and beta arrays, create your own method at parameter time. Not well (or at all) documented because it doesn't seem to work correctly at the moment.

Some tidying of extraneous code as well.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@29 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b

1 files changed, 8 insertions, 0 deletions

diff --git a/src/GenericRK.c b/src/GenericRK.c
index 60edf15..19ffd0b 100644
--- a/src/GenericRK.c
+++ b/src/GenericRK.c
@@ -150,6 +150,10 @@ void MoL_GenericRKAdd(CCTK_ARGUMENTS)
       if (scratchstep) 
       {
         /*
+          The following would work if all drivers considered
+          a vector group to have contiguous storage.
+        */
+        /*
           ScratchVar = &ScratchSpace[(var * MoL_Num_Scratch_Levels + 
                                    scratchindex) * totalsize];
         */
@@ -191,6 +195,10 @@ void MoL_GenericRKAdd(CCTK_ARGUMENTS)
       UpdateVar = (CCTK_REAL *)CCTK_VarDataPtrI(cctkGH, 0, 
                                                 EvolvedVariableIndex[var]);
       /*
+        The following would work if all drivers considered
+        a vector group to have contiguous storage.
+      */
+      /*
         ScratchVar = &ScratchSpace[(var * MoL_Num_Scratch_Levels + 
                                  MoL_Intermediate_Steps -
                                  (*MoL_Intermediate_Step)) * totalsize];