- // phiInit
- const mwSize dimPhiInit[] = {p, m, k};
- Real* phiInit = readArray_real("phiInit",dimPhiInit,3);
-
- // rhoInit
- const mwSize dimRhoInit[] = {m, m, k};
- Real* rhoInit = readArray_real("rhoInit",dimRhoInit,3);
-
- // piInit
- Real* piInit = readArray_real("piInit",&k,1);
-
- // gamInit
- const mwSize dimGamInit[] = {n, k};
- Real* gamInit = readArray_real("gamInit",dimGamInit, 2);
-
- // min number of iterations
- Int* pmini = readArray_int("mini",&lengthOne,1);
- Int mini = *pmini;
- free(pmini);
-
- // max number of iterations
- Int* pmaxi = readArray_int("maxi",&lengthOne,1);
- Int maxi = *pmaxi;
- free(pmaxi);
-
- // gamma
- Real* pgamma = readArray_real("gamma",&lengthOne,1);
- Real gamma = *pgamma;
- free(pgamma);
-
- // lambda
- Real* glambda = readArray_real("glambda",&L,1);
-
- // X
- const mwSize dimX[] = {n, p};
- Real* X = readArray_real("X",dimX,2);
-
- // Y
- const mwSize dimY[] = {n, m};
- Real* Y = readArray_real("Y",dimY,2);
-
- // seuil
- Real* pseuil = readArray_real("seuil",&lengthOne,1);
- Real seuil = *pseuil;
- free(pseuil);
-
- // tau
- Real* ptau = readArray_real("tau",&lengthOne,1);
- Real tau = *ptau;
- free(ptau);
-