[Rcpp-devel] Rcpp matrix generation crashes my R session

[Dale Smith]

I would say you should create a simpler example that still crashes. Comment out everything except the first for loop; see if it crashes then. If not, uncomment the next one, etc. Repeat this until you find where it crashes.

Do you have a one-off problem with a loop index? I think you do, but will let you look into this yourself.

On Jun 22, 2013, at 12:10 PM, Marc Jekel <mjekel at uni-bonn.de> wrote:

> Dear Rcpp users,
> 
> I am a beginner, I was inspired by the new book by Eddelbuettel and just started to use Rccp. So far, everything works fine except for one function. The strange thing about it is that I runs perfectly when I call it several times but crashes my R session at some point (usually after the third call). I checked it on two Win 8 PCS to make sure it is nothing related to my PC.
> 
> Here is the code:
> 
> ##########
> 
> #include <Rcpp.h>
> using namespace Rcpp;
> 
> // [[Rcpp::export]]
> NumericMatrix transPCSpairwiseComp_C(
> NumericVector validities,
> NumericMatrix pattern,
> int numbCues
> ){
> 
> NumericMatrix weightsNet(numbCues+3,numbCues+3);
> 
> for(int loop = 0; loop < numbCues; loop++){
> 
>  weightsNet(0,1+loop) = validities[loop];
> }
> 
> for(int loopCol = 0; loopCol < numbCues; loopCol++){
> for(int loop = 0; loop < numbCues; loop++){
> weightsNet(1+loopCol, numbCues + 1 + loop)=
> pattern(loopCol,loop);
> 
> }}
> 
> weightsNet(numbCues+1, numbCues+1) = 0;
> weightsNet(numbCues+1, numbCues+2) = -.2;
> 
> for(int LoopCol = 0; LoopCol < weightsNet.nrow(); LoopCol++){
> for(int LoopRow = 0; LoopRow < weightsNet.nrow(); LoopRow++){
> weightsNet(LoopCol,LoopRow) = weightsNet(LoopCol,LoopRow) + weightsNet(LoopRow,LoopCol);
> }}
> 
> 
> return weightsNet;
> }
> 
> #########
> 
> So it is very basic code, a matrix is created based on three inputs (the vector "validities", the matrix "pattern", and the integer "numbCues"). My impression is that the very first command (creating the matrix weightsNet) is problematic because R also crashes when I delete all other code (although usually not after the third time I call it).
> 
> Thanks for any hint!
> 
> Marc
> 
> 
> 
> _______________________________________________
> Rcpp-devel mailing list
> Rcpp-devel at lists.r-forge.r-project.org
> https://lists.r-forge.r-project.org/cgi-bin/mailman/listinfo/rcpp-devel