rng.h

Go to the documentation of this file.

// include guard
#ifndef TF_RNG_H_VERSION

#define TF_RNG_H_VERSION \
  "TF_RNG_H   V1.1"

#include <tfxx/error.h>
#include <ostream>
#include <string>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>

namespace tfxx {

  namespace numeric {

  /* ---------------------------------------------------------------------- */

    class RNGgaussian {
      public:
        explicit RNGgaussian(const double& std=1.,
                             const double& mean=0.,
                             const char* type="default");
        ~RNGgaussian();
        double operator()() const;
        double value() const;
        void set(const unsigned long int& seed);
        void set();
        std::string type() const {return std::string(MR->type->name); }
        unsigned long int seed() const { return Mseed; }
        double std() const { return Mstd; }
        double mean() const { return Mmean; }
        static void rng_list_types(std::ostream& os);
        static const char* comment_gsl_rng_ranlux;
      private:
        double Mstd, Mmean;
        gsl_rng* MR;
        unsigned long int Mseed;
      public:
        class UTException: public tfxx::error::Exception
      {
        public:
          typedef tfxx::error::Exception TBase;
          UTException(const char* type, 
                      const char* file,
                      const int& line,
                      const char* condition);
          virtual ~UTException() { }
        private:
          const char* Mtype;
      }; // class UTException
    }; // class RNGgaussian
  } // namespace numeric

} // namespace tfxx

#endif // TF_RNG_H_VERSION (includeguard)

/* ----- END OF rng.h ----- */