Bug corrected

This commit is contained in:
Alberto Ramos 2021-11-14 21:21:20 +01:00
parent 6dd69f192a
commit 6fbc516a6a

View file

@ -429,18 +429,17 @@ function krnl_field_tensor!(frc1, frc2, U::AbstractArray{T}, Ubnd, ipl1, ipl2, z
if ru1 == r if ru1 == r
gt1 = Ush[bu1,2] gt1 = Ush[bu1,2]
else else
if SFBC && (it == lp.iL[end])
gt1 = Ubnd[id2]
else
gt1 = U[bu1,id2,ru1] gt1 = U[bu1,id2,ru1]
end end
end
if ru2 == r if ru2 == r
gt2 = Ush[bu2,1] gt2 = Ush[bu2,1]
else
if SFBC && (it == lp.iL[end])
gt2 = Ubnd
else else
gt2 = U[bu2,id1,ru2] gt2 = U[bu2,id1,ru2]
end end
end
l1 = gt1/gt2 l1 = gt1/gt2
l2 = Ush[b,2]\Ush[b,1] l2 = Ush[b,2]\Ush[b,1]