I have heard that FFT algorithms require the number of grid points to be an integer power of two. Does this restriction apply when using WavePacket?
Matlab is based on the FFTW algorithms (=Fastest Fourier Transform in the World) for which no such restrictions apply. However, be aware that FFTW works best if the number of grid points can be decomposed into as small as possible prime factors. Hence, integer powers of two are ideal but by no means mandatory. The only other restriction in WavePacket is that the number of points comprising the spatial grid has to be even so that inversion symmetry can be applied!