% Recall using max-prod FAM...

clear;
load faces;

W = maxC_P(X,X'); 

yGFAM_P=maxC_P(W,Xin);
disp('max-prod GFAM test - concluded');

figure(3); clf;
for i=1:size(yGFAM_P,2)
	subplot(4,6,i); imshow(reshape(yGFAM_P(:,i),s1,s2));
	disp(sprintf('%.4f',norm(yGFAM_P(:,i)-X(:,1))/norm(X(:,1))));
end
print -deps facesGFAM_P.eps
