Add support for evolving complex GFs and (real and complex) GAs.

Only works with ICN or RK2 for now - in fact this commit may break the generic RK methods temporarily.

Note the documentation isn't quite right - there's no longer a seperate function for each different type...


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

1 files changed, 124 insertions, 1 deletions

diff --git a/src/ChangeType.c b/src/ChangeType.c
index 5e229b4..966445e 100644
--- a/src/ChangeType.c
+++ b/src/ChangeType.c
@@ -47,6 +47,34 @@ CCTK_INT MoL_ChangeToSaveAndRestore(CCTK_INT SandRIndex);
 
 CCTK_INT MoL_ChangeToNone(CCTK_INT RemoveIndex);
 
+
+CCTK_INT MoL_ChangeToEvolvedComplex(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex);
+
+CCTK_INT MoL_ChangeToConstrainedComplex(CCTK_INT ConstrainedIndex);
+
+CCTK_INT MoL_ChangeToSaveAndRestoreComplex(CCTK_INT SandRIndex);
+
+CCTK_INT MoL_ChangeToNoneComplex(CCTK_INT RemoveIndex);
+
+
+CCTK_INT MoL_ChangeToEvolvedArray(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex);
+
+CCTK_INT MoL_ChangeToConstrainedArray(CCTK_INT ConstrainedIndex);
+
+CCTK_INT MoL_ChangeToSaveAndRestoreArray(CCTK_INT SandRIndex);
+
+CCTK_INT MoL_ChangeToNoneArray(CCTK_INT RemoveIndex);
+
+
+CCTK_INT MoL_ChangeToEvolvedComplexArray(CCTK_INT EvolvedIndex, 
+                                         CCTK_INT RHSIndex);
+
+CCTK_INT MoL_ChangeToConstrainedComplexArray(CCTK_INT ConstrainedIndex);
+
+CCTK_INT MoL_ChangeToSaveAndRestoreComplexArray(CCTK_INT SandRIndex);
+
+CCTK_INT MoL_ChangeToNoneComplexArray(CCTK_INT RemoveIndex);
+
 /********************************************************************
  *********************     Local Data   *****************************
  ********************************************************************/
@@ -441,7 +469,7 @@ CCTK_INT MoL_ChangeToNone(CCTK_INT RemoveIndex)
 
   CCTK_INT index, usedindex;
   CCTK_INT vartype; /* See the defines at the top of file */
-  CCTK_INT timelevs;
+/*    CCTK_INT timelevs; */
   
   vartype = 0;
   usedindex = -1;
@@ -531,13 +559,107 @@ CCTK_INT MoL_ChangeToNone(CCTK_INT RemoveIndex)
   
 }
 
+/*
+  Currently no-op test functions.
+*/
+
+
+CCTK_INT MoL_ChangeToEvolvedComplex(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
+{
+  CCTK_INT dummy;
+  dummy = EvolvedIndex;
+  dummy = RHSIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToConstrainedComplex(CCTK_INT ConstrainedIndex)
+{
+  CCTK_INT dummy;
+  dummy = ConstrainedIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToSaveAndRestoreComplex(CCTK_INT SandRIndex)
+{
+  CCTK_INT dummy;
+  dummy = SandRIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToNoneComplex(CCTK_INT RemoveIndex)
+{
+  CCTK_INT dummy;
+  dummy = RemoveIndex;
+  return 0;
+}
+
+
+CCTK_INT MoL_ChangeToEvolvedArray(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
+{
+  CCTK_INT dummy;
+  dummy = EvolvedIndex;
+  dummy = RHSIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToConstrainedArray(CCTK_INT ConstrainedIndex)
+{
+  CCTK_INT dummy;
+  dummy = ConstrainedIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToSaveAndRestoreArray(CCTK_INT SandRIndex)
+{
+  CCTK_INT dummy;
+  dummy = SandRIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToNoneArray(CCTK_INT RemoveIndex)
+{
+  CCTK_INT dummy;
+  dummy = RemoveIndex;
+  return 0;
+}
 
 
+CCTK_INT MoL_ChangeToEvolvedComplexArray(CCTK_INT EvolvedIndex, 
+                                         CCTK_INT RHSIndex)
+{
+  CCTK_INT dummy;
+  dummy = EvolvedIndex;
+  dummy = RHSIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToConstrainedComplexArray(CCTK_INT ConstrainedIndex)
+{
+  CCTK_INT dummy;
+  dummy = ConstrainedIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToSaveAndRestoreComplexArray(CCTK_INT SandRIndex)
+{
+  CCTK_INT dummy;
+  dummy = SandRIndex;
+  return 0;
+}
+
+CCTK_INT MoL_ChangeToNoneComplexArray(CCTK_INT RemoveIndex)
+{
+  CCTK_INT dummy;
+  dummy = RemoveIndex;
+  return 0;
+}
+
 /* 
    Fortran wrappers for the above functions. 
    Should be replaced by using function aliasing eventually.
  */
 
+/*
 void CCTK_FCALL CCTK_FNAME(MoL_ChangeToEvolved)(int *ierr, 
                                             CCTK_INT *EvolvedIndex,
                                             CCTK_INT *RHSIndex)
@@ -566,6 +688,7 @@ void CCTK_FCALL CCTK_FNAME(MoL_ChangeToNone)(int *ierr,
   *ierr = MoL_ChangeToNone(*EvolvedIndex);
   return;
 }
+*/
 
 /********************************************************************
  *********************     Local Routines   *************************