int is_state_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] = generator1;
-
+
//FIXME : may also be a parameter
trace->s_probabilist.next_event = 0;
trace->s_probabilist.is_state_trace = is_state_trace;
-
+
if(generator2 != NULL) {
trace->s_probabilist.event_generator[1] = generator2;
} else if(is_state_trace) {
} else {
THROW_IMPOSSIBLE; //That case should have been checked before, anyway...
}
-
+
return trace;
}
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.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(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.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(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;
}
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)
* pow( log(RngStream_RandU01(generator->rng_stream)),
1.0 / generator->s_weibull_parameters.shape );
}
-
+
return generator->next_value;
}
{
if (!trace)
return;
-
+
switch(trace->type) {
case e_trace_list:
xbt_dynar_free(&(trace->s_list.event_list));
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;
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;
trace_event->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;
*value = tmgr_event_generator_next_value(trace->s_probabilist.event_generator[1]);
}
xbt_heap_push(h->heap, trace_event, event_date + event_delta);
-
+
break;
}
};
static double A1p76[3][3] = {
- { 82758667.0, 1871391091.0, 4127413238.0 },
- { 3672831523.0, 69195019.0, 1871391091.0 },
+ { 82758667.0, 1871391091.0, 4127413238.0 },
+ { 3672831523.0, 69195019.0, 1871391091.0 },
{ 3672091415.0, 3528743235.0, 69195019.0 }
};
static double A2p76[3][3] = {
- { 1511326704.0, 3759209742.0, 1610795712.0 },
- { 4292754251.0, 1511326704.0, 3889917532.0 },
+ { 1511326704.0, 3759209742.0, 1610795712.0 },
+ { 4292754251.0, 1511326704.0, 3889917532.0 },
{ 3859662829.0, 4292754251.0, 3708466080.0 }
};
static double A1p127[3][3] = {
- { 2427906178.0, 3580155704.0, 949770784.0 },
+ { 2427906178.0, 3580155704.0, 949770784.0 },
{ 226153695.0, 1230515664.0, 3580155704.0 },
{ 1988835001.0, 986791581.0, 1230515664.0 }
};
RngStream RngStream_CopyStream (const RngStream src)
{
RngStream g;
-
+
if(src == NULL) {
printf ("RngStream_CopyStream: 'src' not initialized\n\n");
exit (EXIT_FAILURE);
return -1; /* FAILURE */
for (i = 0; i < 6; ++i)
g->Cg[i] = g->Bg[i] = g->Ig[i] = seed[i];
- return 0; /* SUCCESS */
+ return 0; /* SUCCESS */
}
/*-------------------------------------------------------------------------*/