To estimate the model several times
without orthogonality constraints, do the following:
T1=[];
for i=1:5
tic
f1=parafac(X,5);
t1=toc;
T1=[T1 t1];
end
To do the same using orthogonality
e.g., in the emission mode do
T2=[];
for i=1:5
tic
f2=parafac(X,5,[],[0 1 0]);
t2=toc;
T2=[T2 t2];
end