stats.h File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
struct	OutputFD

Functions
void	set_stats_params (ParameterSet *ps)
void	open_outputfiles (int RestartSnapNum, struct OutputFD *fds, const char *OutputDir, int BlackHoleOn, int StarformationOn)
void	close_outputfiles (struct OutputFD *fds)
void	write_cpu_log (int NumCurrentTiStep, const double atime, FILE *FdCPU, double ElapsedTime)
void	energy_statistics (FILE *FdEnergy, const double Time, struct part_manager_type *PartManager)

Function Documentation

close_outputfiles()

void close_outputfiles

(

struct OutputFD *

fds

)

This function closes the global log-files.

Definition at line 150 of file stats.c.

{
    if(fds->FdCPU)
        fclose(fds->FdCPU);
    if(fds->FdEnergy)
        fclose(fds->FdEnergy);
    if(fds->FdSfr)
        fclose(fds->FdSfr);
    if(fds->FdBlackHoles)
        fclose(fds->FdBlackHoles);
    if(fds->FdBlackholeDetails)
        fclose(fds->FdBlackholeDetails);
}

References OutputFD::FdBlackholeDetails, OutputFD::FdBlackHoles, OutputFD::FdCPU, OutputFD::FdEnergy, and OutputFD::FdSfr.

Referenced by run().

Here is the caller graph for this function:

energy_statistics()

void energy_statistics	(	FILE *	FdEnergy,
		const double	Time,
		struct part_manager_type *	PartManager
	)

This routine first calls a computation of various global quantities of the particle distribution, and then writes some statistics about the energies in the various particle components to the file FdEnergy.

Definition at line 334 of file stats.c.

{
    if(!FdEnergy)
        return;
 
    struct state_of_system SysState = compute_global_quantities_of_system(Time, PartManager);
 
    message(0, "Time %g Mean Temperature of Gas %g\n",
                Time, SysState.TemperatureComp[0]);
 
    fprintf(FdEnergy,
            "%g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g %g\n",
            Time, SysState.TemperatureComp[0], SysState.EnergyInt, SysState.EnergyPot, SysState.EnergyKin, SysState.EnergyIntComp[0],
            SysState.EnergyPotComp[0], SysState.EnergyKinComp[0], SysState.EnergyIntComp[1],
            SysState.EnergyPotComp[1], SysState.EnergyKinComp[1], SysState.EnergyIntComp[2],
            SysState.EnergyPotComp[2], SysState.EnergyKinComp[2], SysState.EnergyIntComp[3],
            SysState.EnergyPotComp[3], SysState.EnergyKinComp[3], SysState.EnergyIntComp[4],
            SysState.EnergyPotComp[4], SysState.EnergyKinComp[4], SysState.EnergyIntComp[5],
            SysState.EnergyPotComp[5], SysState.EnergyKinComp[5], SysState.MassComp[0],
            SysState.MassComp[1], SysState.MassComp[2], SysState.MassComp[3], SysState.MassComp[4],
            SysState.MassComp[5]);
 
    fflush(FdEnergy);
}

References compute_global_quantities_of_system(), state_of_system::EnergyInt, state_of_system::EnergyIntComp, state_of_system::EnergyKin, state_of_system::EnergyKinComp, state_of_system::EnergyPot, state_of_system::EnergyPotComp, state_of_system::MassComp, message(), PartManager, and state_of_system::TemperatureComp.

Referenced by run().

Here is the call graph for this function:

Here is the caller graph for this function:

open_outputfiles()

void open_outputfiles	(	int	RestartSnapNum,
		struct OutputFD *	fds,
		const char *	OutputDir,
		int	BlackHoleOn,
		int	StarformationOn
	)

This function opens various log-files that report on the status and performance of the simulstion. On restart from restart-files (start-option 1), the code will append to these files.

Definition at line 73 of file stats.c.

{
    const char mode[3]="a+";
    char * buf;
    char * postfix;
    int ThisTask;
    MPI_Comm_rank(MPI_COMM_WORLD, &ThisTask);
    memset(fds, 0, sizeof(struct OutputFD));
    fds->FdCPU = NULL;
    fds->FdEnergy = NULL;
    fds->FdBlackHoles = NULL;
    fds->FdSfr = NULL;
    fds->FdBlackholeDetails = NULL;
    fds->FdHelium = NULL;
 
    if(RestartSnapNum != -1) {
        postfix = fastpm_strdup_printf("-R%03d", RestartSnapNum);
    } else {
        postfix = fastpm_strdup_printf("%s", "");
    }
 
    /* all the processors write to separate files*/
    if(BlackHoleOn && StatsParams.WriteBlackHoleDetails){
        buf = fastpm_strdup_printf("%s/%s%s/%06X", OutputDir,"BlackholeDetails",postfix,ThisTask);
        fastpm_path_ensure_dirname(buf);
        if(!(fds->FdBlackholeDetails = fopen(buf,"a")))
            endrun(1, "Failed to open blackhole detail %s\n", buf);
        myfree(buf);
    }
 
    /* only the root processors writes to the log files */
    if(ThisTask != 0) {
        return;
    }
 
    if(BlackHoleOn) {
        buf = fastpm_strdup_printf("%s/%s%s", OutputDir, "blackholes.txt", postfix);
        fastpm_path_ensure_dirname(buf);
        if(!(fds->FdBlackHoles = fopen(buf, mode)))
            endrun(1, "error in opening file '%s'\n", buf);
        myfree(buf);
    }
 
    buf = fastpm_strdup_printf("%s/%s%s", OutputDir, StatsParams.CpuFile, postfix);
    fastpm_path_ensure_dirname(buf);
    if(!(fds->FdCPU = fopen(buf, mode)))
        endrun(1, "error in opening file '%s'\n", buf);
    myfree(buf);
 
    if(StatsParams.OutputEnergyDebug) {
        buf = fastpm_strdup_printf("%s/%s%s", OutputDir, StatsParams.EnergyFile, postfix);
        fastpm_path_ensure_dirname(buf);
        if(!(fds->FdEnergy = fopen(buf, mode)))
            endrun(1, "error in opening file '%s'\n", buf);
        myfree(buf);
    }
 
    if(StarformationOn) {
        buf = fastpm_strdup_printf("%s/%s%s", OutputDir, "sfr.txt", postfix);
        fastpm_path_ensure_dirname(buf);
        if(!(fds->FdSfr = fopen(buf, mode)))
            endrun(1, "error in opening file '%s'\n", buf);
        myfree(buf);
    }
 
    if(qso_lightup_on()) {
        buf = fastpm_strdup_printf("%s/%s%s", OutputDir, "helium.txt", postfix);
        fastpm_path_ensure_dirname(buf);
        if(!(fds->FdHelium = fopen(buf, mode)))
            endrun(1, "error in opening file '%s'\n", buf);
        myfree(buf);
    }
}

References stats_params::CpuFile, endrun(), stats_params::EnergyFile, fastpm_path_ensure_dirname(), fastpm_strdup_printf(), OutputFD::FdBlackholeDetails, OutputFD::FdBlackHoles, OutputFD::FdCPU, OutputFD::FdEnergy, OutputFD::FdHelium, OutputFD::FdSfr, myfree, stats_params::OutputEnergyDebug, qso_lightup_on(), StatsParams, ThisTask, and stats_params::WriteBlackHoleDetails.

Referenced by run().

Here is the call graph for this function:

Here is the caller graph for this function:

set_stats_params()

void set_stats_params

(

ParameterSet *

ps

)

Definition at line 55 of file stats.c.

{
    int ThisTask;
    MPI_Comm_rank(MPI_COMM_WORLD, &ThisTask);
    if(ThisTask == 0) {
        param_get_string2(ps, "EnergyFile", StatsParams.EnergyFile, sizeof(StatsParams.EnergyFile));
        param_get_string2(ps, "CpuFile", StatsParams.CpuFile, sizeof(StatsParams.CpuFile));
        StatsParams.OutputEnergyDebug = param_get_int(ps, "OutputEnergyDebug");
        StatsParams.WriteBlackHoleDetails = param_get_int(ps,"WriteBlackHoleDetails");
    }
    MPI_Bcast(&StatsParams, sizeof(struct stats_params), MPI_BYTE, 0, MPI_COMM_WORLD);
}

References stats_params::CpuFile, stats_params::EnergyFile, stats_params::OutputEnergyDebug, param_get_int(), param_get_string2(), StatsParams, ThisTask, and stats_params::WriteBlackHoleDetails.

Referenced by read_parameter_file().

Here is the call graph for this function:

Here is the caller graph for this function:

write_cpu_log()

void write_cpu_log	(	int	NumCurrentTiStep,
		const double	atime,
		FILE *	FdCPU,
		double	ElapsedTime
	)

Definition at line 165 of file stats.c.

{
    walltime_summary(0, MPI_COMM_WORLD);
 
    if(FdCPU)
    {
        int NTask;
        MPI_Comm_size(MPI_COMM_WORLD, &NTask);
        fprintf(FdCPU, "Step %d, Time: %g, MPIs: %d Threads: %d Elapsed: %g\n", NumCurrentTiStep, atime, NTask, omp_get_max_threads(), ElapsedTime);
        walltime_report(FdCPU, 0, MPI_COMM_WORLD);
        fflush(FdCPU);
    }
}

Referenced by run().

Here is the caller graph for this function: