RandomGeneration.cpp

Go to the documentation of this file.

// //////////////////////////////////////////////////////////////////////
// Import section
// //////////////////////////////////////////////////////////////////////
// STL
#include <cassert>
#include <sstream>
// Boost
#include <boost/version.hpp>
#if BOOST_VERSION_MACRO >= 103500
#include <boost/math/distributions/normal.hpp>
#endif // BOOST_VERSION_MACRO >= 103500
// StdAir
#include <stdair/basic/RandomGeneration.hpp>
 
namespace stdair {
  
  // //////////////////////////////////////////////////////////////////////
  RandomGeneration::RandomGeneration() : _generator (1) {
  }
 
  // //////////////////////////////////////////////////////////////////////
  RandomGeneration::RandomGeneration (const RandomSeed_T& iSeed)
    : _generator (iSeed) {
  }
 
  // //////////////////////////////////////////////////////////////////////
  RandomGeneration::RandomGeneration (const RandomGeneration& iRandomGeneration)
    : _generator (iRandomGeneration._generator) {
  }
 
  // //////////////////////////////////////////////////////////////////////
  RandomGeneration::~RandomGeneration() {
  }
 
  // //////////////////////////////////////////////////////////////////////
  void RandomGeneration::init (const RandomSeed_T& iSeed) {
    _generator.seed (iSeed);
  }
 
  // //////////////////////////////////////////////////////////////////////
  const std::string RandomGeneration::describe() const {
    std::ostringstream oStr;
    oStr << _generator;
    return oStr.str();
  }
 
  // //////////////////////////////////////////////////////////////////////
  RealNumber_T RandomGeneration::generateUniform01() {
    UniformGenerator_T lGenerator (_generator, boost::uniform_real<>(0, 1));
    return lGenerator();
  }
  
  // //////////////////////////////////////////////////////////////////////
  RealNumber_T RandomGeneration::generateUniform(const RealNumber_T& iMinValue,
                                                 const RealNumber_T& iMaxValue) {
    const Probability_T lVariateUnif01 = generateUniform01();
    const RealNumber_T lVariateUnif =
      iMinValue + lVariateUnif01 * (iMaxValue - iMinValue);
    return lVariateUnif;
  }
 
  // //////////////////////////////////////////////////////////////////////
  RealNumber_T RandomGeneration::generateNormal (const RealNumber_T& mu,
                                                 const RealNumber_T& sigma) {
    
#if BOOST_VERSION_MACRO >= 103500
    const Probability_T lVariateUnif = generateUniform01();
    const boost::math::normal lNormal (mu, sigma);
    const RealNumber_T lRealNumberOfRequestsToBeGenerated =
      boost::math::quantile (lNormal, lVariateUnif);
#else // BOOST_VERSION_MACRO >= 103500
    // TODO: rely on GSL when Boost version smaller than 1.35
    const RealNumber_T lRealNumberOfRequestsToBeGenerated = 0.0;
#endif // BOOST_VERSION_MACRO >= 103500
 
    return lRealNumberOfRequestsToBeGenerated;
        
  }
 
  // //////////////////////////////////////////////////////////////////////
  RealNumber_T RandomGeneration::generateExponential (const RealNumber_T& lambda) {
    ExponentialDistribution_T lExponentialDistribution (lambda);
    
    ExponentialGenerator_T lExponentialDistributionGenerator (_generator,
                                                              lExponentialDistribution); 
 
    // Generate a random variate, expressed in (fractional) day
    const RealNumber_T lExponentialVariateInDays =
      lExponentialDistributionGenerator();
 
    return lExponentialVariateInDays;
  }
 
}