pub fn diagonal_density_matrix( eigenvectors: &[Vec<Complex64>], probabilities: &[f64], ) -> Result<Vec<Complex64>>