######################################################## # 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)
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)