From 1c1a54348bbc000ca5fa1f004d57e0c134c5588a Mon Sep 17 00:00:00 2001
From: "Fernando P.Panadero" <fernando.p@csic.es>
Date: Mon, 25 Nov 2024 11:44:34 +0100
Subject: [PATCH] Bug in plaq_pln

---
 src/YM/YMflow.jl | 14 +++++++-------
 1 file changed, 7 insertions(+), 7 deletions(-)

diff --git a/src/YM/YMflow.jl b/src/YM/YMflow.jl
index 645da0f..335c9d2 100644
--- a/src/YM/YMflow.jl
+++ b/src/YM/YMflow.jl
@@ -316,30 +316,30 @@ Eoft_plaq(U, gp::GaugeParm{T,G,NN}, lp::SpaceParm{N,M,B,D}, ymws::YMworkspace) w
 
 
 function krnl_plaq_pln!(plx, U::AbstractArray{T}, Ubnd, ztw, ipl, lp::SpaceParm{N,M,B,D}) where {T,N,M,B,D}
-    
+
     @inbounds begin
         b = Int64(CUDA.threadIdx().x)
         r = Int64(CUDA.blockIdx().x)
         I = point_coord((b,r), lp)
-        
+
         id1, id2 = lp.plidx[ipl]
         SFBC = ((B == BC_SF_AFWB) || (B == BC_SF_ORBI)) && (id1 == N)
         TWP  = ((I[id1]==1)&&(I[id2]==1))
-        
+
         bu1, ru1 = up((b, r), id1, lp)
         bu2, ru2 = up((b, r), id2, lp)
-        
-        if SFBC && (ru1 != r)
+
+        if SFBC && (point_time((b,r),lp) == lp.iL[end])
             gt = Ubnd[id2]
         else
             gt = U[bu1,id2,ru1]
         end
-        
+
         if TWP
             plx[I] = ztw*tr(U[b,id1,r]*gt / (U[b,id2,r]*U[bu2,id1,ru2]))
         else
             plx[I] = tr(U[b,id1,r]*gt / (U[b,id2,r]*U[bu2,id1,ru2]))
-        end            
+        end
     end
     return nothing
 end