pub fn apply_phase_simd(amplitudes: &mut [Complex64], theta: f64)
Apply a phase rotation to a quantum state vector using SIMD when available
This function applies the phase rotation e^(i*theta) to each amplitude.