triqs::atom_diag::trace_rho_op

#include <triqs/atom_diag/functions.hpp>

Synopsis

template<bool Complex>

typename atom_diag<Complex>::scalar_t trace_rho_op (typename atom_diag<Complex>::block_matrix_t const & density_matrix,

typename atom_diag<Complex>::many_body_op_t const & op,

atom_diag<Complex> const & atom)

Compute Tr (op * density_matrix)

Parameters

• density_matrix Density matrix as a list of diagonal blocks for all invariant subspaces in atom.

• op Operator to be averaged.

• atom Solved diagonalization problem.

Returns

Operator op averaged over the density matrix.