Logo AND Algorithmique Numérique Distribuée

Public GIT Repository
Fix a mistake in the Weibull inverse distribution function
authorJean-Baptiste Hervé <jean-baptiste.herve@esial.net>
Mon, 9 Jul 2012 08:10:30 +0000 (10:10 +0200)
committerJean-Baptiste Hervé <jean-baptiste.herve@esial.net>
Mon, 9 Jul 2012 08:10:30 +0000 (10:10 +0200)
src/surf/trace_mgr.c

index c541762..aa70a21 100644 (file)
@@ -133,8 +133,8 @@ double tmgr_event_generator_next_value(probabilist_event_generator_t generator)
                               / 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)),
+      generator->next_value = generator->s_weibull_parameters.scale
+                              * pow( -log(RngStream_RandU01(generator->rng_stream)),
                                     1.0 / generator->s_weibull_parameters.shape );
   }