+RngStream tmgr_rng_stream_from_id(char* id)
+{
+ unsigned int id_hash;
+ RngStream rng_stream = NULL;
+
+ rng_stream = RngStream_CopyStream(common_rng_stream);
+ id_hash = xbt_dict_hash(id);
+ RngStream_AdvanceState(rng_stream, 0, id_hash);
+
+ return rng_stream;
+}
+
+probabilist_event_generator_t tmgr_event_generator_new_uniform(RngStream rng_stream,
+ double alpha,
+ double beta)
+{
+ probabilist_event_generator_t event_generator = NULL;
+
+ 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->rng_stream = rng_stream;
+
+ tmgr_event_generator_next_value(event_generator);
+
+ return event_generator;
+}
+
+probabilist_event_generator_t tmgr_event_generator_new_exponential(RngStream rng_stream,
+ double lambda)
+{
+ probabilist_event_generator_t event_generator = NULL;
+
+ 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->rng_stream = rng_stream;
+
+ tmgr_event_generator_next_value(event_generator);
+
+ return event_generator;
+}
+
+probabilist_event_generator_t tmgr_event_generator_new_weibull(RngStream rng_stream,
+ double lambda,
+ double k)
+{
+ probabilist_event_generator_t event_generator = NULL;
+
+ 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->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;
+ break;
+ case e_generator_exponential:
+ generator->next_value = -log(RngStream_RandU01(generator->rng_stream))
+ / generator->s_exponential_parameters.lambda;
+ break;
+ case e_generator_weibull:
+ generator->next_value = - generator->s_weibull_parameters.lambda
+ * pow( log(RngStream_RandU01(generator->rng_stream)),
+ 1.0 / generator->s_weibull_parameters.k );
+ }
+
+ return generator->next_value;
+}
+