% IFAM Dual de Godel

clear;
load facesdual;

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

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

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