+ if (parll)
+ {
+ require("bigmemory", quietly=TRUE)
+ require("epclust", quietly=TRUE)
+ synchrones <- bigmemory::attach.big.matrix(synchrones_desc)
+ Xwer_dist <- bigmemory::attach.big.matrix(Xwer_dist_desc)
+ }
+
+ computeCWT = function(index)
+ {
+ ts <- scale(ts(synchrones[index,]), center=TRUE, scale=scaled)
+ totts.cwt = Rwave::cwt(ts, totnoct, nvoice, w0, plot=FALSE)
+ ts.cwt = totts.cwt[,s0log:(s0log+noctave*nvoice)]
+ #Normalization
+ sqs <- sqrt(2^(0:(noctave*nvoice)/nvoice)*s0)
+ sqres <- sweep(ts.cwt,2,sqs,'*')
+ sqres / max(Mod(sqres))
+ }
+