1

次のスニペットを使用して、50枚の画像に対してバッチFFTを実行しようとしています。

pix3 = n*pix1*pix2;
fftwf_complex *in2, *f2h; //input for FFT2
in2 = (fftwf_complex*)fftwf_malloc(sizeof(fftwf_complex) * pix3);
f2h = (fftwf_complex*)fftwf_malloc(sizeof(fftwf_complex) * pix3);

for (i = 0; i < pix3; i++)
{
in2[i][0] = ref20[i]; //ref20 is an array of real values 
in2[i][1] = 0;
}    

for(i = 0; i < n; ++i)
{
  plan_forward2 = fftwf_plan_dft_2d (pix1, pix2, in2, f2h, FFTW_FORWARD, FFTW_ESTIMATE);
  fftwf_execute (plan_forward2);
  in2 += pix1*pix2;
  f2h += pix1*pix2;
}

ただし、最初の画像（つまり、最初のpix1 * pix2要素）のFFTしか取得できません。これを正しく行う方法について何か提案はありますか？

前もって感謝します

4

1 に答える 1

2

マニュアル 4.6 New-array Execute Functionsを注意深く読む必要があります。

ループの外で一度だけ計画を作成する必要があり、それを次のように使用します。

fftwf_execute_dft(plan_forward2, in2, f2h);
于 2013-07-31T03:32:41.653 に答える