Libtensor refused to build without fftw
sdf_sort_1d() now accepts more than one field to sort.
Comment the routines in fftw.h
Slight changes in the FFTW interface, plus unit testing of those functions
Minor changes in approx_eq() with more informative messages when testing.