http://paulbourke.net/miscellaneous/dft/
Для вызова используется следующая функция:
Код
int FFT2D(COMPLEX **c,int nx,int ny,int dir)
Поискал в интернете примеры структур COMPLEX, применил такую:
Код
typedef struct cmplx
{
double real;
double imag;
} complex;
{
double real;
double imag;
} complex;
Создал одномерный масив из complex.
При передаче указателя в функцию компилятор выдает ошибку о несовместимости типов
**complex и *complex
Убрал одну звездочку выдает ошибку при передаче данных
Код
for (i=0;i<nx;i++)
{
real[i] = c[i][j].real;
imag[i] = c[i][j].imag;
}
{
real[i] = c[i][j].real;
imag[i] = c[i][j].imag;
}
Я заменил индекс [i][j] на [i*nx+j].
Хотелось бы применить правильное описание типа COMPLEX, чтобы не кроить исходный текст.