pub fn rq_decomposition<R: RealField + Copy>( orig: Matrix3<R>, ) -> Result<(UnitQuaternion<R>, Matrix3<R>)>
perform RQ decomposition and return results as right-handed quaternion and intrinsics matrix