[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Minor bug in fftw planner
From: |
Sukanta Basu |
Subject: |
Minor bug in fftw planner |
Date: |
Thu, 6 Jun 2013 14:16:05 -0400 |
Hi,
I believe I have detected a minor bug in fftw.
The first time (after starting octave), the following code is run, it
creates minor differences between f1 and f2. The differences can be
seen in the pcolor plots.
%-----------------------------
u = randn(128);
fftw('planner','estimate');
fftw('planner','estimate');
f1 = fft2(u);
fftw('planner','exhaustive');
fftw('planner','exhaustive');
f2 = fft2(u);
figure(1); pcolor(real(f1-f2)); colorbar;
figure(2); pcolor(imag(f1-f2)); colorbar;
%-----------------------------
However, if one repeats the exact same code for the second time (or
more), the differences disappear.
Best,
Sukanta
--
Sukanta Basu
Associate Professor
North Carolina State University
http://www4.ncsu.edu/~sbasu5/
- Minor bug in fftw planner,
Sukanta Basu <=