NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ get_timestep() [2/2]

template<typename T >
void cntr::herm_matrix< T >::get_timestep ( int  tstp,
herm_matrix< T > &  timestep 
) const

Definition at line 2179 of file cntr_herm_matrix_impl.hpp.

2179  {
2180 
2181  assert(tstp >= -1 && tstp <= nt_ && tstp <= g1.nt() && "tstp >= -1 && tstp <= nt_ && tstp <= g1.nt()");
2182  assert(g1.size1() == size1_ && "g1.size1() == size1_");
2183  assert(g1.ntau() == ntau_ && "g1.ntau() == ntau_");
2184  if (tstp == -1) {
2185  memcpy(g1.mat_, mat_, sizeof(cplx) * (ntau_ + 1) * element_size_);
2186  } else {
2187  memcpy(g1.retptr(tstp, 0), retptr(tstp, 0),
2188  sizeof(cplx) * (tstp + 1) * element_size_);
2189  memcpy(g1.tvptr(tstp, 0), tvptr(tstp, 0),
2190  sizeof(cplx) * (ntau_ + 1) * element_size_);
2191  memcpy(g1.lesptr(0, tstp), lesptr(0, tstp),
2192  sizeof(cplx) * (tstp + 1) * element_size_);
2193  }
2194 }
std::complex< T > cplx