convolution_density_matrix() [3/4]

template<typename T , class GG >

void cntr::convolution_density_matrix	(	int	n,
		cdmatrix &	rho,
		GG &	A,
		GG &	Acc,
		GG &	B,
		GG &	Bcc,
		T	beta,
		T	h,
		int	SolveOrder
	)

Returns the result of the contour convolution for a density matrix at a given time-step

Purpose

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

Parameters

n	[int] number of the time step ('t=nh')
rho	[std::complex] Matrix to which the result of the convolution is given
A	[GG] contour Green's function
Acc	[GG] complex conjugate to A
B	[GG] contour Green's function
Bcc	[GG] complex conjugate to B
beta	inversed temperature
h	time interval
SolveOrder	[int] integrator order

Definition at line 3990 of file cntr_convolution_impl.hpp.

                                 {

    int size1 = A.size1();
    int element_size = size1;
    std::complex<T> *rho_ptr;
    rho_ptr = new std::complex<T>[element_size];

    convolution_density_matrix<T, GG>(n, rho_ptr, A, Acc, B, Bcc, integration::I<T>(SolveOrder), beta, h);
    map_ptr2matrix<T>(size1, size1, rho_ptr, rho);

    delete[] rho_ptr;

}