Hi, guys. I posted this by accident at rcpp-dev, although it meant to
be only to r-dev, so don't flame me here please, rcpp guys will
do it there, i am sure :).
I have some pretty large arrays in R and i wanted to do some time
consuming modifications of these arrays in C++ without actually copying
them, just by passing pointers to them. Since i don't know internal data
structures of R, i am not sure it's possible, but i thought it was. Here is
some toy code that i thought should work, but doesn't. Maybe someone could
point out the error i am making
i have the following in the passptr.cpp to multiply array elements by 2
===============================
extern "C"{
void modify(double *mem, int *nr, int *nc){
for(int i=0; i< (*nr)*(*nc); i++)
mem[i]=2*mem[i];
}
}
----------------------------------------------
I compile it into a shared library using
R CMD SHLIB passptr.cpp
load and run from R as follows
--------------------------------
>dyn.load("/home/az05625/testarma/passptr.so")
>m<-matrix(1:10,nr=2)
>.C("modify", as.double(m), as.integer(2), as.integer(5), DUP=FALSE)
>From reading docs i thought that DUP=FALSE would ensure that R matrix is
not copied and is multiplied by 2 in place. However, it's not the case,
matrix m is the same after calling .C("modify"...)
as it was before. Am i calling incorrectly, or is it just impossible to
modify R matrix in place from C++? Would greatly appreciate any pointers.
[[alternative HTML version deleted]]
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
