X-Git-Url: https://git.auder.net/?p=epclust.git;a=blobdiff_plain;f=epclust%2FR%2Fclustering.R;h=9a55495bb6f650a3445c11b281ddef359dd0c4c2;hp=14e1f835a495254b1cde41b44b52b70fa063a2df;hb=c45fd66342e40c8b5387fc6f0059c4d3a9718340;hpb=af3ea94782db30c2fc1bd796d788dab5e1fb79ed

diff --git a/epclust/R/clustering.R b/epclust/R/clustering.R
index 14e1f83..9a55495 100644
--- a/epclust/R/clustering.R
+++ b/epclust/R/clustering.R
@@ -93,6 +93,8 @@ computeSynchrones = function(medoids, getRefSeries,
 {
 	computeSynchronesChunk = function(indices)
 	{
+		if (verbose)
+			cat(paste("--- Compute synchrones for ",length(indices)," lines\n", sep=""))
 		ref_series = getRefSeries(indices)
 		#get medoids indices for this chunk of series
 		for (i in seq_len(nrow(ref_series)))
@@ -118,17 +120,14 @@ computeSynchrones = function(medoids, getRefSeries,
 		m <- synchronicity::boost.mutex()
 
 	indices_workers = .spreadIndices(seq_len(nb_ref_curves), nb_series_per_chunk)
-	for (inds in indices_workers)
-	{
-		if (verbose)
-			cat(paste("--- Compute synchrones for ",length(inds)," lines\n", sep=""))
+	ignored <-
 		if (parll)
-			ignored <- parallel::mcparallel(computeSynchronesChunk(inds))
+		{
+			parallel::mclapply(indices_workers, computeSynchronesChunk,
+				mc.cores=ncores_clust, mc.allow.recursive=FALSE)
+		}
 		else
-			computeSynchronesChunk(inds)
-	}
-	if (parll)
-		parallel::mccollect()
+			lapply(indices_workers, computeSynchronesChunk)
 
 	mat_syncs = matrix(nrow=K, ncol=ncol(medoids))
 	vec_count = rep(NA, K)
@@ -234,18 +233,16 @@ computeWerDists = function(synchrones, ncores_clust=1,verbose=FALSE,parll=TRUE)
 	if (parll)
 		m <- synchronicity::boost.mutex()
 
-	for (i in 1:(n-1))
-	{
+	ignored <-
 		if (parll)
-			ignored <- parallel::mcparallel(computeDistancesLineI(i))
+		{
+			parallel::mclapply(seq_len(n-1), computeDistancesLineI,
+				mc.cores=ncores_clust, mc.allow.recursive=FALSE)
+		}
 		else
-			computeDistancesLineI(i)
-	}
+			lapply(seq_len(n-1), computeDistancesLineI)
 	Xwer_dist[n,n] = 0.
 
-	if (parll)
-		parallel::mccollect()
-
 	mat_dists = matrix(nrow=n, ncol=n)
 	#TODO: avoid this loop?
 	for (i in 1:n)