work/mixpp/bdm/itpp_ext.h

//
// C++ Interface: itpp_ext
//
// Description: 
//
//
// Author: smidl <smidl@utia.cas.cz>, (C) 2008
//
// Copyright: See COPYING file that comes with this distribution
//
//
#include <itpp/itbase.h>

using std::cout;
using std::endl;

namespace itpp {
        Array<int> to_Arr(const ivec &indices);
        
  class Gamma_RNG {
  public:
    Gamma_RNG(double a=1.0, double b=1.0);
    void setup(double a0, double b0) { alpha=a0; beta=b0;}
    double get_setup() const;
    double operator()() { return sample(); }
    vec operator()(int n);
    mat operator()(int h, int w);
  protected:
  private:
    double sample();
    double alpha;
    double beta;
    Random_Generator RNG;
    Normal_RNG NRNG;
    inline double exp_rand(){return -std::log(RNG.random_01());}
    inline double unif_rand(){return RNG.random_01();}
    inline double norm_rand(){return NRNG.sample();}

  };
}

---