Corrected mistake in MD integration for omf2

This commit is contained in:
Alberto Ramos 2021-12-27 19:20:26 +01:00
parent aeadaeee86
commit 6f6773d2ed

View file

@ -115,6 +115,9 @@ function MD!(mom, U, int::IntrScheme{NI, T}, lp::SpaceParm, gp::GaugeParm, ymws:
else else
mom .= mom .+ (int.r[k]*ee) .* ymws.frc1 mom .= mom .+ (int.r[k]*ee) .* ymws.frc1
end end
if k == NI
off = -1
end
k += off k += off
end end
end end