p.setp(nsweep=20000000, nmonte=50, ncarlo=50)
# verify the parameters
p.check()
# generate the solver.ctqmc.in file
p.write()
# destroy the instance
del p
# mpi barrier
comm.Barrier()
# init ctqmc impurity solver
ctqmc.cat_init_ctqmc(comm.rank, comm.size)
# setup hybridization function
# for the first iteration, we just use the default value
if i > 0: ctqmc.cat_set_hybf(size_t, hybf)
# exec ctqmc impurity solver
ctqmc.cat_exec_ctqmc(i + 1)
# get green's function
grnf = ctqmc.cat_get_grnf(size_t)
# get new hybridization function
hybf = do_dmft_loop(mfreq, norbs, grnf)
# get occupation number
p.setp(nsweep = 20000000, nmonte = 50, ncarlo = 50)
# verify the parameters
p.check()
# generate the solver.ctqmc.in file
p.write()
# destroy the instance
del p
# mpi barrier
comm.Barrier()
# init ctqmc impurity solver
ctqmc.cat_init_ctqmc(comm.rank, comm.size)
# setup hybridization function
# for the first iteration, we just use the default value
if i > 0: ctqmc.cat_set_hybf(size_t, hybf)
# exec ctqmc impurity solver
ctqmc.cat_exec_ctqmc(i+1)
# get green's function
grnf = ctqmc.cat_get_grnf(size_t)
# get new hybridization function
hybf = do_dmft_loop(mfreq, norbs, grnf)
# get occupation number