Added an implementation of Kullback-Leibler divergence calculation (kl_divergence)