Function half_norm 

fn half_norm(v: &[Complex64], nx: usize) -> f64

Norm of a Hermitian-symmetric spectrum as QSM.jl sees it: over the kx = 0..=nx/2 half that its real-to-complex FFT stores.