#include "surf_private.h"
#include "xbt/RngStream.h"
#include <math.h>
+#include <unordered_map>
XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_trace, surf, "Surf trace management");
-static xbt_dict_t trace_list = NULL;
+static std::unordered_map<const char *, simgrid::trace_mgr::trace*> trace_list;
-simgrid::trace_mgr::future_evt_set::future_evt_set() {
-}
-
-simgrid::trace_mgr::future_evt_set::~future_evt_set()
-{
- xbt_heap_free(p_heap);
-}
-
-
-/**
- * \brief Create a #tmgr_trace_t from probabilist generators
- *
- * This trace will generate an infinite set of events.
- * It needs two #probabilist_event_generator_t. The date when the event are
- * triggered is directed by date_generator, and will be interpreted as seconds.
- * The value of the event is set by value_generator. The value should be between
- * 0 and 1.
- *
- * \param id The name of the trace
- * \param date_generator The #probabilist_event_generator_t which generates the time
- * between two events
- * \param generator2 The #probabilist_event_generator_t which generates the value
- * of each events.
- * \return The new #tmgr_trace_t
- */
-tmgr_trace_t tmgr_trace_generator_value(const char *id,
- probabilist_event_generator_t date_generator,
- probabilist_event_generator_t value_generator)
-{
- tmgr_trace_t trace = NULL;
-
- trace = xbt_new0(s_tmgr_trace_t, 1);
- trace->type = e_trace_probabilist;
-
- trace->s_probabilist.event_generator[0] = date_generator;
- trace->s_probabilist.event_generator[1] = value_generator;
- trace->s_probabilist.is_state_trace = 0;
-
- return trace;
-}
-
-/**
- * \brief Create a #tmgr_trace_t from probabilist generators
- *
- * This trace will generate an infinite set of events. Value of the events
- * will be alternatively 0 and 1, so this should be used as a state trace.
- *
- * \param id The name of the trace
- * \param date_generator The #probabilist_event_generator_t which generates the time
- * between two events
- * \param first_event_value Set the first event value
- * \return The new #tmgr_trace_t
- */
-tmgr_trace_t tmgr_trace_generator_state(const char *id,
- probabilist_event_generator_t date_generator,
- int first_event_hostIsOn)
-{
- tmgr_trace_t trace = NULL;
-
- trace = xbt_new0(s_tmgr_trace_t, 1);
- trace->type = e_trace_probabilist;
-
- trace->s_probabilist.event_generator[0] = date_generator;
- trace->s_probabilist.event_generator[1] = date_generator;
- trace->s_probabilist.is_state_trace = 1;
- trace->s_probabilist.next_event = first_event_hostIsOn;
-
- return trace;
-}
-
-/**
- * \brief Create a #tmgr_trace_t from probabilist generators
- *
- * This trace will generate an infinite set of events. Value of the events
- * will be alternatively 0 and 1, so this should be used as a state trace.
- *
- * \param id The name of the trace
- * \param avail_duration_generator The #probabilist_event_generator_t which
- * set the duration of the available state, (ie 1 value)
- * \param unavail_duration_generator The #probabilist_event_generator_t which
- * set the duration of the unavailable state, (ie 0 value)
- * \param first_event_value Set the first event value
- * \return The new #tmgr_trace_t
- */
-tmgr_trace_t tmgr_trace_generator_avail_unavail(const char *id,
- probabilist_event_generator_t avail_duration_generator,
- probabilist_event_generator_t unavail_duration_generator,
- int first_event_hostIsOn)
+simgrid::trace_mgr::trace::trace()
{
- tmgr_trace_t trace = NULL;
-
- trace = xbt_new0(s_tmgr_trace_t, 1);
- trace->type = e_trace_probabilist;
-
- trace->s_probabilist.event_generator[0] = unavail_duration_generator;
- trace->s_probabilist.event_generator[1] = avail_duration_generator;
- trace->s_probabilist.is_state_trace = 1;
- trace->s_probabilist.next_event = first_event_hostIsOn;
-
- return trace;
+ event_list = xbt_dynar_new(sizeof(s_tmgr_event_t), NULL);
}
-/**
- * \brief Create a new #probabilist_event_generator_t following the uniform distribution
- *
- * This generator will generate uniformly distributed random values between min and max
- * The id is important : it controls the seed of the generator. So, generators with the
- * same id and the same parameters will generate the same values.
- *
- * \param id The name of the generator
- * \param min The minimal generated value
- * \param max The maximal generated value
- * \return a new #probabilist_event_generator_t
- */
-probabilist_event_generator_t tmgr_event_generator_new_uniform(const char* id,
- double min,
- double max)
+simgrid::trace_mgr::trace::~trace()
{
- probabilist_event_generator_t event_generator = NULL;
- RngStream rng_stream = NULL;
-
- rng_stream = sg_platf_rng_stream_get(id);
-
- event_generator = xbt_new0(s_probabilist_event_generator_t, 1);
- event_generator->type = e_generator_uniform;
- event_generator->s_uniform_parameters.min = min;
- event_generator->s_uniform_parameters.max = max;
- event_generator->rng_stream = rng_stream;
-
- tmgr_event_generator_next_value(event_generator);
-
- return event_generator;
+ xbt_dynar_free(&event_list);
}
-
-
-/**
- * \brief Create a new #probabilist_event_generator_t following the exponential distribution
- *
- * This generator will generate random values following the exponential distribution.
- * The mean value is 1/rate .
- * The id is important : it controls the seed of the generator. So, generators with the
- * same id and the same parameters will generate the same values.
- *
- * \param id The name of the generator
- * \param rate The rate parameter
- * \return a new #probabilist_event_generator_t
- */
-probabilist_event_generator_t tmgr_event_generator_new_exponential(const char* id,
- double rate)
+simgrid::trace_mgr::future_evt_set::future_evt_set()
{
- probabilist_event_generator_t event_generator = NULL;
- RngStream rng_stream = NULL;
-
- rng_stream = sg_platf_rng_stream_get(id);
-
- event_generator = xbt_new0(s_probabilist_event_generator_t, 1);
- event_generator->type = e_generator_exponential;
- event_generator->s_exponential_parameters.rate = rate;
- event_generator->rng_stream = rng_stream;
-
- tmgr_event_generator_next_value(event_generator);
-
- return event_generator;
}
-/**
- * \brief Create a new #probabilist_event_generator_t following the weibull distribution
- *
- * This generator will generate random values following the weibull distribution.
- * The id is important : it controls the seed of the generator. So, generators with the
- * same id and the same parameters will generate the same values.
- *
- * \param id The name of the generator
- * \param scale The scale parameter
- * \param shape The shape parameter
- * \return a new #probabilist_event_generator_t
- */
-probabilist_event_generator_t tmgr_event_generator_new_weibull(const char* id,
- double scale,
- double shape)
-{
- probabilist_event_generator_t event_generator = NULL;
- RngStream rng_stream = NULL;
-
- rng_stream = sg_platf_rng_stream_get(id);
-
- event_generator = xbt_new0(s_probabilist_event_generator_t, 1);
- event_generator->type = e_generator_weibull;
- event_generator->s_weibull_parameters.scale = scale;
- event_generator->s_weibull_parameters.shape = shape;
- event_generator->rng_stream = rng_stream;
-
- tmgr_event_generator_next_value(event_generator);
-
- return event_generator;
-}
-/**
- * \brief Get the next random value of a #probabilist_event_generator_t
- * \param generator The #probabilist_event_generator_t
- * \return the next random value
- */
-double tmgr_event_generator_next_value(probabilist_event_generator_t generator)
+simgrid::trace_mgr::future_evt_set::~future_evt_set()
{
-
- switch(generator->type) {
- case e_generator_uniform:
- generator->next_value = (RngStream_RandU01(generator->rng_stream)
- * (generator->s_uniform_parameters.max - generator->s_uniform_parameters.min))
- + generator->s_uniform_parameters.min;
- break;
- case e_generator_exponential:
- generator->next_value = -log(RngStream_RandU01(generator->rng_stream))
- / generator->s_exponential_parameters.rate;
- break;
- case e_generator_weibull:
- generator->next_value = generator->s_weibull_parameters.scale
- * pow( -log(RngStream_RandU01(generator->rng_stream)),
- 1.0 / generator->s_weibull_parameters.shape );
- }
-
- return generator->next_value;
+ xbt_heap_free(p_heap);
}
-tmgr_trace_t tmgr_trace_new_from_string(const char *id, const char *input,
- double periodicity)
+tmgr_trace_t tmgr_trace_new_from_string(const char *name, const char *input, double periodicity)
{
- tmgr_trace_t trace = NULL;
int linecount = 0;
- s_tmgr_event_t event;
tmgr_event_t last_event = NULL;
- xbt_dynar_t list;
unsigned int cpt;
char *val;
- if (trace_list) {
- trace = (tmgr_trace_t)xbt_dict_get_or_null(trace_list, id);
- if (trace) {
- XBT_WARN("Ignoring redefinition of trace %s", id);
- return trace;
- }
- }
+ xbt_assert(trace_list.find(name) == trace_list.end(), "Refusing to define trace %s twice", name);
+ xbt_assert(periodicity >= 0, "Invalid periodicity %g (must be positive)", periodicity);
- xbt_assert(periodicity >= 0,
- "Invalid periodicity %g (must be positive)", periodicity);
-
- trace = xbt_new0(s_tmgr_trace_t, 1);
- trace->type = e_trace_list;
- trace->s_list.event_list = xbt_dynar_new(sizeof(s_tmgr_event_t), NULL);
-
- list = xbt_str_split(input, "\n\r");
+ tmgr_trace_t trace = new simgrid::trace_mgr::trace();
+ xbt_dynar_t list = xbt_str_split(input, "\n\r");
xbt_dynar_foreach(list, cpt, val) {
+ s_tmgr_event_t event;
linecount++;
xbt_str_trim(val, " \t\n\r\x0B");
- if (val[0] == '#' || val[0] == '\0' || val[0] == '%')
+ if (val[0] == '#' || val[0] == '\0' || val[0] == '%') // pass comments
continue;
if (sscanf(val, "PERIODICITY " "%lg" "\n", &periodicity) == 1)
continue;
- if (sscanf(val, "%lg" " " "%lg" "\n", &event.delta, &event.value) != 2)
- xbt_die("%s:%d: Syntax error in trace\n%s", id, linecount, input);
+ xbt_assert(sscanf(val, "%lg" " " "%lg" "\n", &event.delta, &event.value) == 2,
+ "%s:%d: Syntax error in trace\n%s", name, linecount, input);
if (last_event) {
- if (last_event->delta > event.delta) {
- xbt_die("%s:%d: Invalid trace: Events must be sorted, "
- "but time %g > time %g.\n%s",
- id, linecount, last_event->delta, event.delta, input);
- }
+ xbt_assert(last_event->delta <= event.delta,
+ "%s:%d: Invalid trace: Events must be sorted, but time %g > time %g.\n%s",
+ name, linecount, last_event->delta, event.delta, input);
+
last_event->delta = event.delta - last_event->delta;
} else {
if(event.delta > 0.0){
s_tmgr_event_t first_event;
first_event.delta=event.delta;
first_event.value=-1.0;
- xbt_dynar_push(trace->s_list.event_list, &first_event);
+ xbt_dynar_push(trace->event_list, &first_event);
}
}
- xbt_dynar_push(trace->s_list.event_list, &event);
- last_event = (tmgr_event_t)xbt_dynar_get_ptr(trace->s_list.event_list, xbt_dynar_length(trace->s_list.event_list) - 1);
+ xbt_dynar_push(trace->event_list, &event);
+ last_event = (tmgr_event_t)xbt_dynar_get_ptr(trace->event_list, xbt_dynar_length(trace->event_list) - 1);
}
if (last_event)
last_event->delta = periodicity;
- if (!trace_list)
- trace_list = xbt_dict_new_homogeneous((void (*)(void *)) tmgr_trace_free);
-
- xbt_dict_set(trace_list, id, (void *) trace, NULL);
+ trace_list.insert({xbt_strdup(name), trace});
xbt_dynar_free(&list);
return trace;
tmgr_trace_t tmgr_trace_new_from_file(const char *filename)
{
- tmgr_trace_t trace = NULL;
-
- if ((!filename) || (strcmp(filename, "") == 0))
- return NULL;
-
- if (trace_list) {
- trace = (tmgr_trace_t)xbt_dict_get_or_null(trace_list, filename);
- if (trace) {
- XBT_WARN("Ignoring redefinition of trace %s", filename);
- return trace;
- }
- }
+ xbt_assert(filename && filename[0], "Cannot parse a trace from the null or empty filename");
+ xbt_assert(trace_list.find(filename) == trace_list.end(), "Refusing to define trace %s twice", filename);
FILE *f = surf_fopen(filename, "r");
xbt_assert(f != NULL,
char *tstr = xbt_str_from_file(f);
fclose(f);
- trace = tmgr_trace_new_from_string(filename, tstr, 0.);
+ tmgr_trace_t trace = tmgr_trace_new_from_string(filename, tstr, 0.);
xbt_free(tstr);
return trace;
tmgr_trace_t tmgr_empty_trace_new(void)
{
- tmgr_trace_t trace = NULL;
+ tmgr_trace_t trace = new simgrid::trace_mgr::trace();
s_tmgr_event_t event;
-
- trace = xbt_new0(s_tmgr_trace_t, 1);
- trace->type = e_trace_list;
- trace->s_list.event_list = xbt_dynar_new(sizeof(s_tmgr_event_t), NULL);
-
event.delta = 0.0;
event.value = 0.0;
- xbt_dynar_push(trace->s_list.event_list, &event);
+ xbt_dynar_push(trace->event_list, &event);
return trace;
}
void tmgr_trace_free(tmgr_trace_t trace)
{
- if (!trace)
- return;
-
- switch(trace->type) {
- case e_trace_list:
- xbt_dynar_free(&(trace->s_list.event_list));
- break;
- case e_trace_probabilist:
- THROW_UNIMPLEMENTED;
- break;
- }
- free(trace);
+ delete trace;
}
/** @brief Registers a new trace into the future event set, and get an iterator over the integrated trace */
-tmgr_trace_iterator_t simgrid::trace_mgr::future_evt_set::add_trace(
- tmgr_trace_t trace, double start_time, surf::Resource *resource)
+tmgr_trace_iterator_t simgrid::trace_mgr::future_evt_set::add_trace(tmgr_trace_t trace, double start_time, surf::Resource *resource)
{
tmgr_trace_iterator_t trace_iterator = NULL;
trace_iterator->idx = 0;
trace_iterator->resource = resource;
- if(trace->type == e_trace_list) {
- xbt_assert((trace_iterator->idx < xbt_dynar_length(trace->s_list.event_list)),
- "You're referring to an event that does not exist!");
- }
+ xbt_assert((trace_iterator->idx < xbt_dynar_length(trace->event_list)), "Your trace should have at least one event!");
xbt_heap_push(p_heap, trace_iterator, start_time);
}
/** @brief returns the date of the next occurring event (pure function) */
-double simgrid::trace_mgr::future_evt_set::next_date()
+double simgrid::trace_mgr::future_evt_set::next_date() const
{
if (xbt_heap_size(p_heap))
return (xbt_heap_maxkey(p_heap));
double date, double *value, simgrid::surf::Resource **resource)
{
double event_date = next_date();
- tmgr_trace_iterator_t trace_iterator = NULL;
- tmgr_event_t event = NULL;
- tmgr_trace_t trace = NULL;
- double event_delta;
-
if (event_date > date)
return NULL;
- if (!(trace_iterator = (tmgr_trace_iterator_t)xbt_heap_pop(p_heap)))
+ tmgr_trace_iterator_t trace_iterator = (tmgr_trace_iterator_t)xbt_heap_pop(p_heap);
+ if (trace_iterator == NULL)
return NULL;
- trace = trace_iterator->trace;
+ tmgr_trace_t trace = trace_iterator->trace;
*resource = trace_iterator->resource;
- switch(trace->type) {
- case e_trace_list:
-
- event = (tmgr_event_t)xbt_dynar_get_ptr(trace->s_list.event_list, trace_iterator->idx);
+ tmgr_event_t event = (tmgr_event_t)xbt_dynar_get_ptr(trace->event_list, trace_iterator->idx);
- *value = event->value;
+ *value = event->value;
- if (trace_iterator->idx < xbt_dynar_length(trace->s_list.event_list) - 1) {
- xbt_heap_push(p_heap, trace_iterator, event_date + event->delta);
- trace_iterator->idx++;
- } else if (event->delta > 0) { /* Last element, checking for periodicity */
- xbt_heap_push(p_heap, trace_iterator, event_date + event->delta);
- trace_iterator->idx = 1; /* not 0 as the first event is a placeholder to handle when events really start */
- } else { /* We don't need this trace_event anymore */
- trace_iterator->free_me = 1;
- }
- break;
-
- case e_trace_probabilist:
-
- //FIXME : not tested yet
- if(trace->s_probabilist.is_state_trace) {
- *value = (double) trace->s_probabilist.next_event;
- if(trace->s_probabilist.next_event == 0) {
- event_delta = tmgr_event_generator_next_value(trace->s_probabilist.event_generator[0]);
- trace->s_probabilist.next_event = 1;
- } else {
- event_delta = tmgr_event_generator_next_value(trace->s_probabilist.event_generator[1]);
- trace->s_probabilist.next_event = 0;
- }
- } else {
- event_delta = tmgr_event_generator_next_value(trace->s_probabilist.event_generator[0]);
- *value = tmgr_event_generator_next_value(trace->s_probabilist.event_generator[1]);
- }
- xbt_heap_push(p_heap, trace_iterator, event_date + event_delta);
- XBT_DEBUG("Generating a new event at date %f, with value %f", event_date + event_delta, *value);
-
- break;
+ if (trace_iterator->idx < xbt_dynar_length(trace->event_list) - 1) {
+ xbt_heap_push(p_heap, trace_iterator, event_date + event->delta);
+ trace_iterator->idx++;
+ } else if (event->delta > 0) { /* Last element, checking for periodicity */
+ xbt_heap_push(p_heap, trace_iterator, event_date + event->delta);
+ trace_iterator->idx = 1; /* not 0 as the first event is a placeholder to handle when events really start */
+ } else { /* We don't need this trace_event anymore */
+ trace_iterator->free_me = 1;
}
return trace_iterator;
void tmgr_finalize(void)
{
- xbt_dict_free(&trace_list);
+ for (auto kv : trace_list) {
+ xbt_free((char*)kv.first);
+ delete kv.second;
+ }
}
void tmgr_trace_event_unref(tmgr_trace_iterator_t *trace_event)