I have now completed the modifications to the FFT to handle squaring.
This wasn't hard at all.
Attached to this email is a diff of the changes for Brian. Of course
it won't apply automatically to mpir, but it should be possible to go
through and apply the same changes by hand.
Beware that in one
I have now written tuning code for the FFT.
The code is in the tune directory of the fft module in the program
tune_fft.c. It prints a valid fft_tuning.h to stdout.
It relies on the *nix timing functions, but presumably Windows has a
similar function. The actual timer doesn't matter too much. So
