Improving interpolation accuracy

[lankef]

All iteration steps in the code use pseudospectral (FFT) derivatives and are relatively accurate. However, eval, some steps in coordinate conversion and axis calculation uses cubic spline interpolation, which can be much less accurate. Add the option to use some form of Fourier interpolation in eval, and the option to have high-resolution axis calculation and low-resolution iteration to improve speed.