added algorithms using FFTW, and giving up on direct inverse, I'll
smooth, then recompute