ooura fft库

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++){