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