Removed velocityHalf property and converted to component-wise calculations for the main model loop to decrease allocations, see #84, see #85