pub fn posterior(psi: &Psi, w: &Col<f64>) -> Result<Mat<f64>>
Calculates the posterior probabilities for each support point given the weights