3rd bug: match gamma/F not done yet

parent 75cdb213
......@@ -78,7 +78,7 @@ function Dw!(so, U, si, dpar::DiracParam, lp::SpaceParm{4,6,B,D}, ymws = "ymws",
z2 = exp(2im * pi * lp.ntw[lp.npls-i]/gp.ng)
CUDA.@sync begin
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_field_tensor!(ymws.frc1, ymws.frc2, U, gp.Ubnd, i, lp.npls-i, z1, z2, lp)
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_Dw_impr!(so, ymws.frc1, ymws.frc2, i, lp.npls-i, si, dpar.csw, lp)
CUDA.@cuda threads=lp.bsz blocks=lp.rsz krnl_Dw_impr!(so, ymws.frc1, ymws.frc2, Gamma{i}, Gamma{lp.npls-i}, si, dpar.csw, lp)
end
end
end
......@@ -229,13 +229,13 @@ function krnl_g5Dw!(so, U, si, m0, th, lp::SpaceParm{4,6,B,D}) where {B,D}
end
function krnl_Dw_impr!(so, F1, F2, id1, id2, si, csw, lp::SpaceParm{4,6,B,D}) where {B,D}
function krnl_Dw_impr!(so, F1, F2, Gid1, Gid2, si, csw, lp::SpaceParm{4,6,B,D}) where {B,D}
b = Int64(CUDA.threadIdx().x); r = Int64(CUDA.blockIdx().x)
@inbounds begin
# asignation between /mu/nu index and sigma index is NOT correct
so[b,r] += 0.25*csw*im*( F1*dmul(Gamma{id1},si[b,r]) + F2*dmul(Gamma{id2},si[b,r]))
so[b,r] += 0.25*csw*im*( F1[b,r]*dmul(Gid1,si[b,r]) + F2[b,r]*dmul(Gid2,si[b,r]))
end
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment