Simple function to calculate Kullback-Leibler divergence