Function quaternion_derivative 

pub fn quaternion_derivative(q: DQuat, omega_body_rad_s: DVec3) -> DQuat

The time derivative q̇ = ½ q ⊗ (0, ω_B) of the attitude quaternion, as a (non-unit) quaternion, for a body angular velocity omega_body_rad_s resolved in the body frame.