Update angular velocity in Particle::correction_step
Store integration parameters (including gamma) in Contact class
Re-calculate Contact::_dist, _delta and _a on demand