Add the Classic RK3 method (as a generic method, so use Generic_Type).

Agrees with other RK3's to floating point round off (except at
boundaries) for linear case. Uses more storage and is slower than
standard RK3 so I don't recommend it.

This showed up (so I fixed) a bug with the generic methods when used
with Carpet.


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

1 files changed, 4 insertions, 0 deletions

diff --git a/src/IndexArrays.c b/src/IndexArrays.c
index e0df297..096a745 100644
--- a/src/IndexArrays.c
+++ b/src/IndexArrays.c
@@ -282,6 +282,10 @@ void MoL_SetupIndexArrays(CCTK_ARGUMENTS)
       sprintf(infoline, "Generic method, options:\n %s\n",
               Generic_Method_Descriptor);
     }
+    else if (CCTK_EQUALS(Generic_Type,"Classic RK3")) 
+    {
+      sprintf(infoline, "Classic Runge-Kutta 3");
+    }
     else
     {
       CCTK_WARN(0, "Generic_Type not recognized!");