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")
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 ) print g