NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library
cntr_getset_herm_matrix_timestep_view_inc.hpp
Go to the documentation of this file.
1 
22 template <typename T>
23 inline void get_mat(const int m, std::complex<T> &G_mat, herm_matrix_timestep_view<T> &G,
24  herm_matrix_timestep_view<T> &Gcc){
25  assert(m <= G.ntau());
26  assert(G.ntau() == Gcc.ntau());
27  assert(G.size1() == Gcc.size1());
28  assert(G.size2() == Gcc.size2());
29  assert(G.sig() == Gcc.sig());
30 
31  G_mat = *G.matptr(m);
32 }
33 
55 template <typename T>
56 void get_mat(const int m, cdmatrix &G_mat, herm_matrix_timestep_view<T> &G,
57  herm_matrix_timestep_view<T> &Gcc){
58  assert(m <= G.ntau());
59  assert(G.ntau() == Gcc.ntau());
60  assert(G.size1() == Gcc.size1());
61  assert(G.size2() == Gcc.size2());
62  assert(G.sig() == Gcc.sig());
63  int size1=G.size1(), size2=G.size2();
64  std::complex<T> *mat;
65 
66  mat = G.matptr(m);
67  map_ptr2matrix(size1, size2, mat, G_mat);
68 }
69 
70 
71 
91 template <typename T>
92 inline void get_mat(const int m, std::complex<T> &G_mat, herm_matrix_timestep_view<T> &G){
93  assert(m <= G.ntau());
94 
95  G_mat = *G.matptr(m);
96 }
97 
117 template <typename T>
118 void get_mat(const int m, cdmatrix &G_mat, herm_matrix_timestep_view<T> &G){
119  assert(m <= G.ntau());
120  std::complex<T> *mat;
121  int size1=G.size1(), size2=G.size2();
122 
123  mat = G.matptr(m);
124  map_ptr2matrix(size1, size2, mat, G_mat);
125 }
126 
127 
128 
153 template <typename T>
154 inline void get_les(const int i, const int j, std::complex<T> &G_les, herm_matrix_timestep_view<T> &G,
155  herm_matrix_timestep_view<T> &Gcc){
156  assert(i <= G.tstp() && j <= G.tstp());
157  assert(i == G.tstp() || j == G.tstp());
158  assert(i <= Gcc.tstp() && j <= Gcc.tstp());
159  assert(i == Gcc.tstp() || j == Gcc.tstp());
160  assert(G.tstp() == Gcc.tstp());
161  assert(G.ntau() == Gcc.ntau());
162  assert(G.size1() == Gcc.size1());
163  assert(G.size2() == Gcc.size2());
164  assert(G.sig() == Gcc.sig());
165 
166  if (G.tstp() == j){
167  G_les = *G.lesptr(i);
168  } else if (G.tstp() == i) {
169  G_les = *Gcc.lesptr(j);
170  G_les = -std::conj(G_les);
171  }
172 }
173 
198 template <typename T>
199 void get_les(const int i, const int j, cdmatrix &G_les, herm_matrix_timestep_view<T> &G,
200  herm_matrix_timestep_view<T> &Gcc){
201  assert(i <= G.tstp() && j <= G.tstp());
202  assert(i == G.tstp() || j == G.tstp());
203  assert(i <= Gcc.tstp() && j <= Gcc.tstp());
204  assert(i == Gcc.tstp() || j == Gcc.tstp());
205  assert(G.tstp() == Gcc.tstp());
206  assert(G.ntau() == Gcc.ntau());
207  assert(G.size1() == Gcc.size1());
208  assert(G.size2() == Gcc.size2());
209  assert(G.sig() == Gcc.sig());
210  std::complex<T> *les;
211  int size1=G.size1(), size2=G.size2();
212 
213  if (G.tstp() == j){
214  les = G.lesptr(i);
215  map_ptr2matrix(size1, size2, les, G_les);
216  } else if (G.tstp() == i) {
217  les = Gcc.lesptr(j);
218  map_ptr2matrix(size1, size2, les, G_les);
219  G_les.adjointInPlace();
220  G_les = -G_les;
221  }
222 }
223 
224 
225 
251 template <typename T>
252 inline void get_ret(const int i, const int j, std::complex<T> &G_ret, herm_matrix_timestep_view<T> &G,
253  herm_matrix_timestep_view<T> &Gcc){
254  assert(i <= G.tstp() && j <= G.tstp());
255  assert(i == G.tstp() || j == G.tstp());
256  assert(i <= Gcc.tstp() && j <= Gcc.tstp());
257  assert(i == Gcc.tstp() || j == Gcc.tstp());
258  assert(G.tstp() == Gcc.tstp());
259  assert(G.ntau() == Gcc.ntau());
260  assert(G.size1() == Gcc.size1());
261  assert(G.size2() == Gcc.size2());
262  assert(G.sig() == Gcc.sig());
263 
264  if (G.tstp() == i){
265  G_ret = *G.retptr(j);
266  } else if (G.tstp() == j) {
267  G_ret = *Gcc.retptr(i);
268  G_ret = -std::conj(G_ret);
269  }
270 }
271 
297 template <typename T>
298 void get_ret(const int i, const int j, cdmatrix &G_ret, herm_matrix_timestep_view<T> &G,
299  herm_matrix_timestep_view<T> &Gcc){
300  assert(i <= G.tstp() && j <= G.tstp());
301  assert(i == G.tstp() || j == G.tstp());
302  assert(i <= Gcc.tstp() && j <= Gcc.tstp());
303  assert(i == Gcc.tstp() || j == Gcc.tstp());
304  assert(G.tstp() == Gcc.tstp());
305  assert(G.ntau() == Gcc.ntau());
306  assert(G.size1() == Gcc.size1());
307  assert(G.size2() == Gcc.size2());
308  assert(G.sig() == Gcc.sig());
309  std::complex<T> *ret;
310  int size1=G.size1(), size2=G.size2();
311 
312  if (G.tstp() == i){
313  ret = G.retptr(j);
314  map_ptr2matrix(size1, size2, ret, G_ret);
315  } else if (G.tstp() == j) {
316  ret = Gcc.retptr(i);
317  map_ptr2matrix(size1, size2, ret, G_ret);
318  G_ret.adjointInPlace();
319  G_ret = -G_ret;
320  }
321 }
322 
323 
324 
325 
350 template <typename T>
351 inline void get_tv(const int i, const int m, std::complex<T> &G_tv, herm_matrix_timestep_view<T> &G,
352  herm_matrix_timestep_view<T> &Gcc){
353  assert(i == G.tstp());
354  assert(m <= G.ntau());
355  assert(G.tstp() == Gcc.tstp());
356  assert(G.ntau() == Gcc.ntau());
357  assert(G.size1() == Gcc.size1());
358  assert(G.size2() == Gcc.size2());
359  assert(G.sig() == Gcc.sig());
360 
361  G_tv = *G.tvptr(m);
362 }
363 
388 template <typename T>
389 void get_tv(const int i, const int m, cdmatrix &G_tv, herm_matrix_timestep_view<T> &G, herm_matrix_timestep_view<T> &Gcc){
390  assert(i == G.tstp());
391  assert(m <= G.ntau());
392  assert(G.tstp() == Gcc.tstp());
393  assert(G.ntau() == Gcc.ntau());
394  assert(G.size1() == Gcc.size1());
395  assert(G.size2() == Gcc.size2());
396  assert(G.sig() == Gcc.sig());
397  int size1 = G.size1(), size2 = G.size2();
398  std::complex<T> *tv;
399 
400  tv = G.tvptr(m);
401  map_ptr2matrix(size1, size2, tv, G_tv);
402 }
403 
404 
405 
406 
407 
408 
409 
434 template <typename T>
435 inline void get_vt(const int m, const int i, std::complex<T> &G_vt, herm_matrix_timestep_view<T> &G,
436  herm_matrix_timestep_view<T> &Gcc){
437  assert(i == G.tstp());
438  assert(m <= G.ntau());
439  assert(G.tstp() == Gcc.tstp());
440  assert(G.ntau() == Gcc.ntau());
441  assert(G.size1() == Gcc.size1());
442  assert(G.size2() == Gcc.size2());
443  assert(G.sig() == Gcc.sig());
444 
445  G_vt = *Gcc.tvptr(G.ntau() - m);
446  G_vt = std::complex<T>(-G.sig(),0.0) * std::conj(G_vt);
447 }
448 
473 template <typename T>
474 void get_vt(const int m, const int i, cdmatrix &G_vt, herm_matrix_timestep_view<T> &G,
475  herm_matrix_timestep_view<T> &Gcc){
476  assert(i == G.tstp());
477  assert(m <= G.ntau());
478  assert(G.tstp() == Gcc.tstp());
479  assert(G.ntau() == Gcc.ntau());
480  assert(G.size1() == Gcc.size1());
481  assert(G.size2() == Gcc.size2());
482  assert(G.sig() == Gcc.sig());
483  int size1 = G.size1(), size2 = G.size2();
484  std::complex<T> *vt;
485 
486  vt = Gcc.tvptr(G.ntau() - m);
487  map_ptr2matrix(size1, size2, vt, G_vt);
488  G_vt.adjointInPlace();
489  G_vt = std::complex<T>(-G.sig(),0.0) * G_vt;
490 }
491 
492 
493 
494 
495 
519 template <typename T>
520 inline void get_gtr(const int i, const int j, std::complex<T> &G_gtr, herm_matrix_timestep_view<T> &G,
521  herm_matrix_timestep_view<T> &Gcc){
522  std::complex<T> G_les, G_ret;
523 
524  get_les(i, j, G_les, G, Gcc);
525  get_ret(i, j, G_ret, G, Gcc);
526  G_gtr = G_les + G_ret;
527 }
528 
553 template <typename T>
554 void get_gtr(const int i, const int j, cdmatrix &G_gtr, herm_matrix_timestep_view<T> &G,
555  herm_matrix_timestep_view<T> &Gcc){
556  cdmatrix G_les, G_ret;
557 
558  get_les(i, j, G_les, G, Gcc);
559  get_ret(i, j, G_ret, G, Gcc);
560  G_gtr.noalias() = G_les + G_ret;
561 }
562 
563 
564 
565 
566 
567 
568 
591 template <typename T>
592 inline void get_les(const int i, const int j, std::complex<T> &G_les, herm_matrix_timestep_view<T> &G){
593  get_les(i, j, G_les, G, G);
594 }
595 
618 template <typename T>
619 void get_les(const int i, const int j, cdmatrix &G_les, herm_matrix_timestep_view<T> &G){
620  get_les(i, j, G_les, G, G);
621 }
622 
623 
624 
648 template <typename T>
649 inline void get_ret(const int i, const int j, std::complex<T> &G_ret, herm_matrix_timestep_view<T> &G){
650  get_ret(i, j, G_ret, G, G);
651 }
652 
676 template <typename T>
677 void get_ret(const int i, const int j, cdmatrix &G_ret, herm_matrix_timestep_view<T> &G){
678  get_ret(i, j, G_ret, G, G);
679 }
680 
681 
682 
683 
706 template <typename T>
707 inline void get_tv(const int i, const int m, std::complex<T> &G_tv, herm_matrix_timestep_view<T> &G){
708  get_tv(i, m, G_tv, G);
709 }
710 
733 template <typename T>
734 void get_tv(const int i, const int m, cdmatrix &G_tv, herm_matrix_timestep_view<T> &G){
735  get_tv(i, m, G_tv, G);
736 }
737 
738 
739 
740 
741 
742 
743 
766 template <typename T>
767 inline void get_vt(const int m, const int i, std::complex<T> &G_vt, herm_matrix_timestep_view<T> &G){
768  get_vt(m, i, G_vt, G);
769 }
770 
793 template <typename T>
794 void get_vt(const int m, const int i, cdmatrix &G_vt, herm_matrix_timestep_view<T> &G){
795  get_vt(m, i, G_vt, G);
796 }
797 
798 
799 
800 
801 
823 template <typename T>
824 inline void get_gtr(const int i, const int j, std::complex<T> &G_gtr, herm_matrix_timestep_view<T> &G){
825  get_gtr(i, j, G_gtr, G, G);
826 }
827 
850 template <typename T>
851 void get_gtr(const int i, const int j, cdmatrix &G_gtr, herm_matrix_timestep_view<T> &G){
852  get_gtr(i, j, G_gtr, G, G);
853 }
get_tv
void get_tv(const int i, const int m, std::complex< T > &G_tv, herm_matrix_timestep_view< T > &G, herm_matrix_timestep_view< T > &Gcc)
Returns the left-mixing component of a general contour function at given times. ...
Definition: cntr_getset_herm_matrix_timestep_view_inc.hpp:351
get_vt
void get_vt(const int m, const int i, std::complex< T > &G_vt, herm_matrix_timestep_view< T > &G, herm_matrix_timestep_view< T > &Gcc)
Returns the right-mixing component of a general contour function at given times. ...
Definition: cntr_getset_herm_matrix_timestep_view_inc.hpp:435
get_mat
void get_mat(const int m, std::complex< T > &G_mat, herm_matrix_timestep_view< T > &G, herm_matrix_timestep_view< T > &Gcc)
Returns the Matsubara component of a general contour function at given times.
Definition: cntr_getset_herm_matrix_timestep_view_inc.hpp:23
get_gtr
void get_gtr(const int i, const int j, std::complex< T > &G_gtr, herm_matrix_timestep_view< T > &G, herm_matrix_timestep_view< T > &Gcc)
Returns the greater component of a general contour function at given times.
Definition: cntr_getset_herm_matrix_timestep_view_inc.hpp:520
get_les
void get_les(const int i, const int j, std::complex< T > &G_les, herm_matrix_timestep_view< T > &G, herm_matrix_timestep_view< T > &Gcc)
Returns the lesser component of a general contour function at given times.
Definition: cntr_getset_herm_matrix_timestep_view_inc.hpp:154
get_ret
void get_ret(const int i, const int j, std::complex< T > &G_ret, herm_matrix_timestep_view< T > &G, herm_matrix_timestep_view< T > &Gcc)
Returns the retarded component of a general contour function at given times.
Definition: cntr_getset_herm_matrix_timestep_view_inc.hpp:252