Sorry, I've found - we really need to do all in the same coordinate system.
But now, I have some other question - probably, more interesting. The torque acting on the rigid body (expressed in the body-fixed frame) is given by: t = -0.5*S(q)*(dU/dq), gde U - potential, q - is a quaternion, ans...