KL function

Kullback-Leibler Divergence