2nd bug

src/Dirac/Dirac.jl

@@ -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, csw, 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)
                     end
             end
         end