int *MakeRandomIndices(int nx, int nmin, int nmax) { int nrand = nmax-nmin; int *random = (int *) malloc(nx*sizeof(int)); for(int i=0;i<nx;i++) random[i] = rng_integer(nrand) + nmin; return random; }
size_t Runif::runif(size_t min_val, size_t max_val) const { return min_val + (rng_integer(state) % (max_val - min_val)); }
static double rng_double(size_t &state) { return rng_integer(state)/DOUBLE_DENOMINATOR; }