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, 8 insertions, 0 deletions

diff --git a/src/ParamCheck.c b/src/ParamCheck.c
index 47c573a..5316c9f 100644
--- a/src/ParamCheck.c
+++ b/src/ParamCheck.c
@@ -75,6 +75,13 @@ int MoL_ParamCheck(CCTK_ARGUMENTS)
                      "of scratch levels must be at least the "
                      "number of intermediate steps - 1");
     }
+    if ( (CCTK_Equals(Generic_Type, "Classic RK3"))&&
+         ((!(MoL_Intermediate_Steps == 3))||(!(MoL_Num_Scratch_Levels > 1))) )
+    {
+      CCTK_PARAMWARN("When using the classic RK3 evolver the "
+                     "number of intermediate steps must be 3 "
+                     "and the number of scratch levels at least 2");
+    }
     if (CCTK_Equals(Generic_Type, "Table"))
     {
       options_table =
@@ -121,6 +128,7 @@ int MoL_ParamCheck(CCTK_ARGUMENTS)
                    "number of intermediate steps must be 3");
   }
 
+
   return 0;
   
 }