dpar.rep modified in randomize

src/Dirac/Dirac.jl

@@ -324,22 +324,20 @@ end
 
 
 
-function pfrandomize!(f,lp::SpaceParm,dpar::DiracParam,t::Int64=0) #DiracParam to apply dirac operator later?
+function pfrandomize!(f,lp::SpaceParm,dpar::DiracParam{T,SU3fund},t::Int64=0) where {T} 
 
-    if dpar.rep == SU3fund && lp.ndim == 4
         @timeit "Randomize pseudofermion field" begin
             p = ntuple(i->CUDA.randn(Float64, lp.bsz, 3, lp.rsz,2),4) # complex generation not suported for Julia 1.5.4
             CUDA.@sync begin
                 CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_assign_pf!(f,p,lp,t)
             end
         end
-    end
 
     return nothing
 
 end
 
-function krnl_assign_pf!(f::AbstractArray{T}, p ,lp::SpaceParm, t::Int64) where {T} #only valid for SU3fund for now. Check performance and maybe change it for the tuple gen in the krnl
+function krnl_assign_pf!(f::AbstractArray{T}, p ,lp::SpaceParm, t::Int64) where {T} 
 
     @inbounds begin
         b = Int64(CUDA.threadIdx().x)