Function kl_divergence_py 

pub fn kl_divergence_py(
    p: &Bound<'_, PyArray1<f64>>,
    q: &Bound<'_, PyArray1<f64>>,
) -> PyResult<f64>

Compute Kullback-Leibler divergence between two probability distributions