#!/usr/bin/env python3
#
# Authors: Christoph Lehner 2020
#
# Production code to generate fine-grid basis vectorscoarse-grid eigenvectors using existing
#
import gpt as g
# parameters
fn = g.default.get("--params", "params.txt")
params = g.params(fn, verbose=True)
# load configuration
U = params["config"]
# matrix to use
fmatrix = params["fmatrix"](U)
op = params["op"](fmatrix)
grid = op.grid[0]
# implicitly restarted lanczos
irl = params["method_evec"]
# run
start = g.vspincolor(grid)
start[:] = g.vspincolor([[1, 1, 1], [1, 1, 1], [1, 1, 1], [1, 1, 1]])
start.checkerboard(g.odd) # traditionally, calculate odd-site vectors
try:
basis, ev = g.load("checkpoint", grids=grid)
except g.LoadError:
# Authors: Christoph Lehner 2020
# Lorenzo Barca 2020
import gpt as g
import numpy as np
from gpt.qcd.hadron import spinmatrices as spm
from gpt.qcd.hadron.twopointcontraction.baryonoctet import \
contract_proton, contract_xi_zero, contract_sigma_plus, contract_lambda_naive, contract_lambda
from gpt.qcd.hadron.twopointcontraction.baryondecuplet import \
contract_sigma_plus_star, contract_delta_plus, contract_xi_zero_star, contract_omega
# show available memory
g.mem_report()
parameter_file = g.default.get("--params", "params.txt")
params = g.params(parameter_file, verbose=True)
# load configuration
rng = g.random("test")
U = g.qcd.gauge.random(g.grid([8, 8, 8, 16], g.double), rng)
g.mem_report()
# use the gauge configuration grid
grid = U[0].grid
vol = np.array(grid.fdimensions)
L, Nt = vol[0], vol[-1]
# smear gauge links
params_link_smear = params["APE"]
g.message("Start APE-link smearing...")
U_ape = U