From 6e2cde616e480895169fcc1d3fc3a5fbaa9c7f0c Mon Sep 17 00:00:00 2001
From: Alberto Ramos <alberto.ramos@ific.uv.es>
Date: Sun, 21 Nov 2021 16:28:20 +0100
Subject: [PATCH] Added some minor tests for Spinors

---
 tests/test_su3.jl | 26 ++++++++++++++++++++++++++
 1 file changed, 26 insertions(+)

diff --git a/tests/test_su3.jl b/tests/test_su3.jl
index 47972c5..b11f5ee 100644
--- a/tests/test_su3.jl
+++ b/tests/test_su3.jl
@@ -143,3 +143,29 @@ println(isgroup(g1))
 g1 = unitarize(g1)
 println(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