flow implemented, greens functions and initial conditions still missing
optimization the matsubara sums to be done ( to many check)
add precalc G x S product