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

Set to zero explicitly SF clover fields

Browse source
This commit is contained in:
Alberto Ramos 2021-11-14 22:12:04 +01:00
parent 6fbc516a6a
commit 500305f76b
1 changed files with 4 additions and 0 deletions
Download patch file Download diff file Expand all files Collapse all files

4
src/YM/YMflow.jl
View file

@ -446,6 +446,8 @@ function krnl_field_tensor!(frc1, frc2, U::AbstractArray{T}, Ubnd, ipl1, ipl2, z
if SFBC && (it == lp.iL[end]) if SFBC && (it == lp.iL[end])
frc1[b,1,r] = projalg(Ush[b,1]*l1/Ush[b,2]) frc1[b,1,r] = projalg(Ush[b,1]*l1/Ush[b,2])
frc1[bu1,2,ru1] = zero(frc1[bu1,2,ru1])
frc1[bd,3,rd] = zero(frc1[bu1,2,ru1])
frc1[bu2,4,ru2] = projalg(l2*l1) frc1[bu2,4,ru2] = projalg(l2*l1)
else else
frc1[b,1,r] = projalg(ztw1, Ush[b,1]*l1/Ush[b,2]) frc1[b,1,r] = projalg(ztw1, Ush[b,1]*l1/Ush[b,2])
@ -485,6 +487,8 @@ function krnl_field_tensor!(frc1, frc2, U::AbstractArray{T}, Ubnd, ipl1, ipl2, z
if SFBC && (it == lp.iL[end]) if SFBC && (it == lp.iL[end])
frc2[b,1,r] = projalg(Ush[b,1]*l1/Ush[b,2]) frc2[b,1,r] = projalg(Ush[b,1]*l1/Ush[b,2])
frc1[bu1,2,ru1] = zero(frc1[bu1,2,ru1])
frc1[bd,3,rd] = zero(frc1[bu1,2,ru1])
frc2[bu2,4,ru2] = projalg(l2*l1) frc2[bu2,4,ru2] = projalg(l2*l1)
else else
frc2[b,1,r] = projalg(ztw1, Ush[b,1]*l1/Ush[b,2]) frc2[b,1,r] = projalg(ztw1, Ush[b,1]*l1/Ush[b,2])