X-Git-Url: http://info.iut-bm.univ-fcomte.fr/pub/gitweb/simgrid.git/blobdiff_plain/a9f3bc989ac7b641f38385167dbd3063e209d88e..693f30b46244c152cd79cdf3ad35d4a79b866c9c:/src/surf/trace_mgr.c diff --git a/src/surf/trace_mgr.c b/src/surf/trace_mgr.c index b3db0851e1..c541762c51 100644 --- a/src/surf/trace_mgr.c +++ b/src/surf/trace_mgr.c @@ -36,94 +36,108 @@ XBT_INLINE void tmgr_history_free(tmgr_history_t h) tmgr_trace_t tmgr_trace_new_from_generator(const char *id, probabilist_event_generator_t generator1, - probabilist_event_generator_t generator2) + probabilist_event_generator_t generator2, + int is_state_trace) { tmgr_trace_t trace = NULL; - RngStream rng_stream = NULL; - - rng_stream = sg_platf_rng_stream_get(id); - + trace = xbt_new0(s_tmgr_trace_t, 1); trace->type = e_trace_probabilist; - + trace->s_probabilist.event_generator[0] = generator1; - trace->s_probabilist.event_generator[0]->rng_stream = rng_stream; - tmgr_event_generator_next_value(trace->s_probabilist.event_generator[0]); - + //FIXME : may also be a parameter trace->s_probabilist.next_event = 0; - - if(generator2 == NULL) { + trace->s_probabilist.is_state_trace = is_state_trace; + + if(generator2 != NULL) { + trace->s_probabilist.event_generator[1] = generator2; + } else if(is_state_trace) { trace->s_probabilist.event_generator[1] = generator1; } else { - trace->s_probabilist.event_generator[1] = generator2; - trace->s_probabilist.event_generator[1]->rng_stream = rng_stream; - tmgr_event_generator_next_value(trace->s_probabilist.event_generator[1]); + THROW_IMPOSSIBLE; //That case should have been checked before, anyway... } - + return trace; } -probabilist_event_generator_t tmgr_event_generator_new_uniform(double alpha, - double beta) -{ +probabilist_event_generator_t tmgr_event_generator_new_uniform(const char* id, + double min, + double max) +{ 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.alpha = alpha; - event_generator->s_uniform_parameters.beta = beta; + 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; } -probabilist_event_generator_t tmgr_event_generator_new_exponential(double lambda) -{ +probabilist_event_generator_t tmgr_event_generator_new_exponential(const char* id, + double rate) +{ 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.lambda = lambda; + event_generator->s_exponential_parameters.rate = rate; + event_generator->rng_stream = rng_stream; + + tmgr_event_generator_next_value(event_generator); return event_generator; } -probabilist_event_generator_t tmgr_event_generator_new_weibull(double lambda, - double k) -{ +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.lambda = lambda; - event_generator->s_weibull_parameters.k = k; + 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; } double tmgr_event_generator_next_value(probabilist_event_generator_t generator) { - + switch(generator->type) { case e_generator_uniform: generator->next_value = (RngStream_RandU01(generator->rng_stream) - * (generator->s_uniform_parameters.beta - generator->s_uniform_parameters.alpha)) - + generator->s_uniform_parameters.alpha; + * (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.lambda; + / generator->s_exponential_parameters.rate; break; case e_generator_weibull: - generator->next_value = - generator->s_weibull_parameters.lambda + generator->next_value = - generator->s_weibull_parameters.scale * pow( log(RngStream_RandU01(generator->rng_stream)), - 1.0 / generator->s_weibull_parameters.k ); + 1.0 / generator->s_weibull_parameters.shape ); } - + return generator->next_value; } @@ -241,7 +255,7 @@ XBT_INLINE 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)); @@ -269,7 +283,7 @@ tmgr_trace_event_t tmgr_history_add_trace(tmgr_history_t h, xbt_assert((trace_event->idx < xbt_dynar_length(trace->s_list.event_list)), "You're referring to an event that does not exist!"); } - + xbt_heap_push(h->heap, trace_event, start_time); return trace_event; @@ -302,10 +316,10 @@ tmgr_trace_event_t tmgr_history_get_next_event_leq(tmgr_history_t h, trace = trace_event->trace; *model = trace_event->model; - + switch(trace->type) { case e_trace_list: - + event = xbt_dynar_get_ptr(trace->s_list.event_list, trace_event->idx); *value = event->value; @@ -320,21 +334,25 @@ tmgr_trace_event_t tmgr_history_get_next_event_leq(tmgr_history_t h, trace_event->free_me = 1; } break; - + case e_trace_probabilist: - - //FIXME : Simplist and not tested - //Would only work for failure events for now - *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 = 0; + + //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 = 0; + } else { + event_delta = tmgr_event_generator_next_value(trace->s_probabilist.event_generator[1]); + 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 = 1; + 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(h->heap, trace_event, event_date + event_delta); - + break; }