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