Použití FFTW3 – C / C++ – Fórum – Programujte.com

Ahoj, snažím se napsat Matlabovskou funkci hilbert ()v C++ a narazil jsem na problém při použití fft (používám fftw3). Je to psáno i v dokumentaci, ale po vytvoření plánu a spuštění mi to vrací místo n prvků jen n/2+1, což já potřebuji n - tak jak to vrací matlabovská fce fft(). Je u knihovny FFTW3 nějaká možnost jak toho docílit?

Vstup je jedno-dimenzionální pole realných dat.

Jak to používám:

vector<double> * d = new vector<double>(data, data+size);
const int dsize = d->size();

fftw_plan p;
fftw_complex *out;

out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * d->size());

p = fftw_plan_dft_r2c_1d(2*(dsize/2+1), &(*d)[0], out, FFTW_ESTIMATE);
fftw_execute(p);    

#1 JakubD 

Nemalo by byt vystupne pole vecsie ako vstupne?

Skus nieco taketo:

 vector<double> in(data, data+size);
 vector<fftw_complex> out((in.size()/2+1)*2, 0);


 fftw_plan p;

 p = fftw_plan_dft_r2c_1d(in.size(), in.begin(), out.begin(), FFTW_ESTIMATE);
 fftw_execute(p);   

#2 vitamin
Bohužel, žádná změna. Nejde mi teda použít iterator na begin(), musel jsem tam dát adresu prvního prvku, ale to by snad nemělo mít vliv.