pub fn symmetric_kl(p: &[f64], q: &[f64]) -> f64
Symmetric KL divergence: (D_KL(P||Q) + D_KL(Q||P)) / 2