from triqs.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 h5 import HDFArchive
R = HDFArchive('myfile.h5', 'r')
for n, calculation in list(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 triqs.gf import GfReFreq
from h5 import HDFArchive
from math import pi
R = HDFArchive('myfile.h5', 'r')
from triqs.plot.mpl_interface import oplot, plt
for name, g in list(
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")
