Batched processing for complex FFTs

This PR enables batched processing for complex FFTs, just like the 'axis' keyword for NumPy's FFT. It also adds a simple function 'fftspace' that returns the sample frequencies of a FFT.