NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ 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().

1102  {
1103  assert(tstp_ == tstp);
1104  cplx x1;
1105  if (tstp_ == -1) {
1106  get_mat(ntau_, x1);
1107  rho = -x1;
1108  } else {
1109  get_les(tstp_, tstp_, x1);
1110  rho = std::complex<T>(0.0, sig_) * x1;
1111  }
1112  }
+ Here is the call graph for this function: