[Rd] Crash in R using embedded.

From: Dave Pugmire <dpugmire_at_gmail.com>
Date: Wed, 01 Feb 2012 10:24:23 -0500

I'm new to R, and am trying to embed R into another application. I'm calling gev.fit() from the ismev package, and it is crashing somewhere inside it.
gdb is not catching it, and valgrind is not showing any memory corruption issues.
I suspect it's memory corruption, because it doesn't crash in exactly the same spot each time.
I'm running R 2.12.2 on a 64 bit linux (Ubuntu 10.04)

Also, before I make the gev.fit() call, I save out the R state using "save.image(....)", then load the image file into a standalone R on the console, and call gev.fit(), and it works fine. So, it's something to do with my app, and the R.

I'm looking for hints/ideas/tricks/etc for debugging this. I'm starting to put print statements into the R code, and into the C code in R (src/main/optim.c), but it hasn't shown me the source of the problem.


Here is the code I'm using to initialize:

putenv("R_HOME=/apps/R/install/lib64/R"); const char *R_argv[]= {"RInterfaceTest", "--gui=none", "--no-save", "--no-readline", "--silent"};
Rf_initialize_R(sizeof(R_argv)/sizeof(R_argv[0]), const_cast<char


ParseStatus status;
SEXP cmdSexp, cmdexpr = R_NilValue;
int error;
string rcommand = "f<-file(paste(tempdir(), \"/Routput.txt\", sep = \"\"), open=\"wt+\")\nsink(f)\n";

PROTECT(cmdSexp = allocVector(STRSXP, 1)); SET_STRING_ELT(cmdSexp, 0, mkChar(rcommand.c_str()));

cmdexpr = PROTECT(R_ParseVector(cmdSexp, -1, &status, R_NilValue)); for(int i = 0; i < length(cmdexpr); i++) {

     R_tryEval(VECTOR_ELT(cmdexpr, i),NULL,&error);  }
 UNPROTECT(2);         [[alternative HTML version deleted]]

R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Wed 01 Feb 2012 - 15:38:31 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 01 Feb 2012 - 19:20:13 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive