コード例 #1
0
########################################################

# staggered parameters
p = {
    "mass": 0,
    "mu5": 1,
    "hop": 0,
    "boundary_phases": [1.0, 1.0, 1.0, 1.0],
}

# grid (each dimension must be at least 4 to get correct sum rule)
L = [8, 4, 4, 4]
grid_dp = g.grid(L, g.double)
grid_sp = g.grid(L, g.single)

# SU(2) fundamental
U = g.qcd.gauge.random(grid_sp,
                       g.random("test"),
                       otype=g.ot_matrix_su2_fundamental())
ev = run_test(U)

# SU(2) adjoint
U = g.qcd.gauge.random(grid_sp,
                       g.random("test"),
                       otype=g.ot_matrix_su2_adjoint())
run_test(U)

# SU(3) fundamental
U = g.qcd.gauge.random(grid_sp, g.random("test"))
run_test(U)
コード例 #2
0
ファイル: staggered_complex.py プロジェクト: wettig/gpt
    return 1


#################################################################
# test sum rules for different gauge groups and representations #
#################################################################

# staggered parameters
p = {
    "mass": .897,
    "hop": 1,
    "mu5": complex(1.23, .537),
    "boundary_phases": [1.0, 1.0, 1.0, 1.0],
}

# grid (each dimension must be at least 4 to get correct sum rule)
L = [8, 4, 4, 4]
grid_dp = g.grid(L, g.double)

# SU(2) fundamental
U = g.qcd.gauge.random(grid_dp, g.random("test"), otype=g.ot_matrix_su2_fundamental())
test_sumrule(U, p)

# SU(2) adjoint
U = g.qcd.gauge.random(grid_dp, g.random("test"), otype=g.ot_matrix_su2_adjoint())
test_sumrule(U, p)

# SU(3) fundamental
U = g.qcd.gauge.random(grid_dp, g.random("test"))
test_sumrule(U, p)