density_matrix() [2/2]

template<typename T >

void cntr::herm_matrix_timestep< T >::density_matrix ( int  tstp, cplx &  rho  )

inline

Returns the density matrix at given time step.

Purpose

Returns the scalar-valued density matrix (occupation, that is) 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) \). The return value is formally complex.

Parameters

tstp	The time step at which the density matrix is returned.
rho	Scalar-valued density matrix (complex number)

Definition at line 1102 of file cntr_herm_matrix_timestep_impl.hpp.

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

                                                                             {
         assert(tstp_ == tstp);
         cplx x1;
         if (tstp_ == -1) {
            get_mat(ntau_, x1);
            rho = -x1;
         } else {
            get_les(tstp_, tstp_, x1);
            rho = std::complex<T>(0.0, sig_) * x1;
         }
      }

Here is the call graph for this function: