from pytriqs.base.gf_local import GfReFreq
from pytriqs.base.archive import HDFArchive
from math import pi
R = HDFArchive('myfile.h5', 'r')
from pytriqs.base.plot.mpl_interface import oplot, plt
plt.xrange(-1,1)
plt.yrange(0,7)
for name, g in R.items() : # iterate on the elements of R, like a dict ...
oplot( (- 1/pi * g).imag, "-o", name = name)
p.savefig("./tut_ex3b.png")
from pytriqs.base.gf_local import *
from pytriqs.base.gf_local.descriptors import iOmega_n
g = GfImFreq(indices = [1], beta = 300, n_matsubara = 1000, name = "g")
g <<= inverse( iOmega_n + 0.5 )
print " van plot"
oplot (g, '-o', x_window = (0,3) )
print "plot done"
g<<= inverse( iOmega_n + 0.5 )
print "ok ----------------------"
from pytriqs.base.archive import HDFArchive
R = HDFArchive('myfile.h5', 'r')
for n, calculation in R.items() :
#g = calculation['g']
g <<= inverse( iOmega_n + 0.5 )
print "pokokook"
X,Y = g.x_data_view (x_window = (0,0.2), flatten_y = True )
#fitl = Fit ( X,Y.imag, linear )
g <<= inverse( iOmega_n + 0.5 )
print " van plot"
oplot (g, '-o', x_window = (0,3) )
g <<= inverse( iOmega_n + 0.5 )
