+#TODO: use dbs(),
+ #https://www.r-bloggers.com/debugging-parallel-code-with-dbs/
+ #http://gforge.se/2015/02/how-to-go-parallel-in-r-basics-tips/
+
+#OK ::
+#write(length(indices), file="TOTO")
+#write( computeMedoidsIndices(medoids@address, getRefSeries(indices[1:600])), file="TOTO")
+#stop()
+
+# write(indices, file="TOTO", ncolumns=10, append=TRUE)
+#write("medoids", file = "TOTO", ncolumns=1, append=TRUE)
+#write(medoids[1,1:3], file = "TOTO", ncolumns=1, append=TRUE)
+#write("synchrones", file = "TOTO", ncolumns=1, append=TRUE)
+#write(synchrones[1,1:3], file = "TOTO", ncolumns=1, append=TRUE)
+
+#NOT OK :: (should just be "ref_series") ...or yes ? race problems mutex then ? ?!
+ #get medoids indices for this chunk of series
+ mi = computeMedoidsIndices(medoids@address, getRefSeries(indices[1:600])) #ref_series)
+write("MI ::::", file = "TOTO", ncolumns=1, append=TRUE)
+write(mi[1:3], file = "TOTO", ncolumns=1, append=TRUE)
+
+# #R-equivalent, requiring a matrix (thus potentially breaking "fit-in-memory" hope)
+# mat_meds = medoids[,]
+# mi = rep(NA,nb_series)
+# for (i in 1:nb_series)
+# mi[i] <- which.min( rowSums( sweep(mat_meds, 2, ref_series[i,], '-')^2 ) )
+# rm(mat_meds); gc()