pub fn matrix_element( left: &[Complex64], operator: &(impl LinearOperator + ?Sized), right: &[Complex64], ) -> Result<Complex64>