NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

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

1167  {
1168  assert(tstp == tstp_);
1169  if(tstp_==-1){
1170  get_mat(ntau_,M);
1171  M *= (-1.0);
1172  }else{
1173  get_les_tstp_t(tstp_,M);
1174  M *= std::complex<T>(0.0,1.0*sig_);
1175  }
1176  }
+ Here is the call graph for this function: