% IFAM Dual de Lukasiewicz...

clear;
load facesdual;

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

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

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