X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=src%2Fsources%2FEMGrank.c;fp=src%2Fsources%2FEMGrank.c;h=5ee44e38f0ce38eef711ed2ef979fd88ec63a46c;hb=8e92c49c15bdacebf46190e7c8279bd227873928;hp=2e85ce6fd019f9c836cc26d82210b1fc7d3a1556;hpb=3ec579a0955aca0591a9e5c4d90c50b87f4f4609;p=valse.git

diff --git a/src/sources/EMGrank.c b/src/sources/EMGrank.c
index 2e85ce6..5ee44e3 100644
--- a/src/sources/EMGrank.c
+++ b/src/sources/EMGrank.c
@@ -1,5 +1,6 @@
-#include "EMGrank.h"
+#include <stdlib.h>
 #include <gsl/gsl_linalg.h>
+#include "utils.h"
 
 // Compute pseudo-inverse of a square matrix
 static double* pinv(const double* matrix, int dim)
@@ -82,8 +83,7 @@ void EMGrank(
 	int* Z = (int*)calloc(n, sizeof(int));
 
 	//Initialize phi to zero, because some M loops might exit before phi affectation
-	for (int i=0; i<p*m*k; i++)
-		phi[i] = 0.0;
+	zeroArray(phi, p*m*k);
 
 	while (ite<mini || (ite<maxi && sumDeltaPhi>tau))
 	{