def METHOD_NAME(q1: qt.quaternion, q2: qt.quaternion) -> float:
    r"""Computes the angular distance between two quaternions
    :return: The angular distance between q1 and q2 in radians
    """
    q1_inv = np.conjugate(q1)
    dq = q1_inv * q2
    return 2 * np.arctan2(np.linalg.norm(dq.imag), np.abs(dq.real))