friend pvector I(const pvector& in1)Return value:
- pvector I: Discrete (forward) Fourier Transform (DFT) of in1
Input:
- in1: pvector
Modify the Numerical Recipes FFT routine ``fourn'' by replacing all occurrences of ``float'' with ``double'', and rename the routine ``fournp480''. Then, uncomment the routine ``friend pvector I(..)'' in pvector.C, and debug with the ``Simple test of I'' in the main program.
Numerical Recipes Note: For debugging, keep in mind that Numerical Recipes requires that all FFT's be done with boxes with dimensions that are integer powers of two (1, 2, 4, 8, 16, etc.).
C++ Note: I found it easiest to compile fournp480 as a standard c code, thus it is declared as an extern ``C'' function at the top of pvector.C. You should change this if you decide to compile fournp480 as a C++ code.