home *** CD-ROM | disk | FTP | other *** search
- #include <stdio.h>
-
- void fft_write(ftp)
- FILE *ftp;
- {
- int i;
- double power;
- extern int tser_len,graph_logy;
- extern double *fft;
-
- for(i=0;i<tser_len;i++){
- power = fft[2*i]*fft[2*i]+fft[2*i+1]*fft[2*i+1];
- if(graph_logy==1 && power==0)
- fprintf(ftp,"%d 1.e-10\n",i);
- else
- fprintf(ftp,"%d %.14lg\n",i,power);
- if(graph_logy==1)
- fprintf(ftp,"%d 1.e-10\n",i);
- else
- fprintf(ftp,"%d 0\n",i);
- }
- }
-
-