- //[phiLambda,rhoLambda,piLambda,~,~] = EMGLLF(...
- // phiInit(a,:,:),rhoInit,piInit,gamInit,mini,maxi,gamma,0,X(:,a),Y,tau);
- 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_core(phia,rhoInit,piInit,gamInit,mini,maxi,gamma,0.0,Xa,Y,tau,
+ //Call to EMGLLF
+ 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.,Xa,Y,tau,