pub fn cross_matrix(v: Vec3) -> [[f32; 3]; 3]
Cross product matrix (skew symmetric) for angular dynamics