density_matrix()

template<typename T >

template<class Matrix >

void cntr::herm_matrix< T >::density_matrix	(	int	tstp,
		Matrix &	M
	)		const

Returns the density matrix at given time step.

Purpose

Returns the matrix-valued density matrix at given time step tstp. Setting tstp = -1 returns the equilibrium density matrix \( \rho = -C^\mathrm{M}(\beta) \), while tstp >= 0 returns \( \rho(t) = i \eta C^<(t,t) \). Works for square-matrices only.

Parameters

tstp	The time step at which the density matrix is returned.
M	The density matrix at time step tstp.

Definition at line 1018 of file cntr_herm_matrix_impl.hpp.

References cntr::get_les(), and cntr::get_mat().

                                                             {
  assert(M.rows() == size1_ && M.cols() == size2_);
  assert(tstp >= -1 && tstp <= nt_);
    if (tstp == -1) {
        get_mat(ntau_, M);
        M *= (-1.0);
    } else {
        get_les(tstp, tstp, M);
        M *= std::complex<T>(0.0, 1.0 * sig_);
    }
}

Here is the call graph for this function: