From ea45b0fd14b255c0c5a4515b2f27c68002a567c8 Mon Sep 17 00:00:00 2001
From: jbrv <jean-baptiste.herve@esial.net>
Date: Mon, 11 Jun 2012 17:07:52 +0200
Subject: [PATCH 1/1] Trace manager : add a function to generate values from
 probabilist distributions

---
 src/surf/trace_mgr.c         | 29 ++++++++++++++++++++++++++---
 src/surf/trace_mgr_private.h |  1 +
 2 files changed, 27 insertions(+), 3 deletions(-)

diff --git a/src/surf/trace_mgr.c b/src/surf/trace_mgr.c
index a73e539acc..16dce1e746 100644
--- a/src/surf/trace_mgr.c
+++ b/src/surf/trace_mgr.c
@@ -11,6 +11,7 @@
 #include "trace_mgr_private.h"
 #include "surf_private.h"
 #include "xbt/RngStream.h"
+#include <math.h>
 
 XBT_LOG_NEW_DEFAULT_SUBCATEGORY(surf_trace, surf, "Surf trace management");
 
@@ -61,7 +62,7 @@ probabilist_event_generator_t tmgr_event_generator_new_uniform(RngStream rng_str
   event_generator->s_uniform_parameters.beta = beta;
   event_generator->rng_stream = rng_stream;
 
-  //FIXME Generate a new event date
+  tmgr_event_generator_next_value(event_generator);
   
   return event_generator;
 }
@@ -76,7 +77,7 @@ probabilist_event_generator_t tmgr_event_generator_new_exponential(RngStream rng
   event_generator->s_exponential_parameters.lambda = lambda;
   event_generator->rng_stream = rng_stream;
 
-  //FIXME Generate a new event date
+  tmgr_event_generator_next_value(event_generator);
   
   return event_generator;
 }
@@ -93,11 +94,33 @@ probabilist_event_generator_t tmgr_event_generator_new_weibull(RngStream rng_str
   event_generator->s_weibull_parameters.k = k;
   event_generator->rng_stream = rng_stream;
 
-  // FIXME Generate a new event date
+  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;
+}
+
 tmgr_trace_t tmgr_trace_new_from_string(const char *id, const char *input,
                                         double periodicity)
 {
diff --git a/src/surf/trace_mgr_private.h b/src/surf/trace_mgr_private.h
index 9498be8ac5..e9a97572b2 100644
--- a/src/surf/trace_mgr_private.h
+++ b/src/surf/trace_mgr_private.h
@@ -70,5 +70,6 @@ typedef struct tmgr_history {
   xbt_heap_t heap;
 } s_tmgr_history_t;
 
+double tmgr_event_generator_next_value(probabilist_event_generator_t generator);
 
 #endif                          /* _SURF_TMGR_PRIVATE_H */
-- 
2.20.1