More flexible fourier transforms

[aplavin]

Currently I don't see a way to compute the nonuniform fourier transform at arbitrary points. nufft2 is the closest variant, but it requires a regular grid of x, y coordinates in the transform space. A more flexible one should take any number of arbitrary (x, y) pairs and return the same number of fourier coefficients. Or did I overlook something?