- double* phiLambda = (double*)malloc(lengthA*m*k*sizeof(double));
- double* rhoLambda = (double*)malloc(m*m*k*sizeof(double));
- double* piLambda = (double*)malloc(k*sizeof(double));
- double* LLF = (double*)malloc((maxi+1)*sizeof(double));
- double* S = (double*)malloc(lengthA*m*k*sizeof(double));
- EMGLLF(phia,rhoInit,piInit,gamInit,mini,maxi,gamma,0.0,Xa,Y,tau,
+ Real* phiLambda = (Real*)malloc(lengthA*m*k*sizeof(Real));
+ Real* rhoLambda = (Real*)malloc(m*m*k*sizeof(Real));
+ Real* piLambda = (Real*)malloc(k*sizeof(Real));
+ Real* LLF = (Real*)malloc((maxi+1)*sizeof(Real));
+ Real* S = (Real*)malloc(lengthA*m*k*sizeof(Real));
+ EMGLLF_core(phia,rhoInit,piInit,gamInit,mini,maxi,gamma,0.0,Xa,Y,tau,