#include <stdlib.h>
#include <omp.h>
+#include <math.h>
#include "EMGLLF.h"
#include "utils.h"
// Main job on raw inputs (after transformation from mxArray)
-void selectiontotale(
+void selectiontotale_core(
// IN parameters
const double* phiInit, // parametre initial de moyenne renormalisé
const double* rhoInit, // parametre initial de variance renormalisé
double* pi = (double*)malloc(k*sizeof(double));
double* LLF = (double*)malloc(maxi*sizeof(double));
double* S = (double*)malloc(p*m*k*sizeof(double));
- EMGLLF(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambda[lambdaIndex],X,Y,tau,
+ EMGLLF_core(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambda[lambdaIndex],X,Y,tau,
phi,rho,pi,LLF,S,
n,p,m,k);
free(LLF);