Optimization:
math::cross(mom, mom_dt) = math::cross( mom[ "now" - 1 ] , mom [ "now" ] )
when assuming we use the average momentum instead of a one sided version.