请写一个长度为1024的FFT C语言函数代码

以下是一个使用 C 语言实现的快速傅里叶变换(FFT)的代码，输入的数组长度为 1024。

#include <stdio.h>
#include <math.h>
#define PI 3.14159265358979323846
 
void fft(double *x, double *y, int n) {
    for (int i = 0; i < n; i++) {
        y[i] = 0;
        for (int j = 0; j < n; j++) {
            y[i] += x[j] * cos(2 * PI * j * i