pub fn expectation_z_simd( amplitudes: &[Complex64], qubit: usize, _num_qubits: usize, ) -> f64
Compute expectation value of a Pauli Z operator
This computes ⟨ψ|Z|ψ⟩ where Z is the Pauli Z operator on the given qubit.