diff --git a/src/Dirac/Diracfields.jl b/src/Dirac/Diracfields.jl index 5559ca3..3210460 100644 --- a/src/Dirac/Diracfields.jl +++ b/src/Dirac/Diracfields.jl @@ -151,15 +151,17 @@ function krnl_assign_pf_su3!(f::AbstractArray, p , lp::SpaceParm, t::Int64) b = Int64(CUDA.threadIdx().x) r = Int64(CUDA.blockIdx().x) - if t == 0 + if t == 0 f[b,r] = Spinor(map(x->SU3fund(x[b,1,r,1] + im* x[b,1,r,2], - x[b,2,r,1] + im* x[b,2,r,2], - x[b,3,r,1] + im* x[b,3,r,2]),p)) - elseif point_time((b,r),lp) == t + x[b,2,r,1] + im* x[b,2,r,2], + x[b,3,r,1] + im* x[b,3,r,2]),p)) + elseif point_time((b,r),lp) == t f[b,r] = Spinor(map(x->SU3fund(x[b,1,r,1] + im* x[b,1,r,2], - x[b,2,r,1] + im* x[b,2,r,2], - x[b,3,r,1] + im* x[b,3,r,2]),p)) - end + x[b,2,r,1] + im* x[b,2,r,2], + x[b,3,r,1] + im* x[b,3,r,2]),p)) + else + f[b,r] = 0.0*f[b,r] + end end @@ -197,13 +199,15 @@ function krnl_assign_pf_su2!(f::AbstractArray, p , lp::SpaceParm, t::Int64) b = Int64(CUDA.threadIdx().x) r = Int64(CUDA.blockIdx().x) - if t == 0 + if t == 0 f[b,r] = Spinor(map(x->SU2fund(x[b,1,r,1] + im* x[b,1,r,2], - x[b,2,r,1] + im* x[b,2,r,2]),p)) - elseif point_time((b,r),lp) == t + x[b,2,r,1] + im* x[b,2,r,2]),p)) + elseif point_time((b,r),lp) == t f[b,r] = Spinor(map(x->SU2fund(x[b,1,r,1] + im* x[b,1,r,2], - x[b,2,r,1] + im* x[b,2,r,2]),p)) - end + x[b,2,r,1] + im* x[b,2,r,2]),p)) + else + f[b,r] = 0.0*f[b,r] + end end diff --git a/src/YM/YMflow.jl b/src/YM/YMflow.jl index 42ed545..645da0f 100644 --- a/src/YM/YMflow.jl +++ b/src/YM/YMflow.jl @@ -93,7 +93,7 @@ function Base.show(io::IO, int::FlowIntr{N,T}) where {N,T} if N == 0 println(io, " * Euler schem3") elseif N == 1 - println(io, " * One stage scheme. Coefficients3") + println(io, " * One stage scheme. Coefficients") println(io, " stg 1: ", int.e0[1], " ", int.e1[1]) elseif N == 2 println(io, " * Two stage scheme. Coefficients:") diff --git a/src/YM/YMsf.jl b/src/YM/YMsf.jl index abd93d9..cd00ac7 100644 --- a/src/YM/YMsf.jl +++ b/src/YM/YMsf.jl @@ -92,7 +92,7 @@ end """ function setbndfield(U, phi, lp::SpaceParm) -Sets abelian boundary fields with phases `phi[1]` and `phi[2]` to the configuration `U` at time salice ``x_0=0``. +Sets abelian boundary fields with phases `phi[1]` and `phi[2]` to the configuration `U` at time slice ``x_0=0``. """ function setbndfield(U, phi, lp::SpaceParm{N,M,B,D}) where {N,M,B,D}