timebinmgr.h

Go to the documentation of this file.

#ifndef TIMEBINMGR_H
#define TIMEBINMGR_H
/* This file manages the integer timeline,
 * and converts from integers ti to double loga.*/
 
#define TIMEBINS 20
#define TIMEBASE (1u<<TIMEBINS)
#define MAXSNAPSHOTS (1u<<(30-TIMEBINS))
 
#include "types.h"
#include "utils/paramset.h"
 
typedef struct SyncPoint SyncPoint;
 
struct SyncPoint
{
    double a;
    double loga;
    int write_snapshot;
    int write_fof;
    inttime_t ti;
};
 
 
/*Convert an integer to and from loga*/
double loga_from_ti(inttime_t ti);
inttime_t ti_from_loga(double loga);
 
/*Convert changes in loga to and from ti*/
inttime_t dti_from_dloga(double loga, const inttime_t Ti_Current);
double dloga_from_dti(inttime_t dti, const inttime_t Ti_Current);
 
/*Get dloga from a timebin*/
double get_dloga_for_bin(int timebin, const inttime_t Ti_Current);
 
/*Get the dti from the timebin*/
static inline inttime_t dti_from_timebin(int bin) {
    /*Casts to work around bug in intel compiler 18.0*/
    return bin > 0 ? (1u << (unsigned) bin) : 0;
}
 
/* Enforce that an integer timestep is a power
 * of two subdivision of TIMEBASE, rounding down
 * to the first power of two less than the ti passed in.
 * Note TIMEBASE is the maximum value returned.*/
inttime_t round_down_power_of_two(inttime_t ti);
 
inttime_t find_next_outputtime(inttime_t ti_curr);
 
/*Get whatever is the last output number from ti*/
inttime_t out_from_ti(inttime_t ti);
 
int OutputListAction(ParameterSet * ps, const char * name, void * data);
void set_sync_params(int OutputListLength, double * OutputListTimes);
void setup_sync_points(double TimeIC, double TimeMax, double no_snapshot_until_time, int SnapshotWithFOF);
 
SyncPoint *
find_next_sync_point(inttime_t ti);
 
SyncPoint *
find_current_sync_point(inttime_t ti);
 
SyncPoint *
make_unplanned_sync_point(inttime_t ti);
 
#endif