rand.h

Go to the documentation of this file.

#ifndef _RAND_H
#define _RAND_H

#include <QObject>
#include <QDataStream>
#include <QTextStream>

#ifndef Q_MOC_RUN

#if __GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 6)
#pragma GCC diagnostic push
#pragma GCC system_header
#endif

#include <boost/random.hpp>

#if __GNUC__ > 4 || (__GNUC__ == 4 && __GNUC_MINOR__ >= 6)
#pragma GCC diagnostic pop
#endif

#endif

#include <ctime>

namespace LangmuirCore
{

class Random : public QObject
{
private:
    Q_OBJECT
    Q_DISABLE_COPY(Random)

public:
    Random(quint64 seed=0, QObject *parent=0);

    ~Random();

    quint64 seed();

    void seed(quint64 seed);

    double random();

    double range(const double low=0.0, const double high=1.0);

    double normal(const double mean, const double sigma);

    int integer(const int low=0, const int high=1);

    bool metropolis(double energyChange, double inversekT);

    bool metropolisWithCoupling(double energyChange, double inversekT, double coupling);

    bool chooseYes(double percent);

    bool chooseNo(double percent);

    friend QDataStream& operator<<(QDataStream& stream, Random& random);

    friend QDataStream& operator>>(QDataStream& stream, Random& random);

    friend QTextStream& operator<<(QTextStream& stream, Random& random);

    friend QTextStream& operator>>(QTextStream& stream, Random& random);

    friend std::ostream& operator<<(std::ostream& stream, Random& random);

    friend std::istream& operator>>(std::istream& stream, Random& random);

private:
    boost::mt19937 *twister;

    boost::variate_generator<boost::mt19937&, boost::uniform_01<double> >* generator01;

    quint64 m_seed;
};

}
#endif