Fix an error in the internal arrays for when changing from constrained/SandR to evolved.

Add a load more debugging statements (ifdef'd).


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

1 files changed, 19 insertions, 4 deletions

diff --git a/src/ChangeType.c b/src/ChangeType.c
index 966445e..9361f3a 100644
--- a/src/ChangeType.c
+++ b/src/ChangeType.c
@@ -164,6 +164,11 @@ CCTK_INT MoL_ChangeToEvolved(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
         EvolvedVariableIndex[MoLNumEvolvedVariables] = EvolvedIndex;
         RHSVariableIndex[MoLNumEvolvedVariables] = RHSIndex;
         MoLNumEvolvedVariables++;
+#ifdef MOLDEBUG
+        printf("Changing (unknown): vars %d var %d (%s).\n",
+               MoLNumEvolvedVariables, EvolvedIndex,
+               CCTK_VarName(EvolvedIndex));
+#endif
         break;
       }
       
@@ -187,9 +192,14 @@ CCTK_INT MoL_ChangeToEvolved(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
           ConstrainedVariableIndex[index] = ConstrainedVariableIndex[index+1];
         }
         MoLNumConstrainedVariables--;
-        EvolvedVariableIndex[index] = EvolvedIndex;
-        RHSVariableIndex[index] = RHSIndex;
+        EvolvedVariableIndex[MoLNumEvolvedVariables] = EvolvedIndex;
+        RHSVariableIndex[MoLNumEvolvedVariables] = RHSIndex;
         MoLNumEvolvedVariables++;
+#ifdef MOLDEBUG
+        printf("Changing (constrained): vars %d var %d (%s).\n",
+               MoLNumEvolvedVariables, EvolvedIndex,
+               CCTK_VarName(EvolvedIndex));
+#endif
         break;
       }
       
@@ -206,9 +216,14 @@ CCTK_INT MoL_ChangeToEvolved(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex)
           SandRVariableIndex[index] = SandRVariableIndex[index+1];
         }
         MoLNumSandRVariables--;
-        EvolvedVariableIndex[index] = EvolvedIndex;
-        RHSVariableIndex[index] = RHSIndex;
+        EvolvedVariableIndex[MoLNumEvolvedVariables] = EvolvedIndex;
+        RHSVariableIndex[MoLNumEvolvedVariables] = RHSIndex;
         MoLNumEvolvedVariables++;
+#ifdef MOLDEBUG
+        printf("Changing (SandR): vars %d var %d (%s).\n",
+               MoLNumEvolvedVariables, EvolvedIndex,
+               CCTK_VarName(EvolvedIndex));
+#endif
         break;
       }