pub fn hadamard_simd( amplitudes: &mut [Complex64], qubit: usize, _num_qubits: usize, )
Apply a Hadamard gate using SIMD operations
This applies H = (1/√2) * [[1, 1], [1, -1]] to the specified qubit.