/* check if steerable parameters have changed */
static void CheckSteerableParameters (asciiioGH *myGH)
{
int i, num_vars, out_old, times_set;
int verboselength;
char *fullname, *msg;
static int out1D_vars_lastset = -1;
DECLARE_CCTK_PARAMETERS
out_old = myGH->out1D_every;
/* How often to output */
myGH->out1D_every = out_every > 0 ? out_every : -1;
if (out1D_every > 0)
{
myGH->out1D_every = out1D_every;
}
/* Report if frequency changed */
if (myGH->out1D_every != out_old)
{
if (CCTK_Equals (newverbose, "standard") ||
CCTK_Equals (newverbose, "full"))
{
CCTK_VInfo (CCTK_THORNSTRING, "IOASCII_1D: Output every %d iterations",
myGH->out1D_every);
}
}
/* re-parse the 'out1D_vars' parameter if it was changed */
times_set = CCTK_ParameterQueryTimesSet ("out1D_vars", CCTK_THORNSTRING);
if (times_set != out1D_vars_lastset)
{
num_vars = CCTK_NumVars ();
memset (myGH->do_out1D, 0, num_vars);
CCTK_TraverseString (out1D_vars, SetOutputFlag, myGH->do_out1D,
CCTK_GROUP_OR_VAR);
if (CCTK_Equals (newverbose, "standard") ||
CCTK_Equals (newverbose, "full"))
{
/* Count the length of the string */
verboselength = 0;
for (i = 0; i < num_vars; i++)
{
if (myGH->do_out1D[i])
{
fullname = CCTK_FullName (i);
verboselength += strlen(fullname)+1;
free (fullname);
}
}
verboselength += strlen("IOASCII_1D: Output requested for ");
msg = (char *)malloc((verboselength+1)*sizeof(char));
sprintf(msg,"IOASCII_1D: Output requested for ");
for (i = 0; i < num_vars; i++)
{
if (myGH->do_out1D[i])
{
fullname = CCTK_FullName (i);
sprintf (msg, "%s %s",msg,fullname);
free (fullname);
}
}
if (msg)
{
CCTK_INFO (msg);
free (msg);
}
}
/* Save the last setting of 'out1D_vars' parameter */
out1D_vars_lastset = times_set;
}
USE_CCTK_PARAMETERS; }
/* check if steerable parameters have changed */
static void CheckSteerableParameters (iobasicGH *myGH)
{
int i, num_vars, out_old;
int times_set;
char *fullname, *msg, *oldmsg;
static int outScalar_vars_lastset = -1;
DECLARE_CCTK_PARAMETERS
/* How often to output */
out_old = myGH->outScalar_every;
myGH->outScalar_every = out_every > 0 ? out_every : -1;
if (outScalar_every > 0)
{
myGH->outScalar_every = outScalar_every;
}
if (myGH->outScalar_every != out_old)
{
if (CCTK_Equals (newverbose, "standard") ||
CCTK_Equals (newverbose, "full"))
{
CCTK_VInfo (CCTK_THORNSTRING, "Scalar: Output every %d iterations",
myGH->outScalar_every);
}
}
/* re-parse the 'outScalar_vars' parameter if it was changed */
times_set = CCTK_ParameterQueryTimesSet ("outScalar_vars", CCTK_THORNSTRING);
if (times_set != outScalar_vars_lastset)
{
num_vars = CCTK_NumVars ();
memset (myGH->do_outScalar, 0, num_vars);
CCTK_TraverseString (outScalar_vars, SetOutputFlag, myGH->do_outScalar,
CCTK_GROUP_OR_VAR);
if (myGH->outScalar_every &&
(CCTK_Equals (newverbose, "standard") ||
CCTK_Equals (newverbose, "full")))
{
msg = NULL;
for (i = 0; i < num_vars; i++)
{
if (myGH->do_outScalar[i])
{
fullname = CCTK_FullName (i);
if (! msg)
{
Util_asprintf (&msg, "Scalar: Output requested for %s", fullname);
}
else
{
oldmsg = msg;
Util_asprintf (&msg, "%s %s", oldmsg, fullname);
free (oldmsg);
}
free (fullname);
}
}
if (msg)
{
CCTK_INFO (msg);
free (msg);
}
}
/* Save the last setting of 'outScalar_vars' parameter */
outScalar_vars_lastset = times_set;
}
USE_CCTK_PARAMETERS; }
/*@@
@routine PUGH_DisableGroupStorage
@author Tom Goodale
@date 30 Mar 1999
@desc
Disables storage for all variables in the group
indicated by groupname.
@enddesc
@calls CCTK_GroupIndex
CCTK_GroupData
CCTK_FirstVarIndexI
PUGH_DisableGArrayGroupStorage
@var GH
@vdesc Pointer to CCTK grid hierarchy
@vtype cGH *
@vio in
@endvar
@var groupname
@vdesc name of the group to enable storage for
@vtype const char *
@vio in
@endvar
@returntype int
@returndesc
1 if storage for given group was disabled
-1 if group type is invalid
@endreturndesc
@@*/
int PUGH_DisableGroupStorage (cGH *GH, const char *groupname)
{
DECLARE_CCTK_PARAMETERS
int group; /* group index */
cGroup pgroup; /* group information */
int vtypesize, retval;
pGA ***variables;
int first_var, var, level;
int unchanged; /* count how many aren't toggled */
char *temp;
#ifdef DEBUG_PUGH
printf (" PUGH_DisableGroupStorage: request for group '%s'\n", groupname);
fflush (stdout);
#endif
group = CCTK_GroupIndex (groupname);
CCTK_GroupData (group, &pgroup);
/* get global index of first variable in group */
first_var = CCTK_FirstVarIndexI (group);
variables = (pGA ***) PUGH_pGH (GH)->variables;
/* get the group info from its index */
unchanged = 0;
retval = 1;
if (pgroup.grouptype == CCTK_GF || pgroup.grouptype == CCTK_ARRAY)
{
for (var = first_var; var < first_var+pgroup.numvars; var++)
{
for (level = 0; level < pgroup.numtimelevels; level++)
{
unchanged += PUGH_DisableGArrayDataStorage (variables[var][level]);
}
}
}
else if (pgroup.grouptype == CCTK_SCALAR)
{
vtypesize = CCTK_VarTypeSize (pgroup.vartype);
for (var = first_var; var < first_var+pgroup.numvars; var++)
{
for (level = 0; level < pgroup.numtimelevels; level++)
{
temp = (char *) variables[var][level];
if (temp[vtypesize] == PUGH_STORAGE)
{
temp[vtypesize] = PUGH_NOSTORAGE;
}
else
{
unchanged++;
}
}
}
}
else
{
CCTK_WARN (1, "Unknown group type in PUGH_DisableGroupStorage");
retval = -1;
}
/* Report on memory usage */
if (!CCTK_Equals(storage_verbose,"no") && retval >= 0)
{
if (unchanged == 0)
{
/* Memory toggled */
if (pgroup.grouptype == CCTK_GF)
{
totalnumberGF -= pgroup.numvars;
}
else if (pgroup.grouptype == CCTK_ARRAY)
{
totalnumberGA -= pgroup.numvars;
}
totalstorage -= (variables[first_var][0]->extras->npoints *
variables[first_var][0]->varsize *
pgroup.numtimelevels * pgroup.numvars) /
(float) (1024 * 1024);
if (CCTK_Equals(storage_verbose,"yes"))
{
CCTK_VInfo (CCTK_THORNSTRING, "Switched memory off for group '%s'"
" [GFs: %d GAs: %d Total Size: %6.2fMB]",
groupname, totalnumberGF, totalnumberGA, totalstorage);
}
}
else if (unchanged == pgroup.numvars)
{
/* Memory already off */
if (CCTK_Equals(storage_verbose,"yes"))
{
CCTK_VInfo (CCTK_THORNSTRING, "Memory already off for group '%s'", groupname);
}
}
else
{
CCTK_WARN (1, "PUGH_DisableGroupStorage: Inconsistency in group memory assignment");
}
}
USE_CCTK_PARAMETERS; return (retval);
}
/*@@
@routine PUGHi_PrintStorageReport
@author Gabrielle Allen
@date 16th Sept 2001
@desc
Print a report about the use of storage
@enddesc
@@*/
void PUGHi_PrintStorageReport ()
{
CCTK_INFO("Storage statistics");
CCTK_VInfo(CCTK_THORNSTRING, " Maximum storage: %6.2fMB",maxstorage);
CCTK_VInfo(CCTK_THORNSTRING, " [%d Grid Functions, %d Grid Arrays]",numberGF,numberGA);
}
/*@@
@routine PUGH_EnableGroupStorage
@author Tom Goodale
@date 30 Mar 1999
@desc
Enables storage for all variables in the group
indicated by groupname.
@enddesc
@calls CCTK_GroupIndex
CCTK_GroupData
PUGH_EnableScalarGroupStorage
PUGH_EnableGArrayGroupStorage
@var GH
@vdesc Pointer to CCTK grid hierarchy
@vtype cGH *
@vio in
@endvar
@var groupname
@vdesc name of the group to enable storage for
@vtype const char *
@vio in
@endvar
@returntype int
@returndesc
return code of @seeroutine PUGH_EnableScalarGroupStorage or
@seeroutine PUGH_EnableGArrayGroupStorage: <BR>
1 if storage was already enabled, or <BR>
0 if storage was successfully enabled <BR>
-1 if group type is not one of the above <BR>
-2 if an invalid GH pointer was given <BR>
-3 if invalid groupname was given
@endreturndesc
@@*/
int PUGH_EnableGroupStorage (cGH *GH, const char *groupname)
{
DECLARE_CCTK_PARAMETERS
int group; /* group index */
int first_var; /* first variable's index */
cGroup pgroup; /* group information */
int retval;
pGA *GA;
pGH *pughGH;
#ifdef DEBUG_PUGH
printf (" PUGH_EnableGroupStorage: request for group '%s'\n", groupname);
fflush (stdout);
#endif
pughGH = PUGH_pGH (GH);
group = CCTK_GroupIndex (groupname);
if (pughGH && group >= 0)
{
first_var = CCTK_FirstVarIndexI (group);
/* get the group info from its index */
CCTK_GroupData (group, &pgroup);
if (pgroup.grouptype == CCTK_SCALAR)
{
retval = PUGH_EnableScalarGroupStorage (pughGH,
first_var,
pgroup.numvars,
pgroup.numtimelevels);
}
else if (pgroup.grouptype == CCTK_GF || pgroup.grouptype == CCTK_ARRAY)
{
retval = PUGH_EnableGArrayGroupStorage (pughGH,
first_var,
pgroup.numvars,
pgroup.numtimelevels);
if (!CCTK_Equals(storage_verbose,"no") && retval == 0)
{
/* get GA pointer of first var in group */
GA = (pGA *) pughGH->variables[first_var][0];
if (pgroup.grouptype == CCTK_GF)
{
totalnumberGF += pgroup.numvars * pgroup.numtimelevels;
}
else
{
totalnumberGA += pgroup.numvars * pgroup.numtimelevels;
}
totalstorage += (GA->extras->npoints * GA->varsize *
pgroup.numtimelevels * pgroup.numvars) /
(float) (1024*1024);
if (totalstorage > maxstorage)
{
numberGF = totalnumberGF;
numberGA = totalnumberGA;
maxstorage = totalstorage;
}
/* Report on memory usage */
if (CCTK_Equals(storage_verbose,"yes"))
{
CCTK_VInfo (CCTK_THORNSTRING, "Switched memory on for group '%s'"
" [GFs: %d GAs: %d Total Size: %6.2fMB]",
groupname, totalnumberGF, totalnumberGA, totalstorage);
}
}
}
else
{
CCTK_WARN (1, "PUGH_EnableGroupStorage: Unknown group type");
retval = -1;
}
}
else
{
if (! pughGH)
{
CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING,
"PUGH_EnableGroupStorage: Error with cctkGH pointer "
"for group %s", groupname);
retval = -2;
}
else
{
CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING,
"PUGH_EnableGroupStorage: Invalid group %s", groupname);
retval = -3;
}
}
USE_CCTK_PARAMETERS; return (retval);
}
