Function ah_op 

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

A^H operator: 1-channel field φ → 2-channel [real(ifft(D·fft(φ))), φ]. Interleaved output [ch0, ch1] per voxel.