boblsturm@0: function [cost] = KLDivCost(A , B)
boblsturm@0: 
boblsturm@0: % dim=size(A);
boblsturm@0: % cost=0;
boblsturm@0: 
boblsturm@0: cost = sum(sum(A.*log10(A./B)-A+B));
boblsturm@0: 
boblsturm@0: % for i=1:dim(1)
boblsturm@0: %     for j=1:dim(2)
boblsturm@0: %         cost = cost + A(i,j)*log10(A(i,j)/B(i,j))-A(i,j)+B(i,j);
boblsturm@0: %     end
boblsturm@0: % end