During a leap-frog integration step all λ i -coordinates are first propagated without constraints to their unconstrained new values λ u i , which do not fulfill the constraints in Equation . To obtain the constrained λ c i values, we first connect the constrained and unconstrained λ i values using the definition of σ k (λ):