+ # Serialization
+ binarize(as.matrix(c(as.double(Re(ts_cwt)),as.double(Im(ts_cwt)))), cwt_file, 1,
+ ",", nbytes, endian)
+ }
+
+ if (parll)
+ {
+ cl = parallel::makeCluster(ncores_clust)
+ synchrones_desc <- bigmemory::describe(synchrones)
+ Xwer_dist_desc <- bigmemory::describe(Xwer_dist)
+ parallel::clusterExport(cl, varlist=c("parll","synchrones_desc","Xwer_dist_desc",
+ "noctave","nvoice","verbose","getCWT"), envir=environment())
+ }
+
+ if (verbose)
+ cat(paste("--- Precompute and serialize synchrones CWT\n", sep=""))
+
+ ignored <-
+ if (parll)
+ parallel::parLapply(cl, 1:n, computeSaveCWT)
+ else
+ lapply(1:n, computeSaveCWT)
+
+ # Function to retrieve a synchrone CWT from (binary) file
+ getSynchroneCWT = function(index, L)
+ {
+ flat_cwt <- getDataInFile(index, cwt_file, nbytes, endian)
+ cwt_length = length(flat_cwt) / 2
+ re_part = as.matrix(flat_cwt[1:cwt_length], nrow=L)
+ im_part = as.matrix(flat_cwt[(cwt_length+1):(2*cwt_length)], nrow=L)
+ re_part + 1i * im_part
+ }
+
+ # Compute distance between columns i and j in synchrones
+ computeDistanceIJ = function(pair)
+ {
+ if (parll)