NESSi
v1.0.2
The NonEquilibrium Systems Simulation Library
|
void cntr::dyson_timestep | ( | int | n, |
herm_matrix< T > & | G, | ||
T | mu, | ||
function< T > & | H, | ||
herm_matrix< T > & | Sigma, | ||
T | beta, | ||
T | h, | ||
const int | SolveOrder | ||
) |
One step Dyson solver (integral-differential form) for a Green's function \(G\)
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)\) at a given timestep 'n', i.e., G^ret(nh,t'<=nh), G^les(t<=nh,nh), G^tv(nt,tau=0..beta). Timestep must be >k, where k is the Integration order 'I'. The timesteps n=0..k must be computed seperately, using the routine "_start", which assumes that the Matsubara component of \(G\) and \(\Sigma(t,t^\prime)\) for \(t,t^\prime\)<=k are given. Here, are given: \(\Sigma(t, t^\prime)\), \(\mu\), and \(H(t)\).
n |
|
&G |
|
mu |
|
&H |
|
&Sigma |
|
beta |
|
h |
|
SolveOrder |
|
Definition at line 1618 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().