exponentialdistribution.h

Go to the documentation of this file.

#ifndef EXPONENTIALDISTRIBUTION_H

#define EXPONENTIALDISTRIBUTION_H

#include "probabilitydistribution.h"
#include <cmath>

class ExponentialDistribution : public ProbabilityDistribution
{
public:
        ExponentialDistribution(double a, GslRandomNumberGenerator *pRng) : ProbabilityDistribution(pRng)       { assert(a > 0); m_a = a; }

        double pickNumber() const;
        double getA() const                                                                                     { return m_a; }
private:
        double m_a;
};

inline double ExponentialDistribution::pickNumber() const
{
        double z = getRandomNumberGenerator()->pickRandomDouble();
        double y = -std::log(z);
        double x = y/m_a;

        return x;
}

#endif // EXPONENTIALDISTRIBUTION_H