X-Git-Url: https://git.auder.net/?p=epclust.git;a=blobdiff_plain;f=epclust%2Fsrc%2FcomputeMedoidsIndices.cpp;h=f24758411ce9d30b0670c3fb5f26811b41cb5708;hp=693418118f9f3d38d1007c3923ea534c44703ec9;hb=0fe757f750f51e580d2c5a7b7f7df87cc405d12d;hpb=0486fbadb122cb4d78c5d9f248c29800a59eb24e

diff --git a/epclust/src/computeMedoidsIndices.cpp b/epclust/src/computeMedoidsIndices.cpp
index 6934181..f247584 100644
--- a/epclust/src/computeMedoidsIndices.cpp
+++ b/epclust/src/computeMedoidsIndices.cpp
@@ -21,9 +21,9 @@ IntegerVector computeMedoidsIndices(SEXP pMedoids, NumericMatrix ref_series)
 {
 	XPtr<BigMatrix> pMed(pMedoids);
 	MatrixAccessor<double> medoids = MatrixAccessor<double>(*pMed);
-	int nb_series = ref_series.nrow(),
-		K = pMed->nrow(),
-		L = pMed->ncol();
+	int nb_series = ref_series.ncol(),
+		K = pMed->ncol(),
+		L = pMed->nrow();
 	IntegerVector mi(nb_series);
 
 	for (int i=0; i<nb_series ; i++)
@@ -36,7 +36,7 @@ IntegerVector computeMedoidsIndices(SEXP pMedoids, NumericMatrix ref_series)
 			double dist_ij = 0.;
 			for (int k=0; k<L; k++)
 			{
-				double delta = ref_series(i,k) - *(medoids[k]+j);
+				double delta = ref_series(k,i) - *(medoids[j]+k);
 				dist_ij += delta * delta;
 			}
 			if (dist_ij < best_dist)