correlation_energy() [1/2]

template<typename T >

T cntr::correlation_energy	(	int	tstp,
		herm_matrix< T > &	G,
		herm_matrix< T > &	Sigma,
		T	beta,
		T	h,
		int	SolveOrder
	)

Returns the result of the correlation energy at a given time-step from the time-diagonal convolution

Purpose

Calculates the correlation energy given by \( -i/2\) Tr \([G*\Sigma]^<\) at a given time step 'tstp'. The objects 'G' and 'Sigma' are of the class type 'herm_matrix'.

Parameters

tstp	[int] time step
G	[herm_matrix] contour Green's function
Sigma	[herm_matrix] self-energy
beta	inversed temperature
h	time interval
SolveOrder	[int] integration order

Definition at line 533 of file cntr_utilities_impl.hpp.

References cntr::herm_matrix< T >::size1().

                                          {
  int size1=G.size1();
  std::complex<T> trGxSGM;
  std::complex<T> *GxSGM = new std::complex<T>[size1*size1];

  convolution_density_matrix<T,herm_matrix<T>>(tstp, GxSGM, Sigma, G, integration::I<T>(SolveOrder), beta, h);
  trGxSGM = (0.0,0.0);
  for(int i=0; i< size1; i++){
    trGxSGM += GxSGM[i*size1 + i];
  }
  return 0.5*trGxSGM.real();

  delete GxSGM;
}

Here is the call graph for this function: