fpc/packages/fftw/examples/example.pas
joost 66dfc61eaa * Placed fftw_getmem and fftw_freemem in interface section
* Removed dependency on libraries that are not directly called 
 * Fixed compilation of example (part of bug #12881)

git-svn-id: trunk@12505 -
2009-01-04 12:41:56 +00:00

32 lines
639 B
ObjectPascal

program example;
uses fftw_s;
const s=128;
var i,o:Pcomplex_single;
p:fftw_plan_single;
a:cardinal;
begin
fftw_getmem(i,s*sizeof(complex_single));
fftw_getmem(o,s*sizeof(complex_single));
p:=fftw_plan_dft_1d(128,i,o,fftw_forward,[fftw_estimate]);
for a:=0 to 127 do
begin
i[a].re:=(single(a)-64);
i[a].im:=0;
end;
writeln('input:');
for a:=0 to 127 do
writeln('(',i[a].re:8:4,',',i[a].im:8:4,')');
fftw_execute(p);
writeln('output:');
for a:=0 to 127 do
writeln('(',o[a].re:8:4,',',o[a].im:8:4,')');
fftw_destroy_plan(p);
fftw_freemem(i);
fftw_freemem(o);
end.