fn apply_dipole_conv( fft_ws: &mut Fft3dWorkspaceF32, x: &[f32], d_kernel: &[f32], out: &mut [f32], complex_buf: &mut [Complex32], )
Apply dipole convolution: out = real(ifft(D * fft(x)))