There are currently 2 types implementations for pseudospectral ODE solves: real and FFT space. solve_ODE and the looped solve both uses FFT space, but the leading order will tentatively use real space. Can this be cleaned up in the future? Keep it in mind.
`