X-Git-Url: https://git.auder.net/?p=epclust.git;a=blobdiff_plain;f=epclust%2FR%2Fclustering.R;h=14e1f835a495254b1cde41b44b52b70fa063a2df;hp=74d009e45495220304a0087554e619d3c9efc4e4;hb=af3ea94782db30c2fc1bd796d788dab5e1fb79ed;hpb=492cd9e74a79cbcc0ecde55fa3071a44b7e463dc diff --git a/epclust/R/clustering.R b/epclust/R/clustering.R index 74d009e..14e1f83 100644 --- a/epclust/R/clustering.R +++ b/epclust/R/clustering.R @@ -121,10 +121,7 @@ computeSynchrones = function(medoids, getRefSeries, for (inds in indices_workers) { if (verbose) - { - cat(paste("--- Compute synchrones for indices range ", - min(inds)," -> ",max(inds),"\n", sep="")) - } + cat(paste("--- Compute synchrones for ",length(inds)," lines\n", sep="")) if (parll) ignored <- parallel::mcparallel(computeSynchronesChunk(inds)) else @@ -184,14 +181,16 @@ computeWerDists = function(synchrones, ncores_clust=1,verbose=FALSE,parll=TRUE) ts.cwt = totts.cwt[,s0log:(s0log+noctave*nvoice)] #Normalization sqs <- sqrt(2^(0:(noctave*nvoice)/nvoice)*s0) - sqres <- sweep(ts.cwt,MARGIN=2,sqs,'*') + sqres <- sweep(ts.cwt,2,sqs,'*') sqres / max(Mod(sqres)) } if (parll) { cl = parallel::makeCluster(ncores_clust) - parallel::clusterExport(cl, varlist=c("getContribs","K1","verbose"), envir=environment()) + parallel::clusterExport(cl, + varlist=c("synchrones","totnoct","nvoice","w0","s0log","noctave","s0","verbose"), + envir=environment()) } # (normalized) observations node with CWT