NESSi  v1.0.2
The NonEquilibrium Systems Simulation Library

◆ convolution_density_matrix() [4/4]

template<typename T , class GG >
void cntr::convolution_density_matrix ( int  tstp,
cdmatrix &  rho,
GG &  A,
GG &  B,
beta,
h,
int  SolveOrder 
)

Returns the result of the contour convolution for a density matrix

Purpose

Calls convolution routines to compute contour convolution \(\rho=-i (A*B)^<\). The objects A and B are of the class type 'GG'. Works for a scalar and square matrices.

Parameters
tstp[int] given time-step
rho

[std::complex] Matrix to which the result of the convolution is given

A

[GG] contour Green's function

B

[GG] contour Green's function

beta

inversed temperature

h

time interval

SolveOrder

[int] integrator ordder

Definition at line 4073 of file cntr_convolution_impl.hpp.

4074  {
4075  int size1 = A.size1();
4076  int element_size = size1;
4077  std::complex<T> *rho_ptr;
4078  rho_ptr = new std::complex<T>[element_size];
4079 
4080  convolution_density_matrix<T, GG>(tstp, rho_ptr, A, A, B, B, integration::I<T>(SolveOrder), beta, h);
4081  map_ptr2matrix<T>(size1, size1, rho_ptr, rho);
4082 
4083  delete[] rho_ptr;
4084 }