#!/usr/bin/env python3
import qlat as q
import sys, os
q.begin()
q.qremove_all_info("results")
q.qmkdir_info("results")
rs = q.RngState("seed")
total_site = [4, 4, 4, 8]
geo = q.Geometry(total_site, 1)
q.displayln_info("geo.show() =", geo.show())
psel = q.PointSelection([[0,0,0,0], [0,1,2,0]])
psel.save("results/psel.txt")
psel = q.PointSelection()
psel.load("results/psel.txt")
q.displayln_info("psel.list() =", psel.list())
n_per_tslice = 16
fsel = q.FieldSelection(geo.total_site(), n_per_tslice, rs)
q.displayln_info("fsel info =", fsel.geo().show(), fsel.n_per_tslice(), fsel.prob())
fsel.save("results/fsel.field")
def begin_with_gpt():
grid = mk_grid()
size_node = grid.mpi
coor_node = grid.processor_coor
id_node = q.index_from_coordinate(coor_node, size_node)
q.begin(id_node, size_node)
"results/gf_info/traj={}.lat".format(traj),
"results/wilson_flow_energy_info/traj={}.lat".format(traj), gf,
0.1, 5, 2)
q.save_gm_force_magnitudes_list(
"results/gm_force_info/traj={}.lat".format(traj))
@q.timer_verbose
def main():
total_site = [4, 4, 4, 8]
ga = q.GaugeAction(2.13, -0.331)
test_hmc(total_site, ga)
ga = q.GaugeAction(5.5, 0.0)
test_hmc(total_site, ga)
size_node_list = [[1, 1, 1, 1], [1, 1, 1, 2], [1, 1, 1, 4], [1, 2, 2, 2],
[2, 2, 2, 2], [2, 2, 2, 4]]
q.begin(sys.argv, size_node_list)
# q.show_machine()
q.qremove_all_info("results")
main()
q.timer_display()
q.end()