pub fn cartesian_to_cylindrical(x: f64, y: f64, z: f64) -> (f64, f64, f64)
Returns (rho, phi, z) where rho is the radial distance in the xy-plane and phi is the azimuthal angle from +x.