C++ (Cpp) RNG::rtgamma_rate Examples

Example #1

File: BridgeWrapper.cpp Project: jwindle/BayesBridge

//------------------------------------------------------------------------------
void rrtgamma_rate(double *x, double *scale, double *rate, double *right_t, int *num)
{
  RNG r;

  #ifdef USE_R
  GetRNGstate();
  #endif

  for(int i=0; i < *num; ++i){
    #ifdef USE_R
    if (i%SAMPCHECK==0) R_CheckUserInterrupt();
    #endif

    x[i] = r.rtgamma_rate(scale[i], rate[i], right_t[i]);
  }

  #ifdef USE_R
  PutRNGstate();
  #endif
}