ENH Implement Jensen-Shannon Divergence

Adds two functions: one computes a pairwise divergence, the other a
matrix.

This is used to compute divergences between distributions. Interface is similar to that of ``entropy()``.