Fmunu added in secon index

parent 3646f825
...@@ -291,7 +291,7 @@ function krnl_Dw_impr!(so, F1, F2, Gid1, Gid2, si, csw, lp::SpaceParm{4,6,B,D}) ...@@ -291,7 +291,7 @@ function krnl_Dw_impr!(so, F1, F2, Gid1, Gid2, si, csw, lp::SpaceParm{4,6,B,D})
@inbounds begin @inbounds begin
so[b,r] += 0.25*csw*im*( F1[b,r]*dmul(Gid1,si[b,r]) - F2[b,r]*dmul(Gid2,si[b,r])) so[b,r] += 0.25*csw*im*( (F1[b,1,r]+F1[b,2,r]+F1[b,3,r]+F1[b,4,r])*dmul(Gid1,si[b,r]) - (F2[b,1,r]+F2[b,2,r]+F2[b,3,r]+F2[b,4,r])*dmul(Gid2,si[b,r]))
end end
...@@ -304,7 +304,7 @@ function krnl_g5Dw_impr!(so, F1, F2, Gid1, Gid2, si, csw, lp::SpaceParm{4,6,B,D} ...@@ -304,7 +304,7 @@ function krnl_g5Dw_impr!(so, F1, F2, Gid1, Gid2, si, csw, lp::SpaceParm{4,6,B,D}
@inbounds begin @inbounds begin
so[b,r] += 0.25*csw*im*dmul(Gamma{5},( F1[b,r]*dmul(Gid1,si[b,r]) - F2[b,r]*dmul(Gid2,si[b,r]))) so[b,r] += 0.25*csw*im*dmul(Gamma{5},( (F1[b,1,r]+F1[b,2,r]+F1[b,3,r]+F1[b,4,r])*dmul(Gid1,si[b,r]) - (F2[b,1,r]+F2[b,2,r]+F2[b,3,r]+F2[b,4,r])*dmul(Gid2,si[b,r])))
end 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