lattice/latticegpu.jl
1
0
Fork 0
mirror of https://igit.ific.uv.es/alramos/latticegpu.jl.git synced 2025-06-29 21:39:27 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

changed time-boundary of Eoft_clover and Eoft_plaq to match openQCD

Browse source
This commit is contained in:
Nicolas Lang 2025-05-20 16:58:19 +02:00
parent 349ff2405f
commit 58685f63a2
1 changed files with 10 additions and 14 deletions
Download patch file Download diff file Expand all files Collapse all files

20
src/YM/YMflow.jl
View file

@ -430,28 +430,24 @@ function Eoft_clover(Eslc, U, gp::GaugeParm, lp::SpaceParm{4,M,B,D}, ymws::YMwor
end
Etmp .= reshape(Array(CUDA.reduce(+, ymws.rm;dims=tp)),lp.iL[end])/V3
for it in 1:lp.iL[end]
if (B == BC_OPEN) && (it == 1 || it == lp.iL[end])
Eslc[it,ipl1] = 0.0
else
Eslc[it,ipl1] = Etmp[it]/8
end
## tentative (weight 1/2 on spatial planes on boundary)
#if (B == BC_OPEN) && (ipl1 >= 4)
# Eslc[end,ipl1] = Etmp[end]/16
# Eslc[1,ipl1] = Etmp[1]/16
#end
end
CUDA.@sync begin
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_add_et!(ymws.rm, ymws.frc2, lp)
end
Etmp .= reshape(Array(CUDA.reduce(+, ymws.rm;dims=tp)),lp.iL[end])/V3
for it in 1:lp.iL[end]
if (B == BC_OPEN) && (it == 1 || it == lp.iL[end])
Eslc[it,ipl2] = 0.0
else
Eslc[it,ipl2] = Etmp[it]/8
end
## tentative (weight 1/2 on spatial planes on boundary)
#if (B == BC_OPEN) && (ipl2 >= 4)
# Eslc[end,ipl2] = Etmp[end]/16
# Eslc[1,ipl2] = Etmp[1]/16
#end
end
return nothing
end