pub fn diagonal_ensemble( eigenvalues: &[f64], eigenvectors: &[Vec<Complex64>], initial: &[Complex64], ) -> Result<DiagonalEnsemble>