from pytriqs.gf.local import GfReFreq from pytriqs.archive import HDFArchive from math import pi R = HDFArchive('myfile.h5', 'r') from pytriqs.plot.mpl_interface import oplot, plt for name, g in R.items() : # iterate on the elements of R, like a dict ... oplot( (- 1/pi * g).imag, "-o", name = name) plt.xlim(-1,1) plt.ylim(0,7) p.savefig("./tut_ex3b.png")
from pytriqs.gf.local import GfReFreq from pytriqs.archive import HDFArchive from math import pi R = HDFArchive('myfile.h5', 'r') from pytriqs.plot.mpl_interface import oplot, plt for name, g in R.items(): # iterate on the elements of R, like a dict ... oplot((-1 / pi * g).imag, "-o", name=name) plt.xlim(-1, 1) plt.ylim(0, 7) p.savefig("./tut_ex3b.png")
g = GfImFreq(indices = [1], beta = 300, n_points = 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.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
from pytriqs.gf.descriptors import iOmega_n g = GfImFreq(indices=[1], beta=300, n_points=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.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