Discussion of the Gear-Gupta-Leimkuhler method for multibody systems including clearance and impacts