Small fix so that everything doesnt explote

src/Dirac/Diracflw.jl

@@ -30,7 +30,7 @@ function flw(U, psi, int::FlowIntr{NI,T}, ns::Int64, eps, gp::GaugeParm, dpar::D
                 ymws.mom .= int.e0[k].*ymws.mom .+ int.e1[k].*ymws.frc1
                 U .= expm.(U, ymws.mom, 2*eps)
 
-                Nablanabla!(dws.sAp, U, psi, dpar, dws, lp)
+                Nablanabla!(dws.sp, U, psi, dpar, dws, lp)
 
                 dws.sAp .= int.e0[k].*dws.sAp .+ int.e1[k].*dws.sp
                 psi .= psi + 2*eps*dws.sAp