if w.shape[0] == 1 and w.shape[1] > 1 or w.shape[0] > 1 and w.shape[1] == 1: G = pyfftw.interfaces.numpy_fft.fft(Lo, cfslength, planner_effort='FFTW_ESTIMATE ...

Consider a particle moving in a harmonic oscillator potential, with the associated dimensionless Hamiltonian given by: $H^0 = \frac{p^2}{2}+\frac{1}{2} x^2$ with the ...