From 50f07f40e603202207fdf66947c59da61b6d45d1 Mon Sep 17 00:00:00 2001 From: "Fernando P. Panadero" Date: Tue, 21 Nov 2023 14:37:08 +0100 Subject: [PATCH] Constructor for dws fixed --- src/Dirac/Dirac.jl | 2 +- test/dirac/test_fp_fa.jl | 4 ++-- test/dirac/test_solver_plw.jl | 2 +- test/dirac/test_solver_rand.jl | 2 +- 4 files changed, 5 insertions(+), 5 deletions(-) diff --git a/src/Dirac/Dirac.jl b/src/Dirac/Dirac.jl index 832b3e9..8fd19af 100644 --- a/src/Dirac/Dirac.jl +++ b/src/Dirac/Dirac.jl @@ -72,7 +72,7 @@ struct DiracWorkspace{T} end end - return new{T}(sr,sp,sAp,st,csw,cs) + return new{T}(sr,sp,sAp,st,csw) end end diff --git a/test/dirac/test_fp_fa.jl b/test/dirac/test_fp_fa.jl index ee931ae..6610bb4 100644 --- a/test/dirac/test_fp_fa.jl +++ b/test/dirac/test_fp_fa.jl @@ -14,7 +14,7 @@ lp = SpaceParm{4}(size,(4,4,4,4),1,(0,0,0,0,0,0)); exptheta = exp.(im.*theta./lp.iL); dpar = DiracParam{Float64}(SU3fund,m,0.0,exptheta,1.0); -dws = DiracWorkspace(SU3fund{Float64},Float64,lp); +dws = DiracWorkspace(SU3fund,Float64,lp); U = fill!(vector_field(SU3{Float64},lp),one(SU3{Float64})); psi = scalar_field(Spinor{4,SU3fund{Float64}},lp); @@ -66,7 +66,7 @@ function fA_test(;theta = (0.5,0.7,1.0,0.0), m = 1.3, size = (8,8,8,16),prec = 1 exptheta = exp.(im.*theta./lp.iL); dpar = DiracParam{Float64}(SU3fund,m,0.0,exptheta,1.0); - dws = DiracWorkspace(SU3fund{Float64},Float64,lp); + dws = DiracWorkspace(SU3fund,Float64,lp); U = fill!(vector_field(SU3{Float64},lp),one(SU3{Float64})); psi = scalar_field(Spinor{4,SU3fund{Float64}},lp); diff --git a/test/dirac/test_solver_plw.jl b/test/dirac/test_solver_plw.jl index 14100d0..195a1a3 100644 --- a/test/dirac/test_solver_plw.jl +++ b/test/dirac/test_solver_plw.jl @@ -8,7 +8,7 @@ function Dwpw_test(;p=0,s=1,c=1) lp = SpaceParm{4}((16,16,16,16), (4,4,4,4), 0, (0,0,0,0,0,0)) gp = GaugeParm{Float64}(SU3{Float64}, 6.0, 1.0) dpar = DiracParam{Float64}(SU3fund,1.3,0.0,(1.0,1.0,1.0,1.0),0.0) -dws = DiracWorkspace(SU3fund{Float64},Float64,lp); +dws = DiracWorkspace(SU3fund,Float64,lp); p==0 ? p = Int.(round.(lp.iL.*rand(4),RoundUp)) : nothing U = fill!(vector_field(SU3{Float64},lp),one(SU3{Float64})) diff --git a/test/dirac/test_solver_rand.jl b/test/dirac/test_solver_rand.jl index 91a477d..c06de0a 100644 --- a/test/dirac/test_solver_rand.jl +++ b/test/dirac/test_solver_rand.jl @@ -10,7 +10,7 @@ using CUDA, LatticeGPU, TimerOutputs gp = GaugeParm{Float64}(SU3{Float64}, 6.0, 1.0) ymws = YMworkspace(SU3, Float64, lp) dpar = DiracParam{Float64}(SU3fund,2.3,0.0,(1.0,1.0,1.0,1.0),0.0) - dws = DiracWorkspace(SU3fund{Float64},Float64,lp); + dws = DiracWorkspace(SU3fund,Float64,lp); randomize!(ymws.mom, lp, ymws) U = exp.(ymws.mom)