STAT_UniformDistribution.h

Go to the documentation of this file.

#ifndef STAT_UniformDistribution_H
#define STAT_UniformDistribution_H

#include "STAT_Distribution.h"
#include <boost/math/distributions/uniform.hpp>
#include <boost/random/uniform_01.hpp>
#include <boost/random/uniform_real.hpp>
#include <boost/random/uniform_int.hpp>
#include <boost/random/variate_generator.hpp>

DEFINE_SPTR(STAT_UniformDistribution);
class STAT_UniformDistribution : public STAT_Distribution {
    SP_OBJECT(STAT_UniformDistribution);
    // ATTRIBUTES
  
private:

  
  
    // ASSOCIATIONS
    boost::math::uniform_distribution<tReal> *mDistributionReal;
    boost::math::uniform_distribution<int> *mDistributionInt;
    
    boost::uniform_01<tReal> *mRandomDistribution01;
    boost::uniform_real<tReal> *mRandomDistributionReal;
    boost::uniform_int<int> *mRandomDistributionInt;
    
    boost::random::variate_generator<RNGType&,
                                     boost::uniform_01<tReal> > *mGenerator01;
    boost::random::variate_generator<RNGType&,
                                     boost::uniform_real<tReal> > *mGeneratorReal;
    boost::random::variate_generator<RNGType&,
                                     boost::uniform_int<int> > *mGeneratorInt;
    
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    STAT_UniformDistribution();
  
  
  
    // DESTRUCTORS 
  
  
    virtual ~STAT_UniformDistribution(void);
private:
    void freeAttributes();
  
  
  
public:
    //----------------
    // NEW constructors
    // ---------------

    inline static SP::STAT_UniformDistribution New(const tReal& lower,const tReal& upper) {
        SP::STAT_UniformDistribution p=New();
        p->setProperties(lower,upper);
        return p;
    }
    inline static SP::STAT_UniformDistribution New(const int& lower,const int& upper) {
        SP::STAT_UniformDistribution p=New();
        p->setProperties(lower,upper);
        return p;
    }
    inline static SP::STAT_UniformDistribution New() {
        SP::STAT_UniformDistribution p(new STAT_UniformDistribution(),
                                     STAT_UniformDistribution::Delete());
        p->setThis(p);
        return p;
    }
    virtual SP::STAT_Distribution NewInstance() const {
        return New();
    }
    
    
    virtual void copy(const STAT_Distribution& d) {
        tReal a,b;
        tReal m=d.mean();
        tReal v=d.variance();
        b=m+sqrt(3*v);
        a=m-sqrt(3*v);
        setProperties(a,b);
    }
    // --------------------
    // ATTRIBUTES setting
    // --------------------

    void setProperties(const tReal& lower,const tReal& upper);
    void setProperties(const int&  lower,const int& upper);
    
    // -----------------------
    // USUAL Functions
    // ------------------------

    
    virtual tReal pdf(const tReal& x) const;
    
    virtual tReal cdf(const tReal& x) const;

    virtual tReal cdfc(const tReal& x) const;

    virtual void range(tReal& xmin,tReal& xmax) const;
    
    virtual tReal mean() const;

    virtual tReal variance() const;
    virtual tReal standardDeviation() const;
    
    virtual tReal random() const;
    
    
    virtual void computeSample(const int& n,CORE_RealArray& sample) const;
      
    
      
protected:
    /* \brief return the growth rate of the distribution
     * @param adjust: is true rescale the parameters of the distribution suach as its growth rate is rho
     * 
     */
    virtual tReal computeGrowthRate(const tReal& v,
                                    const tBoolean& adjust,
                                    const tReal& rho);
  
    
    
    
  
   
};

#endif