- double* phi = (double*)malloc(p*m*k*sizeof(double));
- double* rho = (double*)malloc(m*m*k*sizeof(double));
- 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,
+ Real* phi = (Real*)malloc(p*m*k*sizeof(Real));
+ Real* rho = (Real*)malloc(m*m*k*sizeof(Real));
+ Real* pi = (Real*)malloc(k*sizeof(Real));
+ Real* LLF = (Real*)malloc(maxi*sizeof(Real));
+ Real* S = (Real*)malloc(p*m*k*sizeof(Real));
+ EMGLLF_core(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambda[lambdaIndex],X,Y,tau,