Skip to content
Snippets Groups Projects
Commit 6b4188a4 authored by Henrik (Grubba) Grubbström's avatar Henrik (Grubba) Grubbström
Browse files

Fixed a few bugs in Math.Transforms.FFT.rFFT() and Math.Transforms.FFT.rIFFT().

Rev: src/modules/Math/transforms.cmod:1.6
parent 0c15352b
No related branches found
No related tags found
No related merge requests found
...@@ -86,11 +86,12 @@ PIKECLASS FFT ...@@ -86,11 +86,12 @@ PIKECLASS FFT
THIS->r_out = malloc(sizeof(fftw_real)*n); THIS->r_out = malloc(sizeof(fftw_real)*n);
if (THIS->r_rc_plan!=n) if (THIS->r_rc_plan_size != n)
{ {
if (THIS->r_rc_plan>0) if (THIS->r_rc_plan_size > 0)
rfftw_destroy_plan(THIS->r_rc_plan); rfftw_destroy_plan(THIS->r_rc_plan);
THIS->r_rc_plan = rfftw_create_plan(n, FFTW_REAL_TO_COMPLEX, FFTW_ESTIMATE); THIS->r_rc_plan =
rfftw_create_plan(n, FFTW_REAL_TO_COMPLEX, FFTW_ESTIMATE);
THIS->r_rc_plan_size = n; THIS->r_rc_plan_size = n;
}; };
rfftw_one(THIS->r_rc_plan, THIS->r_in, THIS->r_out); rfftw_one(THIS->r_rc_plan, THIS->r_in, THIS->r_out);
...@@ -196,9 +197,10 @@ PIKECLASS FFT ...@@ -196,9 +197,10 @@ PIKECLASS FFT
if (THIS->r_cr_plan_size!=n) if (THIS->r_cr_plan_size!=n)
{ {
if (THIS->r_cr_plan>0) if (THIS->r_cr_plan_size > 0)
rfftw_destroy_plan(THIS->r_cr_plan); rfftw_destroy_plan(THIS->r_cr_plan);
THIS->r_cr_plan = rfftw_create_plan(n, FFTW_COMPLEX_TO_REAL, FFTW_ESTIMATE); THIS->r_cr_plan =
rfftw_create_plan(n, FFTW_COMPLEX_TO_REAL, FFTW_ESTIMATE);
THIS->r_cr_plan_size = n; THIS->r_cr_plan_size = n;
} }
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment