Суть вопроса такова. Есть такой код:
#include "stdafx.h"
#define PI 3.1415926535897932384626433832795
#define FR 3
#define SAMPLES 1024
int _tmain(int argc, _TCHAR* argv[])
{
base signal[SAMPLES];
for (int i = 0; i < SAMPLES; ++i)
{
signal[i] = (sin(2*PI*i*FR), 0);
}
FFTCalculate(signal, SAMPLES, false);
FILE* f_pointer = fopen("fft_result.txt", "w");
for (int i = 0; i < SAMPLES; ++i)
{
char num = (char)signal[i].real();
int out = fprintf(f_pointer, &num);
}
int out = fclose(f_pointer);
return 0;
}
в нем вызывается функция FFTCalculate(signal, SAMPLES, false). Ниже приведен код функции:
void FFTCalculate(base signal[], int n, bool invert)
{
int log_N;
double x = frexp((double)n, &log_N);
calc_rev(n, log_N);
for(int i = 0; i < n; ++i)
if (i < rev[i])
swap(signal[i], signal[rev[i]]);
for(int len = 2; len <= n; len<<=1)
{
double ang = 2*PI/len * (invert?-1:1);
int len2 = len>>1;
base wlen (cos(ang), sin(ang));
wlen_pw[0] = base(1, 0);
for (int i=0; i<len2; ++i)
wlen_pw[i] = wlen_pw[i-1] * wlen;
for (int i=0; i<n; ++i)
{
base t;
base *pu = signal + i;
base *pv = signal + i + len2;
base *pu_end = signal + i + len2;
base *pw = wlen_pw;
for(; pu != pu_end; ++pu, ++pv, ++pu){
t = *pv * *pw;
*pv = *pu - t;
*pu += t;
}
}
}
if(invert)
for (int i=0; i<n; ++i)
signal[i] /= n;
}
При запуске программы появляется ошибка, связанная с выходом за пределы массива. Начал дебажить и выяснил, что когда вызывается функция FFTCalculate(signal, SAMPLES, false)
, то ей почему то передаются не signal, SAMPLES = 1024 и false
, а sugnal, 1245452 и true
см. скриншот ниже. В чем дело?