“…Therefore, there are various numerical methods developed for their solution -cf. [2,3,11,13,25,[28][29][30]. In the last decades, the numerical methods preserving the geometric invariants along the flows, such as symplectic and Lie group structures, and quantities showing that the exact solution evolves on a manifold of the dimension smaller than that of n , have been widely studied [5,8,9,27,32].…”