[Rcpp-devel] Accessing list of list members to permute them

Le 30 mars 2015 ? 11:05, ogami musashi <uragami at hotmail.com> a ?crit :

Hello

I'm a Rcpp newbie, so sorry if the question is trivial.

I have an R object which is a list of 1000 elements. Each elements is a result from a discrete wavelet transform (from package wmtsa) with 5 elements. One of them have a 50000 rows and 16 colums.

The first level (1000) are in fact random signals of 50 000 days. On each of them i ran a DWT that separated each signal into 16 components (of decreasing frequency) for each day thus resulting in a 50000 days * 16 components matrix for each of the 1000 signals.

I would like to sort it like that: for each component (a list of 16 elements) having for each day (50000 rows) the 1000 signals (1000 columns)


here are the objects in R:

dwtlist: 1000 elements * 5 elements* (50000*16) matrix
sortdwtlist:16 elements*(50000*1000) matrix

the code in R:

for(i in 1:16){

   for(j in 1:50000){

       for (k in 1:1000){

           sortdwtlist[[i]][j,k]<-dwtlist[[k]][[1]][j,i]
}
}
}


This takes about 27 minutes to complete.

I thus tried to put in C++ using Rcpp (i'm using sourceCpp):

#include <Rcpp.h>
using namespace Rcpp;

// [[Rcpp::export]]

List rearrangelist(List x){
   int ni=15;
   int nj=49999;
   int nk=999;
   List output;

for (int i= 0 ; i<ni ; i++){
   for (int j=0 ; j<nj ; j++) {
       for (int k=0 ; j<nk ; k++){
           output[i](j,k)=x[k][0](j,i);
       }
   }
}
return output;
}


But it doesn't work. error stops at the line of permutation, so i'm sure there's something with the types inside the list.
Can someone help me?

Thank you!