function D = relent(p,q)
% Compute the relative entropy or Kullback-Leibler
% informational divergence sum_i p(i) log2(p(i)/q(i)).
D = entropy2(q./p,p);