pub fn rotation_from_axis_angle<T: Float>(omega: &V3<T>, theta: T) -> M33<T>
Compute the matrix exponential for omega theta(exponential coordinates): so(3) —> SO(3) with the Rodriguez formula