diff options
-rw-r--r-- | src/InitialCopy.c | 31 |
1 files changed, 31 insertions, 0 deletions
diff --git a/src/InitialCopy.c b/src/InitialCopy.c index c4cfbf4..650f248 100644 --- a/src/InitialCopy.c +++ b/src/InitialCopy.c @@ -144,6 +144,37 @@ void MoL_InitialCopy(CCTK_ARGUMENTS) MoL_LinearCombination(cctkGH, EvolvedVariableIndex[var], rl, tl, 0.0, srcs, tls, facts, nsrc); } + + for (var = 0; var < MoLNumEvolvedVariablesSlow; var++) + { + const int nsrc = 1; + const int srcs[1] = {EvolvedVariableIndexSlow[var]}; + const int tls[1] = {1}; + const CCTK_REAL facts[1] = {1.0}; + + StorageOn = CCTK_QueryGroupStorageI(cctkGH, + CCTK_GroupIndexFromVarI(EvolvedVariableIndexSlow[var])); + + if (StorageOn < 0) + { + CCTK_VWarn(1,__LINE__,__FILE__,CCTK_THORNSTRING,"Warning for index %i", + EvolvedVariableIndexSlow[var]); + CCTK_WARN(0, "The index passed does not correspond to a GF."); + } + else if (StorageOn == 0) { +#ifdef MOLDEBUG + printf("Aargh! Vars %d var %d index %d name %s\n", + MoLNumEvolvedVariablesSlow, var, EvolvedVariableIndexSlow[var], + CCTK_VarName(EvolvedVariableIndexSlow[var])); +#endif + CCTK_VWarn(1,__LINE__,__FILE__,CCTK_THORNSTRING,"Warning for GF %s", + CCTK_VarName(EvolvedVariableIndexSlow[var])); + CCTK_WARN(0, "The grid function does not have storage assigned."); + } + + MoL_LinearCombination(cctkGH, EvolvedVariableIndexSlow[var], rl, tl, 0.0, + srcs, tls, facts, nsrc); + } /* Set up the array sizes */ |