Hello everyone,
I am a new user of Rcpp, I use it to rewrite a whole R program.
at the end of the program, my "final function is the following"
#include <RcppArmadillo.h>
#include <Rcpp.h>
#define _USE_MATH_DEFINES
#include <math.h>
using namespace Rcpp;
// [[Rcpp::depends("RcppArmadillo")]]
List main_function_C(arma::vec X0,arma::vec XF,double t0,
double tF, arma::mat Gamma,
NumericVector piks,arma::mat muks, List Cks,
int max_iter = 500){
List bounds = bounds_fun_C(piks,Cks,Gamma);
double sup_bound = 0.5 *(as<double>(bounds["alpha_bar"])
+ as<double>(bounds["delta_max"])
- as<double>(bounds["delta_min"]));
bool accept = false;
int kappa = 0;
arma::mat omegas(2,2);//arbitrary size, will change in loop
NumericVector Psi(5);//arbitrary size, will change in loop
for(int i = 0; i < max_iter; ++i){
kappa = rpois(1,(tF-t0)*sup_bound)[0];
if(kappa > 0){
Psi = stl_sort(runif(kappa,t0,tF));
NumericVector Upsilon = runif(kappa, 0, sup_bound);
omegas = rubb_it_C(X0, XF, t0, tF, Psi);
NumericVector phi_omega(kappa);
for(int i = 1; i < kappa+2; i++){
phi_omega[i-1] = phi_C(arma::trans(omegas.row(i)),
piks,muks,Cks,Gamma);
}
accept = all_C((phi_omega < Upsilon));
}//end if kappa >0
else{
accept = true;
}
if(accept){
break;
}
}//end for
arma::mat res(kappa+2,3);//result matrix,
if(kappa > 0){
res.submat(0,0,kappa+1,1) = omegas;
res(0,2) = t0;
for(int i =1; i < kappa+1;++i){
res(i,2) = Psi[i-1];
//Psi is of size kappa the index then goes till kappa-1
}
res(kappa+1,2) = tF;
}
else{
res(0,arma::span(0,1)) = trans(X0);
res(1,arma::span(0,1)) = trans(XF);
res(0,2) = t0;
res(kappa+1,2) = tF;
}
return List::create(Named("skel") = res,
Named("accepted") = accept);
}
This function calls other funtion that I have written, all of them
are attached on the cpp_funcs.cpp file.
The sourceCpp performs well and I can, from R obtain the following result
sourceCpp("cpp_funcs.cpp")
X0 <- c(0.5,0.5)
XF <- c(1,1)
t0 <- 0
tF <- 1
Gam <- diag(0.1,2)
ncomp=2
pis <- c(0.5,0.5)
mu <- matrix(c(-1,1,
1,1),ncol=2,nrow=ncomp)
cov <- list(diag(0.5,2),diag(1,2))
set.seed(123)
test <- main_function_C(X0,XF,t0,tF,Gam,pis,mu,cov)
test
#$skel
[,1] [,2] [,3]
#[1,] 0.5000000 0.500000 0.0000000
#[2,] -0.1261731 1.313880 0.4089769
#[3,] 0.5564577 1.069211 0.7883051
#[4,] 1.0000000 1.000000 1.0000000
#$accepted
#[1] TRUE
PROBLEM:
When I run the exact same code as above a large number of times, as this
for(i in 1:1000){
## same code as above
}
The r session aborts all the time, giving no error message but "R
encountered a fatal error".
I do not encounter this problem with the intermediate functions in
the attached file (and, of course, I'm not asking for a debugging
of those), I only have it with this one
Is this a problem with my code? with how the loop is written?The List object?
Is this a problem of the memory?
It does it either on Rstudio or Rgui
I work on windows 7, with R version 3.0.2
Rcpp 0.11.2 and RcppArmadillo 0.4.450.1.0
I also attach a R debugging file giving parameters for all
intermediate functions.
Sorry for the length of the function, I'm pretty sure the problem
comes from the main function, but it sadly depends on others (which,
as far as I know, don't pose problems)
Thanks in advance for any help,
Pierre Gloaguen
