lattice/latticegpu.jl
1
0
Fork 0
mirror of https://igit.ific.uv.es/alramos/latticegpu.jl.git synced 2025-05-14 19:23:42 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

Added some minor tests for Spinors

Browse source
This commit is contained in:
Alberto Ramos 2021-11-21 16:28:20 +01:00
parent b03e4469cc
commit 6e2cde616e
1 changed files with 26 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

26
tests/test_su3.jl
View file

@ -143,3 +143,29 @@ println(isgroup(g1))
g1 = unitarize(g1) g1 = unitarize(g1)
println(g1) println(g1)
println(isgroup(g1)) println(isgroup(g1))
println("## Fundamental")
psi = rand(SU3fund{Float64})
println(psi)
s = Spinor{4,SU3fund{Float64}}((rand(SU3fund{Float64}),rand(SU3fund{Float64}),rand(SU3fund{Float64}),rand(SU3fund{Float64})))
println(s)
for n in 1:4
t1 = pmul(Pgamma{n,1}, s)
t2 = g1\gpmul(Pgamma{n,1}, g1, s)
println("Direction (+)$n: ", LatticeGPU.norm(t1-t2))
t1 = pmul(Pgamma{n,-1}, s)
t2 = g1\gpmul(Pgamma{n,-1}, g1, s)
println("Direction (-)$n: ", LatticeGPU.norm(t1-t2))
end
for n in 1:4
t1 = pmul(Pgamma{n,1}, s)
t2 = g1*gdagpmul(Pgamma{n,1}, g1, s)
println("Direction (+)$n: ", LatticeGPU.norm(t1-t2))
t1 = pmul(Pgamma{n,-1}, s)
t2 = g1*gdagpmul(Pgamma{n,-1}, g1, s)
println("Direction (-)$n: ", LatticeGPU.norm(t1-t2))
end