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

Corrected bug for improved SF simulations

Browse source
This commit is contained in:
Alberto Ramos 2021-12-12 01:28:45 +01:00
parent 46a1c996e5
commit aeadaeee86
1 changed files with 10 additions and 2 deletions
Download patch file Download diff file Expand all files Collapse all files

12
src/YM/YMact.jl
View file

@ -41,7 +41,11 @@ function krnl_impr!(plx, U::AbstractArray{T}, c0, c1, Ubnd::NTuple{NB,T}, cG, zt
if r2 == r if r2 == r
gb = Ush[b2,2] gb = Ush[b2,2]
else else
gb = U[b2,id2,r2] if SFBC && (it == lp.iL[end]-1)
gb = Ubnd[id2]
else
gb = U[b2,id2,r2]
end
end end
(b2, r2) = up((b1,r1), id2, lp) (b2, r2) = up((b1,r1), id2, lp)
@ -279,7 +283,11 @@ function krnl_force_impr_pln!(frc1, frc2, U::AbstractArray{T}, c0, c1, Ubnd, cG,
if r2 == r if r2 == r
gb = Ush[b2,2] gb = Ush[b2,2]
else else
gb = U[b2,id2,r2] if SFBC && (it == lp.iL[end]-1)
gb = Ubnd[id2]
else
gb = U[b2,id2,r2]
end
end end
(b2, r2) = up((b1,r1), id2, lp) (b2, r2) = up((b1,r1), id2, lp)