new way to compute quaternion Shat_E ala RRA 2000