function weightsLDA = getLDAweights(matrixData,clabels)

% Feature extraction based on LDA
%  
% Inputs:
%
% matrixData - d*t*n tensor containing n matrices of size d*t, where n is the number of samples.
% clabels - n-dimensional vector of binary coded category labels
%
% Ouput:
% 
% weightsLDA - d*t matrix of t-dimensional directions estimated with LDA
%
%
% Jukka-Pekka Kauppi
% jukka-pekka.kauppi@helsinki.fi
% University of Helsinki, Department of Computer Science
% 23.8.2012

siz = size(matrixData);
d = siz(1);
t = siz(2);
n = siz(3);

weightsLDA = zeros(d,t);

for m = 1:d
  X = squeeze(matrixData(m,:,:))'; % n*t matrix
  weightsLDA(m,:) = computeLDA(X,clabels);  
end