From 500305f76b3f5e7f23bb9a5e6dbc7094ec736cbf Mon Sep 17 00:00:00 2001 From: Alberto Ramos Date: Sun, 14 Nov 2021 22:12:04 +0100 Subject: [PATCH] Set to zero explicitly SF clover fields --- src/YM/YMflow.jl | 4 ++++ 1 file changed, 4 insertions(+) diff --git a/src/YM/YMflow.jl b/src/YM/YMflow.jl index 5286c5f..3287e29 100644 --- a/src/YM/YMflow.jl +++ b/src/YM/YMflow.jl @@ -446,6 +446,8 @@ function krnl_field_tensor!(frc1, frc2, U::AbstractArray{T}, Ubnd, ipl1, ipl2, z if SFBC && (it == lp.iL[end]) 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) else 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]) 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) else frc2[b,1,r] = projalg(ztw1, Ush[b,1]*l1/Ush[b,2])