pub fn kl_divergence(p: &[f64], q: &[f64]) -> f64
KL divergence D_KL(P || Q) = sum_x P(x) log(P(x)/Q(x)).