Rng.h

Go to the documentation of this file.

//                                                         //
// Copyright (c) 2003-2014 by The University of Queensland //
// Centre for Geoscience Computing                         //
// http://earth.uq.edu.au/centre-geoscience-computing      //
//                                                         //
// Primary Business: Brisbane, Queensland, Australia       //
// Licensed under the Open Software License version 3.0    //
// http://www.opensource.org/licenses/osl-3.0.php          //
//                                                         //


#ifndef ESYS_LSMRNG_H
#define ESYS_LSMRNG_H

#include <boost/random.hpp>

namespace esys
{
  namespace lsm
  {
    template <typename TmplRng=boost::mt19937>
    class UniformRng
    {
    public:
      typedef TmplRng Rng;
      typedef boost::uniform_real<> UniformReal;
      typedef boost::variate_generator<Rng &, UniformReal> Generator;

      UniformRng(double minRn, double maxRn);

      double operator()();

      void seed();

      template <typename Tmpl>
      void seed(Tmpl &s);

      template <typename TmplIt>
      void seed(TmplIt begin, TmplIt end);

    private:
      Rng         m_rng;
      UniformReal m_uniform;
      Generator   m_generator;
    };

    typedef UniformRng<> DefaultUniformRng;

    namespace rng
    {
      extern DefaultUniformRng s_zeroOneUniform;
    }
  }
}

#include "Foundation/Rng.hpp"

#endif