fn rodrigues_rotate(v: [f64; 3], k: [f64; 3], theta: f64) -> [f64; 3]
Rodrigues’ rotation formula: rotate vector v around unit axis k by angle θ.