on line 215, subset residIneq vector to matcch how it's defined below equation (9)
fixed dimension mismatch issue when finding next rho for inequality constraint event (escaping boundary)