disp init;

load a3faces.mat;
load a3testfaces.mat;

randn('state',sum(100*clock));

[d, n_train] = size(xx_train);
[d, n_test] = size(xx_test);

mu = mean(xx_train, 2);
s = cov(xx_train', 1);

xm_train = zeros(size(xx_train));
for i = 1:n_train
    xm_train(:,i) = xx_train(:,i) - mu;
end
xm_test = zeros(size(xx_test));
for i = 1:n_test
    xm_test(:,i) = xx_test(:,i) - mu;
end

k = 2; % number of latent variables