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

[ogami musashi]

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!