X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=src%2Fsources%2FEMGrank.c;h=0791892429ac28a2ee4e854ce76cafad9c0b73f0;hb=09ab3c164abb566764e86a175b5973241e708fd6;hp=2e85ce6fd019f9c836cc26d82210b1fc7d3a1556;hpb=552b00e200e8a990c1247989d29e98d4ae8679f3;p=valse.git

diff --git a/src/sources/EMGrank.c b/src/sources/EMGrank.c
index 2e85ce6..0791892 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)
@@ -37,7 +38,7 @@ static double* pinv(const double* matrix, int dim)
 }
 
 // TODO: comment EMGrank purpose
-void EMGrank(
+void EMGrank_core(
 	// IN parameters
 	const double* Pi, // parametre de proportion
 	const double* Rho, // parametre initial de variance renormalisé
@@ -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))
 	{