def zmobius(U, params):
params = copy.deepcopy(params) # save current parameters
assert "Ls" not in params
params["Ls"] = len(params["omega"])
return fine_operator("zmobius",
U,
params,
otype=gpt.ot_vector_spin_color(4, 3))
def wilson_clover(U, params):
params = copy.deepcopy(params) # save current parameters
if "kappa" in params:
assert "mass" not in params
params["mass"] = 1.0 / params["kappa"] / 2.0 - 4.0
del params["kappa"]
return fine_operator("wilson_clover",
U,
params,
otype=gpt.ot_vector_spin_color(4, 3))
def wilson_clover(U, params):
params = copy.deepcopy(params) # save current parameters
if "kappa" in params:
assert "mass" not in params
params["mass"] = 1.0 / params["kappa"] / 2.0 - 4.0
del params["kappa"]
if "cF" not in params:
params["cF"] = 1.0 # default to 1.0 (i.e., doing nothing)
if "use_legacy" not in params:
params["use_legacy"] = False # default to new, faster implementation
if params["use_legacy"]:
assert params["boundary_phases"][
-1] != 0.0 # only new op supports open bc
if params["boundary_phases"][-1] != 0.0:
assert params["cF"] == 1.0 # forbid usage of cF without open bc
return fine_operator("wilson_clover",
U,
params,
otype=gpt.ot_vector_spin_color(4, 3))
def mobius(U, params):
params = copy.deepcopy(params) # save current parameters
return fine_operator("mobius",
U,
params,
otype=gpt.ot_vector_spin_color(4, 3))