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
