tsa
power spectral density using padding to length n done by fft
currently returns two-sided according to fft frequencies, use first half
ArmaFft.spddirect