density_matrix() [1/2]

template<typename T >

template<class Matrix >

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

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	given time step
M	The density matrix at time step tstp.

Definition at line 1167 of file cntr_herm_matrix_timestep_impl.hpp.

References cntr::get_mat().

                                                                     {
         assert(tstp == tstp_);
         if(tstp_==-1){
            get_mat(ntau_,M);
            M *= (-1.0);
         }else{
            get_les_tstp_t(tstp_,M);
            M *= std::complex<T>(0.0,1.0*sig_);
         }
      }

Here is the call graph for this function: