home *** CD-ROM | disk | FTP | other *** search
- /*
- ### procedure for computing fft ###
- */
- #include <stdio.h>
-
- void fft_compute(ftp)
- FILE *ftp;
- {
- int i,isign=1; /* should be made to be external later */
- double *fft_tmp;
- extern int tser_len;
- extern double *fft,*tser;
-
- for(i=0;i<tser_len;i++){
- fft[i*2]=tser[i];
- fft[i*2+1]=0;
- }
- fft_tmp = fft - 1;
- four1(fft_tmp,tser_len,isign);
- fft = fft_tmp + 1;
- }
-