RunAlgorithmHOSVD.m

function percentage = RunAlgorithmHOSVD(train, test, gt)
    %TODO: optimize - lots of redundant operations 
    %input:
        %train: train dataset
        %test: test dataset
        %gt: ground truth; indexes of the correct person

    n_samples = size(test, 1);
    A = prepareTensor3D(train);
    [S, F, G, H] = HOSVD(A);
    C = ten_mat_mult(ten_mat_mult(S, F, 1), G, 2);
  
    count = 0;
    for(i = 1:n_samples)
       z = test(i, :)';
       person = HOSVDmethod(C, H, z);
       
       if(person == gt(i))
           count = count + 1;
       end
       fprintf('Current percentage after iter %s is %s %% \n', num2str(i), num2str(count/i * 100));
    end
    percentage = count / n_samples;
    fprintf('Final percentage of correct recognitions: %s %% \n', num2str(percentage * 100));
end