Function apply_ry_simd 

pub fn apply_ry_simd(
    state: &mut Vec<Complex64>,
    theta: f64,
    target: usize,
    n_qubits: usize,
)

Apply the RY(theta) rotation gate to qubit target using SIMD.

RY(θ) = [[cos(θ/2), −sin(θ/2)], [sin(θ/2), cos(θ/2)]]