Function a_op 

fn a_op(sus: &[[f64; 2]], dk: &[f64], grid: &Grid) -> Vec<f64>

A operator: 2-channel susceptibility → 1-channel field real(ifft(D·fft(χ33) + fft(χ13/23))). Input is interleaved [ch0, ch1].