Singularity-free computation of quaternions from rotation matrices in E4 and E3

Trabajo presentado en la 7th Conference on Applied Geometric Algebras in Computer Science and Engineering (AGACSE), celebrada en Campinas (Brasil), del 23 al 27 de julio de 2018A real orthogonal matrix representing a rotation in E 4 can be decomposed into the commutative product of a left-isoclinic and a right-isoclinic rotation matrix. The double quaternion representation of rotations in E 4 follows directly from this decomposition. In this paper, it is shown how this decomposition can be performed without divisions. This avoids the common numerical issues attributed to the computation of quaternions from rotation matrices. The map from the 4×4 rotation matrices to the set of double unit quaternions is a 2-to-1 covering map. Thus, this map cannot be smoothly inverted. As a consequence, it is erroneously assumed that all inversions should necessarily contain singularities that arise in the form of quotients where the divisor can be arbitrarily small. This misconception is herein clarified. When particularized to three dimensions, it is shown how the resulting formulation outperforms, from the numerical point of view, the celebrated Shepperd’s method.This work was partially supported by the Spanish Ministry of Economy and Competitiveness through project DPI2014-57220-C2-2-P, DPI2017-88282-P, MDM-2016-0656, and the USA National Science Foundation under Grant No. 1208385.Peer reviewe