- 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,
+ float* phiLambda = (float*)malloc(lengthA*m*k*sizeof(float));
+ float* rhoLambda = (float*)malloc(m*m*k*sizeof(float));
+ float* piLambda = (float*)malloc(k*sizeof(float));
+ float* LLF = (float*)malloc((maxi+1)*sizeof(float));
+ float* S = (float*)malloc(lengthA*m*k*sizeof(float));
+ EMGLLF_core(phia,rhoInit,piInit,gamInit,mini,maxi,gamma,0.0,Xa,Y,tau,