setting the seed in standalone code using Rlib

David Faden, dfaden at iastate.edu
AIM: pitulx
--

#include <assert.h>
#include <R.h>
#include <Rinternals.h>
#include <R_ext/Parse.h>
#include <Rmath.h>
#include <stdio.h>

void setSeed(unsigned int seed)
{
#ifndef MATHLIB_STANDALONE

  /**
   * Following the example in Writing R Extensions 5.10
   */

  char cmd[256];
  SEXP cmdSexp;
  SEXP cmdExp;
  ParseStatus status;

  sprintf(cmd, "set.seed(%u)", seed);

  PROTECT(cmdSexp = allocVector(STRSXP, 1));
  SET_STRING_ELT(cmdSexp, 0, mkChar(cmd));
  cmdExp = PROTECT(R_ParseVector(cmdSexp, -1, &status, R_NilValue));

  assert(status == PARSE_OK);

  eval(VECTOR_ELT(cmdExp, 0), R_GlobalEnv);

  UNPROTECT(2);

#else

  set_seed(seed, seed);

#endif

}

/* Test */

extern int Rf_initEmbeddedR(int argc, const char *argv[]);

int main(int numArgs, char** args)
{
  const char *argv[]= {"setseed", "--gui=none", "--silent"};
  const int argc = 3;

  Rf_initEmbeddedR(argc, argv);

  setSeed((unsigned int)42);
  printf("%f\n", rnorm(0.0, 1.0));
  return 0;
}