kl_batch.m

The missing values are again left out when calculating the Kullback-Leibler divergence.

sizedata = prod(size(fs));

% the reconstruction error
fs_x = fs - x;
if nargin > 4
  fs_x(find(missing)) = probdist(0, 0);
  sizedata = sizedata - sum(sum(missing));
end;