Wrote a script for testing the `Fourier.phase_shifted_coefficients`
function. It crashes because I apparently don't understand `numpy` array
broadcasting. Should work after some minor debugging.