On the Adjoint Formulation of Design Sensitivity Analysis of Multibody Dynamics. The Piecewise Adjoint Method for Design Sensitivity Analysis (DSA) of Non-centroidal Multibody Dynamics with Euler Parameters for Orientation