dyson_start()

template<typename T >

void cntr::dyson_start	(	herm_matrix< T > &	G,
		T	mu,
		function< T > &	H,
		herm_matrix< T > &	Sigma,
		T	beta,
		T	h,
		const int	SolveOrder
	)

Start-up procedure for solving the Dyson equation of the integral-differential form for a Green's function \(G\)

Purpose

One solves the Dyson equation of the following form: \( [ id/dt + \mu - H(t) ] G(t,t^\prime) - [\Sigma*G](t,t^\prime) = \delta(t,t^\prime)\) for a hermitian matrix \(G(t, t^\prime)\) for the first k timesteps (given by the integrator class 'I'). One assumes that the Matsubara component of \(G\) and \(\Sigma(t,t^\prime)\) for \(t,t^\prime\)<=k are given.

Parameters

&G	[herm_matrix<T>] solution
mu	[T] chemical potential
&H	[function<T>] time-dependent function
&Sigma	[herm_matrix<T>] self-energy
beta	[double] inverse temperature
h	[double] time interval
SolveOrder	[int] integrator order

Definition at line 1502 of file cntr_dyson_impl.hpp.

References cntr::herm_matrix< T >::nt(), cntr::herm_matrix< T >::ntau(), cntr::function< T >::ptr(), and cntr::herm_matrix< T >::size1().

                                                    {
    int size1 = G.size1();
    assert(G.size1() == Sigma.size1());
    assert(G.ntau() == Sigma.ntau());
    assert(G.nt() >= SolveOrder);
    assert(Sigma.nt() >= SolveOrder);
    if (size1 == 1) {
        dyson_start_ret<T, herm_matrix<T>, 1>(G, mu, H.ptr(0), Sigma, integration::I<T>(SolveOrder), h);
        dyson_start_tv<T, herm_matrix<T>, 1>(G, mu, H.ptr(0), Sigma, integration::I<T>(SolveOrder), beta, h);
        dyson_start_les<T, herm_matrix<T>, 1>(G, mu, H.ptr(0), Sigma, integration::I<T>(SolveOrder), beta, h);
    } else {
        dyson_start_ret<T, herm_matrix<T>, LARGESIZE>(G, mu, H.ptr(0), Sigma, integration::I<T>(SolveOrder), h);
        dyson_start_tv<T, herm_matrix<T>, LARGESIZE>(G, mu, H.ptr(0), Sigma, integration::I<T>(SolveOrder), beta, h);
        dyson_start_les<T, herm_matrix<T>, LARGESIZE>(G, mu, H.ptr(0), Sigma, integration::I<T>(SolveOrder), beta, h);
    }
}

Here is the call graph for this function: