N_Legendre_Coeffs=50, # Number of Legendre coefficients
Random_Generator_Name='mt19937', # Name of the random number generator
Use_Segment_Picture=True, # Use the segment picture
Measured_Operators={ # Operators to be averaged
'Nimp': N('up', 1) + N('down', 1)
},
Global_Moves=[ # Global move in the QMC
(0.05, lambda (a, alpha, dag): ({
'up': 'down',
'down': 'up'
}[a], alpha, dag))
],
)
# Initialize the non-interacting Green's function S.G0
for spin, g0 in S.G0:
g0 <<= inverse(iOmega_n - e_f - V**2 * Wilson(D))
# Run the solver. The result will be in S.G
S.Solve()
# Save the results in an hdf5 file (only on the master node)
from pytriqs.Base.Archive import HDF_Archive
import pytriqs.Base.Utility.MPI as MPI
if MPI.IS_MASTER_NODE():
Results = HDF_Archive("solution.h5", 'w')
Results["G"] = S.G
Results["Gl"] = S.G_Legendre
Results["Nimp"] = S.Measured_Operators_Results['Nimp']
