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

diff --git a/src/Registration.c b/src/Registration.c
index 2dd4767..2d851c3 100644
--- a/src/Registration.c
+++ b/src/Registration.c
@@ -735,8 +735,8 @@ CCTK_INT MoL_RegisterEvolvedReal(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
                CCTK_VarName(RHSIndex)); 
   }
 
-  numtimelevs1 = CCTK_NumTimeLevelsFromVarI(EvolvedIndex);
-  numtimelevs2 = CCTK_NumTimeLevelsFromVarI(RHSIndex);
+  numtimelevs1 = CCTK_MaxTimeLevelsVI(EvolvedIndex);
+  numtimelevs2 = CCTK_MaxTimeLevelsVI(RHSIndex);
   
   if ( (numtimelevs1 < 0) || (numtimelevs2 < 0) ) 
   {
@@ -864,7 +864,7 @@ CCTK_INT MoL_RegisterConstrainedReal(CCTK_INT ConstrainedIndex)
                CCTK_VarName(ConstrainedIndex)); 
   }
   
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(ConstrainedIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(ConstrainedIndex);
 
   if (numtimelevs < 1) {
 
@@ -962,7 +962,7 @@ CCTK_INT MoL_RegisterSaveAndRestoreReal(CCTK_INT SandRIndex)
                CCTK_VarName(SandRIndex)); 
   }
   
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(SandRIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(SandRIndex);
 
   if (numtimelevs < 1) {
 
@@ -1164,8 +1164,8 @@ CCTK_INT MoL_RegisterEvolvedComplex(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
                CCTK_VarName(RHSIndex)); 
   }
 
-  numtimelevs1 = CCTK_NumTimeLevelsFromVarI(EvolvedIndex);
-  numtimelevs2 = CCTK_NumTimeLevelsFromVarI(RHSIndex);
+  numtimelevs1 = CCTK_MaxTimeLevelsVI(EvolvedIndex);
+  numtimelevs2 = CCTK_MaxTimeLevelsVI(RHSIndex);
   
   if ( (numtimelevs1 < 0) || (numtimelevs2 < 0) ) 
   {
@@ -1279,7 +1279,7 @@ CCTK_INT MoL_RegisterConstrainedComplex(CCTK_INT ConstrainedIndex)
                CCTK_VarName(ConstrainedIndex)); 
   }
 
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(ConstrainedIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(ConstrainedIndex);
 
   if (numtimelevs < 1) {
 
@@ -1373,7 +1373,7 @@ CCTK_INT MoL_RegisterSaveAndRestoreComplex(CCTK_INT SandRIndex)
                CCTK_VarName(SandRIndex)); 
   }
 
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(SandRIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(SandRIndex);
 
   if (numtimelevs < 1) {
 
@@ -1574,8 +1574,8 @@ CCTK_INT MoL_RegisterEvolvedArray(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
                CCTK_VarName(RHSIndex)); 
   }
 
-  numtimelevs1 = CCTK_NumTimeLevelsFromVarI(EvolvedIndex);
-  numtimelevs2 = CCTK_NumTimeLevelsFromVarI(RHSIndex);
+  numtimelevs1 = CCTK_MaxTimeLevelsVI(EvolvedIndex);
+  numtimelevs2 = CCTK_MaxTimeLevelsVI(RHSIndex);
   
   if ( (numtimelevs1 < 0) || (numtimelevs2 < 0) ) 
   {
@@ -1683,7 +1683,7 @@ CCTK_INT MoL_RegisterConstrainedArray(CCTK_INT ConstrainedIndex)
                CCTK_VarName(ConstrainedIndex)); 
   }
 
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(ConstrainedIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(ConstrainedIndex);
 
   if (numtimelevs < 1) {
 
@@ -1771,7 +1771,7 @@ CCTK_INT MoL_RegisterSaveAndRestoreArray(CCTK_INT SandRIndex)
                CCTK_VarName(SandRIndex)); 
   }
 
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(SandRIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(SandRIndex);
 
   if (numtimelevs < 1) {
 
@@ -1975,8 +1975,8 @@ CCTK_INT MoL_RegisterEvolvedComplexArray(CCTK_INT EvolvedIndex,
                CCTK_VarName(RHSIndex)); 
   }
 
-  numtimelevs1 = CCTK_NumTimeLevelsFromVarI(EvolvedIndex);
-  numtimelevs2 = CCTK_NumTimeLevelsFromVarI(RHSIndex);
+  numtimelevs1 = CCTK_MaxTimeLevelsVI(EvolvedIndex);
+  numtimelevs2 = CCTK_MaxTimeLevelsVI(RHSIndex);
   
   if ( (numtimelevs1 < 0) || (numtimelevs2 < 0) ) 
   {
@@ -2090,7 +2090,7 @@ CCTK_INT MoL_RegisterConstrainedComplexArray(CCTK_INT ConstrainedIndex)
                CCTK_VarName(ConstrainedIndex)); 
   }
   
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(ConstrainedIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(ConstrainedIndex);
 
   if (numtimelevs < 1) {
 
@@ -2184,7 +2184,7 @@ CCTK_INT MoL_RegisterSaveAndRestoreComplexArray(CCTK_INT SandRIndex)
                CCTK_VarName(SandRIndex)); 
   }
   
-  numtimelevs = CCTK_NumTimeLevelsFromVarI(SandRIndex);
+  numtimelevs = CCTK_MaxTimeLevelsVI(SandRIndex);
 
   if (numtimelevs < 1) {