getSeries = function(inds) getDataInFile(inds, series_file, nbytes, endian)
}
- # Serialize all computed wavelets contributions onto a file
+ # Serialize all computed wavelets contributions into a file
contribs_file = paste(bin_dir,"contribs",sep="") ; unlink(contribs_file)
index = 1
nb_curves = 0
inds, getContribs, K1, nb_series_per_chunk, ncores_clust, verbose, parll)
if (WER=="mix")
{
+
+
+
+
+#TODO: getSeries(indices_medoids) BAD ; il faudrait une big.matrix de medoids en entree
+ #OK en RAM il y en aura 1000 (donc 1000*K1*17519... OK)
+ #...mais du coup chaque process ne re-dupliquera pas medoids
+
+
medoids2 = computeClusters2(getSeries(indices_medoids),
K2, getSeries, nb_curves, nb_series_per_chunk, ncores_clust, verbose, parll)
binarize(medoids2, synchrones_file, nb_series_per_chunk, sep, nbytes, endian)
contribs_file, nb_series_per_chunk, nbytes, endian)
}
+
+
+#TODO: if ntasks==1, c'est deja terminé
+
# Run step2 on resulting indices or series (from file)
if (verbose)
cat("...Run final // stage 1 + stage 2\n")