% IFAM Dual de Godel

clear;
load facesdual;

M = maxJ_P(X,X'); theta=max(X,[],2);

% Para efeitos computacionais!
Theta=[]; 
for i=1:size(Xin,2)
	Theta=[Theta,theta];
end
	
yIFAMd_P=min(minD_P(M,Xin),Theta);

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