function [ d ] = mD_IS( A, B )
%MD_IS matrix difference using Itakura-Saito cost function
%
%   (c) Joachim Thiemann 2013    joachim.thiemann@ieee.org

rat = A ./ B;
is = rat - log(rat) - 1;
d = sum(is(:));

end
