ooura fft库的官网为:http://www.kurims.kyoto-u.ac.jp/~ooura/fft.html
DFT及其反变换
不用FFT库,用原始的方法实现DFT
/* test of fft*g.c */
#include <math.h>
#include <stdio.h>
#include <time.h>
#define PI 3.1415926
#define N (128)
typedef struct{
float real;
float image;
}complex;
double error_check(complex *a,complex *b,size_t n);
#define rdtscll(val) __asm__ __volatile__("rdtsc" : "=A" (val))
static void complex_dft(complex *src,complex *dest,size_t n){
for(size_t i = 0; i < n;i++){