Displaying 2 results from an estimated 2 matches for "gsl_rng_default".
2009 Dec 22
0
Rcpp 0.7.0: R/C++ integration now even easier
...snippet:
## use Rcpp to pass down a parameter for the seed, and a vector size
gslrng <- '
int seed = RcppSexp(s).asInt();
int len = RcppSexp(n).asInt();
gsl_rng *r;
gsl_rng_env_setup();
std::vector<double> v(len);
r = gsl_rng_alloc (gsl_rng_default);
gsl_rng_set (r, (unsigned long) seed);
for (int i=0; i<len; i++) {
v[i] = gsl_rng_get (r);
}
gsl_rng_free(r);
return RcppSexp(v).asSexp();
'
## turn into a function that R can call
## compileargs redundant on Debian/Ubuntu as...
2009 Dec 22
0
Rcpp 0.7.0: R/C++ integration now even easier
...snippet:
## use Rcpp to pass down a parameter for the seed, and a vector size
gslrng <- '
int seed = RcppSexp(s).asInt();
int len = RcppSexp(n).asInt();
gsl_rng *r;
gsl_rng_env_setup();
std::vector<double> v(len);
r = gsl_rng_alloc (gsl_rng_default);
gsl_rng_set (r, (unsigned long) seed);
for (int i=0; i<len; i++) {
v[i] = gsl_rng_get (r);
}
gsl_rng_free(r);
return RcppSexp(v).asSexp();
'
## turn into a function that R can call
## compileargs redundant on Debian/Ubuntu as...