“…Transformations applied to systems with holonomic constraints usually aim at decoupling the dynamics into a "tangent" (independent of the contact force) and "normal" (contact force dependent) dynamics, using various annihilators of ∇h(q) [15,16,27,29]. Systems with a singular mass matrix are descriptor variable systems, or implicit systems in the quasi-linear form G(ẋ, x, t, λ) = E(x)ẋ + H (x, t, λ), where…”