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

pfrandomize modification and minor typos.

Browse source
This commit is contained in:
Fernando P.Panadero 2024-09-09 11:21:41 +02:00
parent 70829c8a20
commit 5fc81739f8
3 changed files with 18 additions and 14 deletions
Download patch file Download diff file Expand all files Collapse all files

28
src/Dirac/Diracfields.jl
View file

@ -151,15 +151,17 @@ function krnl_assign_pf_su3!(f::AbstractArray, p , lp::SpaceParm, t::Int64)
b = Int64(CUDA.threadIdx().x) b = Int64(CUDA.threadIdx().x)
r = Int64(CUDA.blockIdx().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], 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,2,r,1] + im* x[b,2,r,2],
x[b,3,r,1] + im* x[b,3,r,2]),p)) x[b,3,r,1] + im* x[b,3,r,2]),p))
elseif point_time((b,r),lp) == t 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], 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,2,r,1] + im* x[b,2,r,2],
x[b,3,r,1] + im* x[b,3,r,2]),p)) x[b,3,r,1] + im* x[b,3,r,2]),p))
end else
f[b,r] = 0.0*f[b,r]
end
end end
@ -197,13 +199,15 @@ function krnl_assign_pf_su2!(f::AbstractArray, p , lp::SpaceParm, t::Int64)
b = Int64(CUDA.threadIdx().x) b = Int64(CUDA.threadIdx().x)
r = Int64(CUDA.blockIdx().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], 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)) x[b,2,r,1] + im* x[b,2,r,2]),p))
elseif point_time((b,r),lp) == t 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], 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)) x[b,2,r,1] + im* x[b,2,r,2]),p))
end else
f[b,r] = 0.0*f[b,r]
end
end end

2
src/YM/YMflow.jl
View file

@ -93,7 +93,7 @@ function Base.show(io::IO, int::FlowIntr{N,T}) where {N,T}
if N == 0 if N == 0
println(io, " * Euler schem3") println(io, " * Euler schem3")
elseif N == 1 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]) println(io, " stg 1: ", int.e0[1], " ", int.e1[1])
elseif N == 2 elseif N == 2
println(io, " * Two stage scheme. Coefficients:") println(io, " * Two stage scheme. Coefficients:")

2
src/YM/YMsf.jl
View file

@ -92,7 +92,7 @@ end
""" """
function setbndfield(U, phi, lp::SpaceParm) 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} function setbndfield(U, phi, lp::SpaceParm{N,M,B,D}) where {N,M,B,D}