manual with words on fortran
FFT seems OK. To be tested by Xudong. Mathematica file to ease understanding
include choice of algorithm: bruteforce or fourierspace
manual
optimization enabled