- // phiInit
- const int dimPhiInit[] = {p, m, k};
- float* phiInit = readArray_real("phiInit",dimPhiInit,3);
-
- // rhoInit
- const int dimRhoInit[] = {m, m, k};
- float* rhoInit = readArray_real("rhoInit",dimRhoInit,3);
-
- // piInit
- float* piInit = readArray_real("piInit",&k,1);
-
- // gamInit
- const int dimGamInit[] = {n, k};
- float* 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
- float* pgamma = readArray_real("gamma",&lengthOne,1);
- float gamma = *pgamma;
- free(pgamma);
-
- // lambda
- float* glambda = readArray_real("glambda",&L,1);
-
- // X
- const int dimX[] = {n, p};
- float* X = readArray_real("X",dimX,2);
-
- // Y
- const int dimY[] = {n, m};
- float* Y = readArray_real("Y",dimY,2);
-
- // seuil
- float* pseuil = readArray_real("seuil",&lengthOne,1);
- float seuil = *pseuil;
- free(pseuil);
-
- // tau
- float* ptau = readArray_real("tau",&lengthOne,1);
- float tau = *ptau;
- free(ptau);
-