NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ set_submatrix() [2/2]

template<typename T >
void cntr::herm_matrix_timestep< T >::set_submatrix ( int  tstp,
std::vector< int > &  i1,
std::vector< int > &  i2,
herm_matrix_timestep< T > &  g,
std::vector< int > &  j1,
std::vector< int > &  j2 
)

Sets a (sub-) matrix of this contour object at a given time step to a (sub-) matrix of a given herm_matrix_timestep.

Purpose

Sets the (sub-) matrix of a two-time contour object \(C(t,t')\) the (sub-) of a given herm_matrix_timestep \(g\) according to \( C_{i_1(k) i_2(k)}(t,t') = g_{j_1(k) j_2(k)}(t,t') \) with \(k\) denoting an iterator of a subspace (at a given time step).

Parameters
tstp

Time step tstp.

i1

Vector of row indices of C.

i2

Vector of column indices of C.

g

The herm_matrix_timestep from which the (sub-) matrix is copied.

j1

Vector of row indices of g.

j2

Vector of row indices of g.

Definition at line 2400 of file cntr_herm_matrix_timestep_impl.hpp.

References cntr::herm_matrix_timestep< T >::tstp().

2401  {
2402  assert(tstp == tstp_);
2403  assert(tstp == g.tstp());
2404  assert(i1.size()==i2.size() && i1.size()==j1.size() && j1.size()==j2.size());
2405  assert(size1_*size2_==i1.size());
2406  for(int k1=0;k1<i1.size();k1++){
2407  this->set_matrixelement(tstp,i1[k1],i2[k1],g,j1[k1],j2[k1]);
2408  }
2409 }
void set_matrixelement(int i1, int i2, herm_matrix_timestep_view< T > &g, int j1, int j2)
Set the matrix element of for each component from
+ Here is the call graph for this function: