We were losing half the modes in the power computation.
This was because, since Nmesh became unsigned, n - Nmesh was positive
and very large instead of being negative.

Make the storage order macro an inline function instead to avoid this.
Abort if Nmesh is not even, as I'm not totally sure what FFTW expects for its output order in this case