lattice/latticegpu.jl
1
0
Fork 0
mirror of https://igit.ific.uv.es/alramos/latticegpu.jl.git synced 2025-05-14 19:23:42 +02:00
Code Issues Projects Releases Packages Wiki Activity Actions

Added zeuthen term for SF AFW-B boundary conditions

Browse source
This commit is contained in:
Alberto Ramos 2021-10-25 19:38:07 +02:00
parent 5336189ad5
commit e9670762f8
1 changed files with 16 additions and 3 deletions
Download patch file Download diff file Expand all files Collapse all files

17
src/YM/YMflow.jl
View file

@ -28,6 +28,9 @@ end
function krnl_add_zth!(frc, frc2::AbstractArray{TA}, U::AbstractArray{TG}, lp::SpaceParm{N,M,B,D}) where {TA,TG,N,M,B,D} function krnl_add_zth!(frc, frc2::AbstractArray{TA}, U::AbstractArray{TG}, lp::SpaceParm{N,M,B,D}) where {TA,TG,N,M,B,D}
b, r = CUDA.threadIdx().x, CUDA.blockIdx().x b, r = CUDA.threadIdx().x, CUDA.blockIdx().x
it = point_time((b, r), lp)
SFBC = ((B == BC_SF_AFWB) || (B == BC_SF_ORBI) )
Ush = @cuStaticSharedMem(TG, D) Ush = @cuStaticSharedMem(TG, D)
Fsh = @cuStaticSharedMem(TA, D) Fsh = @cuStaticSharedMem(TA, D)
@ -53,8 +56,18 @@ function krnl_add_zth!(frc, frc2::AbstractArray{TA}, U::AbstractArray{TG}, lp::S
Ud = U[bd,id,rd] Ud = U[bd,id,rd]
end end
frc2[b,id,r] = (5/6)*Fsh[b] + (1/6)*(projalg(Ud\Y*Ud) + if SFBC
projalg(Ush[b]*X/Ush[b])) if (it > 1) && (it < lp.iL[end])
frc2[b,id,r] = (5/6)*Fsh[b] + (1/6)*(projalg(Ud\Y*Ud) +
projalg(Ush[b]*X/Ush[b]))
elseif (it == lp.iL[end]) && (id < N)
frc2[b,id,r] = (5/6)*Fsh[b] + (1/6)*(projalg(Ud\Y*Ud) +
projalg(Ush[b]*X/Ush[b]))
end
else
frc2[b,id,r] = (5/6)*Fsh[b] + (1/6)*(projalg(Ud\Y*Ud) +
projalg(Ush[b]*X/Ush[b]))
end
end end
return nothing return nothing