Add a norm keyword arg to fft functions

This fixes #2142.
division is expensive `*= 1/ sqrt(n)` is faster and the rounding difference should be negligible for fft