I've started to use the FFT software you developed.  It works very well,
congratulations.  However, I did come across one bug when using arrays
of interleaved complex numbers: the rescaling was only being applied to
half of the array.  To fix this I made the following changes to

FFTN() in fftn.c

330c330
<       for (i = 0; i < iSign * nTotal; i += iSign)
---
>       for (i = 0; i < nTotal; i += iSign)

________________________________________________________________________
Chris Pudney              Biomedical Confocal Microscopy Research Centre
cpudney@alphapharm.pharm.uwa.edu.au           Department of Pharmacology
http://www.cs.uwa.edu.au/~chrisp         University of Western Australia
PH (618) 9346 4571  FX (618) 9346 3469      Nedlands, WA 6907, AUSTRALIA