From 58685f63a243c06ab0a3b7209ccbc3225b4c448c Mon Sep 17 00:00:00 2001 From: Nicolas Lang Date: Tue, 20 May 2025 16:58:19 +0200 Subject: [PATCH] changed time-boundary of Eoft_clover and Eoft_plaq to match openQCD --- src/YM/YMflow.jl | 24 ++++++++++-------------- 1 file changed, 10 insertions(+), 14 deletions(-) diff --git a/src/YM/YMflow.jl b/src/YM/YMflow.jl index 8004653..92c8e2c 100644 --- a/src/YM/YMflow.jl +++ b/src/YM/YMflow.jl @@ -430,29 +430,25 @@ 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] - Eslc[it,ipl1] = Etmp[it]/8 + if (B == BC_OPEN) && (it == 1 || it == lp.iL[end]) + Eslc[it,ipl1] = 0.0 + else + Eslc[it,ipl1] = Etmp[it]/8 + end 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 - 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] - Eslc[it,ipl2] = Etmp[it]/8 + if (B == BC_OPEN) && (it == 1 || it == lp.iL[end]) + Eslc[it,ipl2] = 0.0 + else + Eslc[it,ipl2] = Etmp[it]/8 + end 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 - return nothing end