1 function[] = testEMGLLF()
7 dimensions = dlmread(strcat(testFolder,'dimensions'), delimiter);
14 phiInit = reshape(dlmread(strcat(testFolder,'phiInit'), delimiter), p, m, k);
15 rhoInit = reshape(dlmread(strcat(testFolder,'rhoInit'), delimiter), m, m, k);
16 piInit = transpose(dlmread(strcat(testFolder,'piInit'), delimiter));
17 gamInit = reshape(dlmread(strcat(testFolder,'gamInit'), delimiter), n, k);
18 mini = int64(dlmread(strcat(testFolder,'mini'), delimiter));
19 maxi = int64(dlmread(strcat(testFolder,'maxi'), delimiter));
20 gamma = dlmread(strcat(testFolder,'gamma'), delimiter);
21 lambda = dlmread(strcat(testFolder,'lambda'), delimiter);
22 X = reshape(dlmread(strcat(testFolder,'X'), delimiter), n, p);
23 Y = reshape(dlmread(strcat(testFolder,'Y'), delimiter), n, m);
24 tau = dlmread(strcat(testFolder,'tau'), delimiter);
27 [phi,rho,pi,LLF,S] = EMGLLF(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,lambda,X,Y,tau);
29 %get all stored outputs
30 ref_phi = reshape(dlmread(strcat(testFolder,'phi'), delimiter), p, m, k);
31 ref_rho = reshape(dlmread(strcat(testFolder,'rho'), delimiter), m, m, k);
32 ref_pi = dlmread(strcat(testFolder,'pi'), delimiter);
33 ref_LLF = dlmread(strcat(testFolder,'LLF'), delimiter);
34 ref_S = reshape(dlmread(strcat(testFolder,'S'), delimiter), p, m, k);
36 %check that output correspond to stored output
38 checkOutput('phi',phi,ref_phi,tol);
39 checkOutput('rho',rho,ref_rho,tol);
40 checkOutput('pi',pi,ref_pi,tol);
41 checkOutput('LLF',LLF,ref_LLF,tol);
42 checkOutput('S',S,ref_S,tol);