C code to process FFT in IGOR